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Finishing the Job 


By L. A. HAWKINS 


General Electric Company, Schenectady, New York 


OM MON to all loyal citizens today must be 

the resolution to contribute, each to the 
limit of his ability, to the winning of this war 
which imperils our most cherished institutions. 
The greater the ability to contribute, the greater 
must be the inner compulsion to all-out effort, 
and since, as Dr. Conant has said, this is a 
physicist’s war, especial obligation rests on them 
to devote to the war problems their continuous 
intensive effort of brain and hand, laying aside 
for the duration all vocational and intellectual 
interests, however alluring, which would detract 
from the effectiveness or quantity of their war 
work. 

But such concentration should not bar a 
thoughtful response to such a presentation of 
vital post-war problems as is given in Dr. 
Bridgman’s ‘‘Challenge to Physicists’? in the 
April issue of this journal. Bitter experience has 
taught us that we may win a war but lose the 
peace. When this war is won, the work will be 
only half done. However intensive our present 
war effort, surely we should give thought to the 
consolidation and lasting protection of the 
temporary security that victory will bring. 

If we must look forward to a mere respite of 
barely two decades, such as the victory of 1918 
brought us, to be followed in endless cycle by 
new wars more terrible still, we might well 
despair of the future. But the physicist knows 
that past evils need not be forever endured, and 


that new conceptions, adequately acted upon. 
may completely change the course of events. 

If a serious defect develops in the operation 
of a complicated apparatus, and if a repair is 
made and the apparatus continues to function 
for a time, only to develop similar but more 
serious trouble, it requires no special brilliance 
for a physicist to deduce that the repair was 
inadequate, that it did not reach the seat of the 
trouble, or that it was insufficient in its execution. 

So, too, the physicist, looking back to 1918, 
and the decisive victory then attained at the 
cost of so much blood and treasure, and reflecting 
on how short a time after that victory the world 
remained ‘‘safe for democracy,” is forced to 
conclude that winning the war was not enough, 
that it was only an emergency repair job which 
did not reach the real cause of the trouble. 
He must conclude that we left the job only half 
done, when, after the armistice, we yielded to 
war weariness, discarded all international re- 
sponsibilities and commitments, and shut our 
eyes to the sores which were left festering in 
Europe and which invited all sorts of malignant 
growths, while we reveled in a post-war business 
boom. 

He sees that victory alone will not safeguard 
and perpetuate peace and security. He will 
agree with Dr. Bridgman that new conceptions 
and a new morale are needed, and that it is his 
duty to try to make others see the need. 
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Geological Applications of Nuclear Physics 


By CLARK GOODMAN 


Massachusetts Institute of Technology, Cambridge, Massachusetts 


INTRODUCTION 


ULLY a decade before nuclear phenomena 

were recognized as such, radioactivity was 
improving our knowledge of the earth’s heat. 
Since this beginning, nuclear physics has aided 
to some extent in understanding each of the 
three major branches of geological study, i.e., 
the history, the structure, and the composition 
of the earth. In the discussion which follows, 
these applications are briefly reviewed and pos- 
sible future developments are outlined. It should 
be emphasized at the outset, however, that the 
benefits have not been unidirectional. Nuclear 
physicists have acquired new knowledge of their 
own field through collaborative researches with 
geologists. 


RADIOACTIVITY OF TERRESTRIAL MATERIALS 
Distribution 


Most geological applications of nuclear physics 
depend either directly or indirectly on the radio- 
active properties of certain naturally occurring 
atomic species. The physical properties, esti- 
mated abundance and modes of occurrence of the 
geologically more important of these elements are 
given in Table I. The only commercially im- 
portant ones in this group are potassium, radium, 
thorium, and to a lesser extent, uranium. Except 
for potassium, each of these elements averages 
less than 0.001 percent in the crust of the earth. 
However, because of the energetic radiations 
which they emit, the radioactive elements can be 
studied with considerably greater accuracy than 
many of the more abundant elements. 

The primordial source of all rocks appears to 
be a substratum of basaltic* material lying a few 
tens of kilometers below the surface of the earth. 
The radioactive elements presumably are dis- 
tributed fairly uniformly throughout this layer 
but at concentrations substantially lower than 


* The schematic classification of igneous rocks given in 
Fig. 1 is included for the benefit of non-geologists. 
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are found in acidic surface rocks. This decrease 
of radioactive content with depth and increasing 
basicity is the result of differentiation from the 
source material. Upon cooling of any molten 
mass, crystallization begins and chemical separa- 
tion takes place in accordance with the existing 
pressure, temperature, and composition condi- 
tions.” Certain mineral species, such as allanite, 
zircon, and sphene, seem to have crystal struc- 
tures admitting a greater number of the large 
uranium and thorium atoms.’ Accordingly, if 
conditions are favorable to the formation of these 
minerals, which are usually only minor con- 
stituents in igneous rocks, the resulting rock may 
show localized concentrations of the radioactive 
elements. On the other hand, the more common 
rock-forming minerals, such as quartz, feldspar, 
pyroxene, amphibole, micas, and magnetite, are 
less tolerant of foreign radioactive atoms. During 
crystallization only relatively few uranium or 
thorium atoms are admitted, and these are 
found randomly distributed throughout the struc- 
tures of the crystals.t 

The rejected radioactive atoms, together with 
other incompatible, stable atoms, thus become 
concentrated in the uncrystallized portion of 
magmatic solutions. In some instances a rela- 
tively small residuum may contain the bulk of 
the uranium and thorium originally distributed 
throughout a volume several orders of magnitude 
larger. These concentrated solutions may cause 
replacement of other minerals or may crystallize 
directly into minerals such as uraninite, thorite, 
thorianite, and monazite. Such minerals contain 
the parent elements as primary constituents 
($0.1 percent) and are generally classified as 
“radioactive minerals.’ Because of these genetic 
relationships, radioactive minerals are often 
found in coarsely crystalline intrusive rock bodies 
known as pegmatites. 





+ Piggot (reference 4) observed a localized concentration 
of radium on the surface of micas. However, too few meas- 
urements of this effect have been made to draw any con- 
clusions as to its generality. 


JOURNAL OF APPLIED PHYSICS 








tion 
\eas- 
con- 


ICS 


Properties 


Regardless of their chemical combination or 
physical environment in the earth, the radio- 
active elements spontaneously disintegrate at 
known, systematic rates. All experimental at- 
tempts to alter this disintegration by changes of 
pressure and temperature have failed. On the 
basis of present concepts of the nature of atomic 
nuclei, no effect is to be expected. However, 
geological phenomena furnish the most convin- 
cing evidence of the constancy of radioactive 
decay during the lifetime of the earth. The 
Geiger-Nuttall law relates the range of a-rays 
with the disintegration constant \ of the emitting 


nucleus. Alpha-rays produce a pleochroic dis- 
coloration in certain minerals, micas and fluorite 
in particular. This effect results in a micro- 
scopically visible sphere, or halo in the central 
section, surrounding minute inclusions of a-active 
material, usually a tiny crystal of zircon or 
allanite. In numerous direct photometric meas- 
urements of pleochroic haloes,’ many of which 
have required geological periods of time for their 
formation, no variation in range, and hence no 
variation in \, have been observed for any radio- 
active element. 

After an a-ray has lost its energy by absorp- 
tion, it settles down as a normal, well-behaved 


TABLE I. Physical properties, abundance and occurrence of the geologically important radioactive elements. 


Active Radia- Energy Decay constant 


Element isotope tions in Mev A in sec.~! 
Potassium Kk? B 1.2 1.4+0.3 x 107!" 
K?0 oY 2.0 4x 10719 
Rubidium Rb%? B 0.13 3.541 X107'9 
Samarium Sm" a 2.59 2.510719 
Radon Rn a 5.49 2.10 10-6 
Tn a 6.28 1.27107 
An a 6.82 0.177 
Radium Ra a 4.87 1.4 107" 
ThxX a 5.78 2.20 x 10-6 
Thorium Th a 4.31 1.58 x 10718 
RdTh a 5.47 1.16 107 
lo a 4.66 3X 107% 
Protactinium Pa a S.1 7.0 107% 
Uranium Us a 4.13 4.87 x 1078 
U2 a 4.4 3.140.1 107!" 





tMean abundance of 
active isotope 


Half-period 
T g/g rock 


Occurrence 





2.2+0.4 x 10-6 





1.6+0.3 x 10%y Rocks, evaporites, 


5X 10!y and sea-water. 
6.342 X 10!y ca 6X 1076 Lepidolite and 
other Li-bearing 
minerals. 
0.9 10"y ca 1076 Samarskite, cerite 
and monazite. 
3.82d ca 3X 1078 Genetically associ- 
54.5s ca 7X 10-2 ated with uranium, 
3.92s ca 2x 107-*4 thorium, and ac- 


tino-uranium, re- 
spectively, in rocks, 
radioactive miner- 
als, atmosphere, 
mineral waters, and 
sea water. 


1600y 0.5+0.1 x 10-" 


Genetically associ- 
3.64d ca 4X 1078 


ated with uranium 
and thorium, re- 
spectively, also 
found separately in 
certain mineral 
springs. 


1.40 « 10!y 5.6+1xX10-* Rocks, monazite, 
1.9y ca 8X 107'6 thorite, and thoria- 
8X 10*y ca 3X10 nite. Ionium gen- 
etically associated 

with uranium. 


3.1 104y ca 4X 107% Genetically associ- 
ated with actino- 
uranium (U**). 
4.52 10%y 1.4+0.3 x 1076 Rocks, —_uraninite 
7.140.3 x 108y 1.0 10-8 pitchblende and 


carnotite. 





t Abundances obtained as follows: K*° from average chemical analyses of intermediate igneous rocks, see Daly (reference 1) and Grout (reference 
12), and isotopic abundance by Nier (reference 13), Rb’? and Sm"8 from Hevesy (reference 14) with isotopic abundances by Nier (reference 13) 


and Dempster (reference 15), respectively; Th and Ra from Evans and Goodman's (reference 7) averages; others calculated from known genetic re- 
lations using indicated decay constants. 
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helium atom, presumably lodged by self-enforced 
admittance in the structure of the mineral. The 
heat produced by absorption of the a-rays con- 


stitutes the major source of terrestrial radio- 
active heat.* After the entire series of trans- 
formations has taken place, the residual products 
of the three main radioactive families are the 
stable isotopes of lead Pb***, Pb*°7, and Pb?®*. 
The genetic relations for the series are shown 
diagrammatically in Fig. 2. 

Even though the parent elements are dis- 
tributed as individual atoms in a rock or mineral, 
the number within only one milligram is more 
than sufficient to allow application of statistical 
averages. It requires about one million years to 
establish radioactive equilibrium in the UI series, 
and a much shorter time suffices for the AcU and 
Th series. Once this condition obtains, the rate 
of production of heat, helium, and radiogenic 
lead by the entire series is known from a deter- 


* The 6-rays and the A-rays, which have much greater 
ranges of penetration than a-rays and which cause much 
less distortion of the crystal structure, contribute about 11 
percent of the heat (reference 7). 





mination of the concentration of any one member 
or from direct a-count measurements. In practice 
the chemically inert, gaseous emanations, Rn, 
An, and Tn, are usually the most convenient 
individual members to separate and measure. 
Because the electrical effects produced by the 
minute contents of radioactive elements in ter- 
restrial materials are so small, a number of 
analytical precautions must be taken to insure 
reliable results. While satisfactory techniques 
now appear to be established,” * these have 
evolved only after four decades of research by 
numerous investigators. During this time hun- 
dreds of measurements were made, most of 
which are now of qualitative value only. Because 
of these inaccuracies, a number of the earlier 
interpretations have undergone rather complete 
revision.*!° The National Research Council's 
Committee on Standards of Radioactivity in 
collaboration with the National Bureau of Stand- 
ards recently has made available" dilute radium 
standards and standard rock samples which 
should aid materially in avoiding errors in future 
terrestrial radioactivity measurements. 
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Fic. 1. Classification of igneous rocks. Igneous rocks are classified on the basis of chemical composition into four major 
groups. The proportions of acidic and basic constituents (geologically speaking) are designated by 2A and LB, respec- 
tively. The center of each rectangular area coincides with the world-wide averages (reference 1) for the rock types in- 
cluded. The horizontal and vertical spread has been chosen to include at least 90 percent of the rocks of each group with 
some overlapping to indicate the lack of a sharp boundary. The basaltic substratum in the earth is considered to have an 
average composition corresponding to about the center of the basic group. 
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Fic. 2. The radioactive series. Both UI and AcU are iso- 
topes of uranium and thus occur in a definite proportion 
in all uranium-bearing minerals. Because AcU decays more 
rapidly than UI, the AcU/UI ratio has decreased systemat- 
ically during geological time. The diagram includes the 
recently discovered (reference 6) branch product Ack. 


GEOCHRONOLOGY 
General Principles 


Radioactivity provides the only quantitative 
method of measuring geological time. This appli- 
cation aids the geologist in correlating igneous 
rocks, in measuring the span of paleontological 
and diastrophic changes, and in determining the 
rate of erosion and deposition. In some instances 
such knowledge is of considerable economic value. 

If the age of a rock or mineral is defined as the 
time ¢ that has elapsed since final solidification 
or crystallization, it can be expressed as the 
ratio: 


Age= Amount of disintegration product/ 
Rate of production of product (1) 


in which the numerator represents the amount of 
radiogenic lead, Pb*, or helium, He, formed and 
the denominator measures the rate at which the 
disintegration product is being formed in the 
material. Because the parent elements decay so 
slowly, this rate is essentially constant for most 
terrestrial samples and is directly proportional to 
the concentration of these elements. 

In the application of Eq. (1) to rocks and 
minerals, three basic requirements must be ful- 
filled. These have been discussed in detail else- 
where® and are summarized for the present 
purpose in Fig. 3. The age methods and the 
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materials to which they are applicable fall into 
two groups as the result of these requirements. 
In general, the lead method is applicable only 
to radioactive minerals and the helium method 
only to selected rocks and minerals. The reason 
for this division is apparent from Fig. 4. 


Helium Method 


Only recently has sufficient attention been 
given to the requirements outlined in Fig. 3. 
Accordingly, most of the helium ages reported in 
the literature are now known to be unreliable.* !8 
The major sources of this difficulty have been the 
analytical techniques and the recently demon- 
strated leakage of helium from certain common 
rock minerals.'® Reliable measurements on a 
routine basis now appear possible. However, the 
retentivity of minerals for helium has been found 
to be extremely variable and is of paramount 
importance in the application of this method. 

Although it has not been definitely established 
that magnetite retains all its helium, there are a 
number of indications which suggest this possi- 
bility. From tests on several common mineral 
types, magnetite has been found to possess the 
highest specific retentivity. Helium age ratios for 
a series of magnetites ranging in geologic age 
from pre-Cambrian to mid-Tertiary with few 
exceptions show the proper sequence.'? Within 
the small experimental uncertainty of the meas- 
urements the same helium age ratio was obtained 
for a series of Triassic magnetites from widely 
separated localities and containing a broad range 


1. Raa 
2. At t 


Rate 


0, He=Pb=0 
3.From t=0 to t=t: 


= PB 
= Th 


He 
Urs 


4 

















Fic. 3. Requirements for age measurements by radioac- 
tivity. The observed radioactivity, Raa, must accurately 
represent the rate of production of He or Pb* in the ma- 
terial. Hence, requirement (1) implies both reliable meas- 
urement and representative sampling. (2) At the time of 
solidification or crystallization, t=0, the rock or mineral 
contained no helium or radiogenic lead. (3) No addition or 
subtraction of these disintegration products or their sources 
has occurred during the geological history of the material, 
i.e., from ¢=0 to the present time, t=¢. 
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of helium and radioactive contents.'® In this 
connection it is particularly pertinent that the 
proportions of radium and thorium in these 
specimens varied considerably. 

Accordingly, as a tentative referent, magnetite 
has been adopted for helium age studies, and 
considerable progress has been made already in 
the establishment of a geologic time scale based 
on measurements of this mineral. The results to 
date are summarized in Fig. 5 and are compared 
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Fic. 4. Representative radioactive content and decay prod- 
ucts of terrestrial materials. The amounts of radiogenic lead 
(Pb*) and helium (He) generated in 100 million years by 
the representative uranium (U) and thorium (Th) contents 
and the amounts of ordinary lead (Pb) to be expected are 
indicated. The lead content of such a radioactive mineral 
can be separated quantitatively by the usual chemical 
techniques and measured gravimetrically. The lead content 
of ordinary rocks is considerably beyond the range of 
chemical methods. Similarly the uranium and thorium 
content of radioactive minerals can be determined by 
chemical means, but the content of these elements in or- 
dinary rocks requires sensitive physical detection methods 
for measurement. The helium content, corresponding to 
chemically separable quantities of radiogenic lead, is so 
great that some loss is inevitable. Accordingly, helium 
ratios on radioactive minerals can represent only minimal 
age values and cannot be compared directly with corre- 
sponding lead ratios. It is conceivable that a material may 
be found with an intermediate composition, such as the 
hypothetical rock, that would be suitable for both the lead 
and helium age methods. Thus far a direct comparison of 
this type has not been made. 
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with the best available lead age ratios on radio- 
active minerals. 


Lead Method 


In the case of radioactive minerals, the denom- 
inator of Eq. (1) is based on chemical analyses 
of U and Th together with the decay constants 
of these elements.* While the lead can also be 
separated and determined chemically, it is neces- 
sary to identify the proportions of ordinary lead, 
Pb, and radiogenic lead, Pb*, since only the 
latter may be used as the numerator of the 
age equation. Three methods have been used: 
1. atomic weight determination, 2. mass spectro- 
graphic determination of the abundance of the 
lead isotopes, 204, 206, 207 and 208, and 3. a de- 
termination similar to (2) based on the hyperfine 
structure of the lead spectrum. The mass spectro- 
graphic method yields the highest accuracy and 
offers a number of other advantages.f 

Assuming reliable analytical results, what are 
the limitations imposed by the fundamental re- 
quirements discussed above as applied to the lead 
method? If any alteration, such as the leaching 
out of the parent elements, U and Th, has 
occurred, the calculated rate of production of 
lead will not be representative of the rate which 
existed during the lifetime of the mineral. Be- 
cause the internal radiations shatter radioactive 
minerals, these materials are particularly sus- 
ceptible to leaching and replacement by hydro- 
thermaf agencies. The extent to which alteration 
affects the validity of the measurements has not 
been thoroughly established. Some minerals with 


.considerable evidence of alteration show con- 


sistent age measurements, while others with 
much less apparent change yield obviously un- 
reliable ratios. Accordingly, the practice has been 
to select minerals, whenever possible, which show 
the least alteration, but to base the final accept- 
ance of the results on their compatibility with 
other lead age measurements rather than on any 
absolute criterion. This procedure is an unfortu- 
nate one, but appears to be the only method 


* The researches of Kovarik and Adams (reference 19) 
are largely responsible for the reduction of the uncertainty 
in the knowledge of Av and Arn. 

+ Following Aston’s initial investigations (reference 20) 
Nier (reference 21) has made a systematic study of a 
number of leads separated from radioactive minerals by 
Baxter and others. 
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Fic. 5. Lead and helium time scales. While the sequence of geologic periods is well-established from paleontologic evi- 
dence, radioactive methods provide the only quantitative estimate of their length. The time scale shown is based largely 
on lead ratio ages, see Eq. (1). Support for the essential correctness of this distribution is provided by the helium ages on 
numerous magnetites and a few sulfides shown on the right (references 16 and 17). The observed age for the upper group 
of these minerals is indicated immediately below on the upper (post Proterozoic) scale of abscissas, while the correspond- 
ing geologic assignments are shown on the scale of ordinates. It is seen that with the exception of the Blagodat specimen, 
which is a questionable Devonian, the agreement with the lead ratio time scale is quite good. The isotope ratio, based on 
the systematic variation of AcD/RaG with time, in some instances is quite different from the lead ratio on the same radio- 
active mineral (reference 21). More consistent agreement between these two independent lead methods is observed for 
older minerals (reference 21). These data have not been included because of the geologic uncertainty in comparing helium 


ages in the earlier Pre-Cambrian. 


available at present. A similar procedure is 
necessary in the study of magnetites. However, 
the mutual compatibility of these helium age 
measurements can be determined more reliably 
than similar comparisons for radioactive minerals, 
because there are far more magnetite samples 
available and their geologic sequence is more 
clearly established. 


GEOPHYSICS 
Well Logging 


Schlumberger and associates”? have estab- 
lished practical methods of measuring the elec- 
trical characteristics of subsurface formations. 
The two most useful electrical parameters are the 
impedance or resistivity and the natural poten- 
tial. From empirical study, it has been found that 
oil bearing sand often can be distinguished from 
salt-water bearing sand, and the common sedi- 
mentary rocks, such as limestone, shale, and 
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sandstone, generally can be identified without 
expensive subsurface sampling. While the elec- 
trical methods fill an important need in petroleum 
geology, they cannot be used on cased wells. 
Some nuclear radiations, on the other hand, can 
penetrate casing.* 

For many years ionization chambers and 
quantum counters have been used to determine 
the y-ray activity of terrestrial materials, but the 
extension of these techniques to well logging is 
a recent development. Spak (1937)* first reported 
this application using a Geiger-Miiller counter. 
He credits Gorkov with making the original 
suggestion some three years earlier. Independ- 
ently, Howell and his collaborators (1935)* 
applied high pressure ionization chambers for 
drill hole y-ray measurements. These earlier 
techniques have been ingeniously improved by 


* Essentially all completed wells contain one or more 
liners of jointed steel pipe called casing. 
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Howell and Frosch*® and by Scherbatskoy, 
Westby, and Fearon.” *? Duplicable y-logs, both 
of cased and uncased wells, can now be recorded 
continuously at logging speeds up to thirty feet 
per minute. Examples of the types of records 
obtained are given in Fig. 6. 

Because the radioactive content of sedimentary 
rocks varies roughly with type—for example pure 
sandstones and limestones are low while shales are 
relatively high—there is a correlation between 
electrical logs and y-logs. Thus this newer tech- 
nique inherits considerable geological information 
that has accrued from the extensive application 
of electrical methods. Some idea of the degree of 
correlation between the methods can be obtained 
from Fig. 7. It is not to be implied that y-logs 
replace electrical logs. Instead the methods are 
supplementary and provide the geophysicist with 
a number of different parameters from which to 
deduce the stratigraphy of subsurface formations. 

In addition to this more general application, 
y-logs afford a means of locating cement behind 
oil well casing that is not subject to the limita- 
tions of the older temperature method. It is 
customary to cement the casing when reducing 
the size of the drill hole or when shutting off 
undesirable water flowage from the penetrated 
formation. Cement is pumped down and flows 
around the casing and into the surrounding for- 
mations for a generally unknown distance and 
height. Howell and Frosch** have shown that the 
inclusion of only about 0.5 percent carnotite, or 
any other equally radioactive mineral, in the 
cement allows the position and relative thickness 
of the cement behind the casing to be determined 
with considerable accuracy. Examples of this 
application are shown in Fig. 8. 


Earth Heat 
As Slichter?® has recently shown, the observed 
surface heat flows establish definite upper limits 
to the total heat generation which may exist 


within about 100 or 200 kilometers of the surface : 


of the earth. On the basis of early measurements, 
it appeared that the rate of generation of heat 
by the radioactive elements was so great that a 
stable crust could hardly exist. Accordingly, 
volcanism, diastrophism, and other dynamic 
processes were considered by some investigators 
to be manifestations of the transfer of this 
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Fic. 6. Electrical and gamma-ray logs. The upper traces 
indicate impedance and potential measurements on a por- 
tion of an uncased well. The lower traces are y-logs in the 
same well made before and after casing. (After Howell.) 


excess energy to the surface where it could be 
radiated into space. A rapid decrease in radio- 
active content with depth was also postulated; 
an effect supported to some extent by experi- 
mental evidence but not to a sufficient degree to 
be entirely convincing. 

As indicated in the preceding section of age 
measurements, many of the earlier terrestrial 
radioactivity measurements are of qualitative 
value only and are subject to a significant down- 
ward revision. The evidence which led to this 
change has been presented in detail elsewhere.” ' 
The immediate geophysical importance of this 
revision is that it allows the existence of a stable 
crust of a thickness indicated by seismic studies 
and leaves sufficient additional radioactive ma- 
terial for the underlying mantle and core not to 
be entirely devoid of these elements. 

Because of the large thermal inertia of the 
earth and the only approximate knowledge of the 
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heat flux to the surface, thermal conditions 
below about 200 kilometers are shrouded in 
mystery. In fact it appears impossible without 
the aid of as yet unknown devices to determine 
whether the earth is heating or cooling at depth. 


Exploration 


Numerous attempts have been made to corre- 
late surface observations of radioactivity with 
subsurface features. These measurements have 
been of two types: (1) determinations of varia- 
tions in y-ray intensity, and (2) surveys of the 
radioactivity of soil gases. 

A few feet of overburden suffice to reduce the 
y-ray intensity from underlying rocks to an 
undetectable level. Hence, observed variations on 
the surface are likely the result of local changes 
in the radioactivity of the soil rather than the 
result of structural changes below. Some success 
has been reported” in the location of radioactive 
mineral deposits by y-ray measurements. How- 
ever, the detector must be quite close to the ore 
body before a significant increase is observable. 

The free radon content of soil gases is likewise 
a function of the features of the soil rather than 
of the underlying rocks. For example, Clark and 
Botset *° observed a rather close correlation 
between the radon and the heavy mineral content 
of soil. Unless a subsurface fault exists in ma- 
terials of widely different radioactive content, no 
significant change in the radon content of the 
overlying soil is to be expected. In general few 
geologically informative data have accrued from 
surface radioactivity observations although meas- 
urements of this type have been of considerable 
value in studies of the cosmic radiation. 


GEOCHEMISTRY 
General Statement 


The less abundant elements can often be used 
as geochemical tracers of the reactions which 
have taken place in the formation of terrestrial 
materials. Quantitative spectroscopic methods 
are finding increased application in the study of 
the less abundant stable elements. The radio- 
active elements can be measured when present 
in spectroscopically and chemically undetectable 
amounts. Applications of these two physical 
techniques have already added considerably to 
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the knowledge of geochemistry. Future researches 
undoubtedly will greatly extend this knowledge. 


Oceanography 


Although oceanographic radioactive relations 
have been studied for over three decades, only 
within the last few years have sufficiently accu- 
rate measurements been made to allow reliable 
conclusions to be drawn. An attempt has been 
made to summarize the more useful data in the 
schematic diagram of Fig. 9. The disintegration 
of continental rocks by erosion transfers the 
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Fic. 7. Correlation between electrical logs and gamma-logs. 
The gamma-logs were made while lowering the detection 
apparatus and again while removing it. The two directions 
are seen to give substantially the same results within the 
statistical fluctuations. Numerous correlations are observed 
between the y-logs and the electrical logs, major ones being 
at 2700-2800, 3125-3175, 3400, and 3800 feet. (After 
Howell.) 
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Fic. 8. Gamma-logs reveal radioactive cement. The shaded 


area is proportional to the quantity of cement behind the 


casing in each of the three wells. The marked differences observed serve to indicate the efficacy of the cementing operations. 


(After Howell and Frosch.) 


radioactive elements from the land to the sea. 
The clastic materials are deposited near shore 
with little change in composition. However, the 
elements dissolved out of the rocks or in colloidal 
suspension are subject to precipitation according 
to their individual chemical properties and the 
nature of the environment. 

From the measurements it is inferred that 
thorium and ionium (being isotopic) are precipi- 
tated, together with manganese and ferric iron,* 
in the region of slightly basic conditions near 
shore. Since uranium is not precipitated under 
these conditions, unusually high Th/U ratios 
should obtain in these sediments. Such marine 
basins are considered the most likely regions for 
the formation of source beds of petroleum. The 
evidence, therefore, suggests that petroleum may 
originate in rocks with high Th/U ratios. No 
direct verification of this hypothesis has appeared 
though there is indirect supporting evidence from 
the low Th/U ratio observed for sea water at a 
distance from shore. It is hoped that research 
now in progress will cast some light on this 
interesting problem. 


Mineralization 
As we have seen, alpha-radiation produces dis- 
coloration in certain minerals, biotite and fluorite 


* The limiting concentration of ferric ions is about 
10~*° 0 /oo. 
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in particular. Because the range of natural a-rays 
is only a few tens of micra in these minerals, 
this pleochroic effect is restricted to a micro- 
scopic region surrounding the a-emitter. In the 
course of extensive studies of these ‘‘pleochroic 
haloes,’’ Henderson and co-workers? discovered a 
number of unusual types, presumably formed by 
the radiations from radioactive elements in solu- 
tion moving through veinlets in the minerals. 
Many of these haloes are extinct, i.e., they no 
longer contain active material. However, a per- 
manent, photographic-like record remains to 
provide information concerning movement of 
radioactive solutions that occurred millions of 
years ago. The conclusions to be drawn from 
these vestiges of the past are not unexpected from 
related knowledge, but such haloes provide 
unique, graphic evidence of the microscopic 
processes of solution movement and deserve more 
detailed geologic consideration. 

Beginning with Miigge’s*® original autoradio- 
graphs of radioactive minerals, the photographic 
action of radiations from radioactive substances 
has been used widely to determine the distribu- 
tion of uranium and thorium in terrestrial ma- 
terials. The locus of the active atoms and qualita- 
tive estimates of their concentrations can be 
obtained quite readily by this method. This 
information is of interest in understanding the 
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modes of formation of radioactive minerals and 
the subsequent alteration that may have oc- 
curred. The latter is of value in determining the 
suitability of the mineral for lead age studies 
(see page 280). More recently the photographic 
technique has been extended to ores of ordinarily 
stable elements.*7 ** 

Because of their lack of charge, neutrons are 
able to penetrate the densest matter and coalesce 
with the nuclei of practically all elements. If a 
given element has a sufficiently large capture 
cross section for neutrons and in addition yields 
a radioactive product, it is possible to determine 
the distribution and concentration of the element 
from the radiations emitted. Samples of rock or 
ore are bombarded by stray, slow neutrons from 
the cyclotron, while other research is in progress, 
then removed and placed in direct contact with 
a photographic film. The autoradiograph pro- 
duced by the 8-radiation reveals the locus and 
relative concentration of the active element.*? 
When correlated with petrographic study, new 
information is obtained concerning the genesis of 
minerals. Manganese, gold, tungsten, copper, 
arsenic, sodium, potassium, and phosphorus are 
particularly well suited to this type of investiga- 
tion.*® Examples of the autoradiographs obtained 
are given in Fig. 10. Since many of the elements 
have characteristic rates of decay, the identity 
of the active atoms can often be obtained from 
counter measurements of the material. When 
combined. with similar measurements on known 
concentrations, the method can be used as a 
semi-quantitative assay of the ore. The more 
common elements, such as silicon, oxygen, alumi- 
num, and sulfur, are activated relatively little 
and hence do not interfere with the application 
of these techniques. 

In the lead age method the abundance of lead 
isotopes varies with the proportion of uranium 
and thorium in the mineral, the age of the 
mineral and the amount of ordinary lead present. 
Until recently it was assumed that all ordinary 
lead was composed of a fixed proportion of 
isotopes.* Nier*® has shown that there is a small 
but definite variation in the relative abundance 
of the lead isotopes that may be of considerable 


* This conclusion was based on the apparent constant 
atomic weight of ordinary leads and the lack of variations 
in other heavy elements. Specific definitions of the various 
forms of lead are included in reference 8, pp. 531-532. 
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significance in tracing the origin of primary lead 
ores. Thus far these observations have not re- 
ceived adequate geologic interprétation. 

Still more uncertain is the anomalous occur- 
rence of excess helium in certain minerals, beryls 
and sylvines in particular. Fay, Gluckauf, and 
Paneth*® have demonstrated that this helium is 
not attributable to nuclear disintegration by 
y-rays of the beryllium jn the beryls. Were it not 
for the apparent correlation with age,*! the excess 
helium could be considered as a product of 
magmatic differentiation. It is conceivable that 
the helium in these minerals contains an un- 
usually high proportion of the recently dis- 
covered® stable isotope, He*. As yet no isotopic 
abundance measurements have been made. This 
tantalizing problem remains for future solution. 


Isotopic Variations 


The variation in abundance of lead isotopes in 
ordinary lead has been considered in the pre- 
ceding section. Similar variations also have been 
observed in a number of the lighter elements. 

Following the discovery of deuterium (H?) 
numerous determinations of the density of waters 
from various sources were made. Variations in 
density were commonly observed, some being as 
high as 20 parts per million, i.e., y=20. These 
differences were attributed to supposed varia- 
tions in the H!/H? ratio, presumably due to 
differences in the rates of diffusion in nature. 
However, Dole and his associates** have demon- 
strated conclusively that the observed differences 
are largely the result of faulty laboratory purifi- 
cation and preparation of the water. After elimi- 
nating uncertainties from this source, only minor 
differences in density are observed and these are 
due primarily to isotopic variations in oxygen 
rather than in hydrogen. Geologically the sig- 
nificance of such variations, therefore, is much 
less than was previously supposed. 

The abundance of He* atoms in atmospheric 
helium is only about one per million He‘. For 
this reason discovery of the lighter isotope was 
delayed until application of the cyclotron as a 
high intensity mass spectrometer.” Since He* 
must be lost from the earth’s atmosphere more 
readily than He‘, terrestrial helium would be 
expected to have a higher concentration of the 
lighter isotope than atmospheric helium. While 
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only a few measurements have been made, the 
opposite appears to be the case. Helium from 


natural gases contains only about one-tenth as 
much He’ as atmospheric helium. A likely expla- 
nation is that terrestrial helium shows a greater 
dilution by radiogenic He* of the primordial 
helium, which presumably contained He* as one 
of its components in a somewhat higher propor- 
tion than is now found in helium from either 
source. An alternate possibility is that He’ is 


formed in the atmosphere by an unknown nuclear 
reaction involving cosmic radiation. 

Nier and his associates*® have established vari- 
ations up to 5 percent in the C?/C™ abundance 
ratio of carbon obtained from several natural 
sources of widely different geologic age and geo- 
graphic origin. No trends with age were found. 
Plant spores, algae, petroleum, bituminous shales, 
torbanite, and kerosene shales have relatively 
higher ratios than limestones, zeolitic calcite, 











| 
ry ~ | 
+ 5 Ar + o 
Radioactive Relalions in the Ocean 
. | 
* “ « A | 
i) 
a S RAL | 
m 10°" core R . | 
= SEA LEVEL 46 610 | 
a 3 - 
tele, ° FZ COMPOSITION 
Goa Fe tn, 10) “es “T= GS N PARTS PER LITER 
U } 2 
tom fi ‘ “ ~ 7 
es ee y ALGAE . NaCci- 2723 
x) LUOAC 10 GRa/G 
e270 ; R) - MgCl, --- 3807 
ee \ DEPTH PTH (E-K- M,P-R) e', : 
eee € s . MgSO. 658 
@ oe \ fon mw ' aM aes . 
€,2,%0 J : CaSO, -- 1260 
fee Vo K a Kz $04------ 0863 
0° ' >> mo ‘ 2W« 
e2e3 — j = 4s cacao - 0123 
“oe % E-« ‘ ——_ ‘a » 012 
LPN \ } ee MgB, --- 0076 
fo. the dy PLANKTON 10 “ora/G 
* om, A ‘ (E-K-M,P-R) } 
| 0% Ra G . ) 
| 2o ue | 
} ‘ 4 
| to A 
_ ‘ “ . | 
~ N | 
/ \ | 
CONTINENTAL ROCKS ¢ > ( 
Acidic Ba eormentary an a y * Ly ee 
e — 7? 
Un 6/6 30 | ’ Veale 
’ SHELL 3*10°“oRa/oCa CaCO, 
Ra 4 ‘4 (P-R) 
h ° 3 ‘ 3 
5 4 
~ — j 
2 WLOBIGERINA OOZE = ——" —— RED CLAY —_~——*'t 
= — eer 
Ra WW WO 'G/G Ra IN 10 Go | 
a oe #+-$9 2. | 
30 2 T 
| « 60 | 
PTH DEPTH ) | 
- P-u) “ ° } 
M " M 
| | 
200 + 
eens = — — wateaiesiintth eats 


Fic. 9. Radioactive relations in the ocean. Considerable selecticn of data on the basis of evident reliability 





has been 


made. Individual investigations are designated as follows: (E-G) Evans and Goodman (reference 7), (E-K-M) Evans, Kip, 
and Moberg (reference 31), (F-R) Foyn and Rona (reference 32), (K) Karlik (reference 32), (M) Mauchly (reference 33), 
(P-R) Pettersson and Rona (reference 32), and (P-U) Piggot and Urry (reference 34). The composition of sea water is 
from the Challenger expedition as reported by Posnjak (reference 35). (E-K-M) observed a decrease in excess radon, i.e., 
above equilibrium amount, with depth and an increase in radium content as shown in the upper left graph. (P-R) did 
not confirm the latter effect, although their observed average radium content, 0.7 X 10~'* g/g is in good agreement with 
that of (E-K-M). A slight increase in uranium content with depth was observed by (K) as shown in the adjacent graph. 
U/Ra=2.84 10° at equilibrium. The available evidence indicates the lack of radioactive equilibrium in sea water and 
that radium is being precipitated in excess of uranium. The lower graphs represent the interesting observations of (P-U) 


: on bottom core samples collected by the ingenious mechanism developed by Piggot. These results clearly show the lack 


of radioactive equilibrium near the top of these sediments. Two extremes are observed; the usual case very likely is a 
combination of these two. Excess precipitation of ionium (and hence likewise its isotope thorium) relative to radium and 
uranium occurs in globigerina oozes, while excess radium relative to ionium and of the latter relative to uranium occurs in 
red clays. The variation of the content of these elements with depth in the cores follows closely the decay curves. Thus a 
constant rate of deposition is indicated, and a means of estimating this rate likewise is made available. From the above 
measurements, it is inferred that thorium and ionium are precipitated (together with manganese and ferric iron) in the 
region of slightly basic conditions near shore. Since uranium is not precipitated under similar circumstances, unusually 
high Th/U ratios (>>4) should characterize these sediments. Algae and plankton possess the ability to concentrate 
radium more than 100-fold above the concentration of the water in which they grow. Conversely, calcareous shells select 
calcium in preference to radium over 1000-fold. These interesting biochemical mechanisms await explanation. The 
emanation content of air over the ocean apparently is due to the exhalation of radon from sea water, which has a Henry’s 


law coefficient, a, of 0.16. 
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marine shell, and sea water. Since the C?/C® 
ratio would be expected to remain substantially 
the same in a given life cycle, these observations 
are of value in unravelling the genesis of certain 
terrestrial materials. The data support independ- 
ent evidence for the geologic origins of petroleums 
and limestones. 

Early mass spectrographic measurements*® 
were able to show that there were no large scale 
variations in the abundance of oxygen isotopes in 
various natural substances. However, these ob- 
servations were not capable of detecting a differ- 
ence of less than about 2.5 percent in the O!8/O'* 
ratio, this ratio being very close to 1/500. Be- 
cause the oxygen in carbonate rocks and iron 
oxide ores is easily recovered in the form of 
water and density measurements on the water 
can be made with high precision, Dole and 
Slobod*? were able to study the isotopic abun- 
dance of oxygen in these materials from widely 
different geologic horizons. No significant differ- 
ence was observed between the oldest rocks 
(Grenville marble) and the youngest (clam shells), 
from which it is concluded tentatively that 
“there has been no change in the isotopic com- 
position of the oxygen in the oceans (from which 
these rocks were precipitated) over a_ period 
greater than a billion years.’’ This is in agree- 
ment with Nier and Gulbransen’s observations*® 
on the lack of any significant variation in the 
carbon isotopes from similar sources. The oxygen 
in certain cuprite ores®® and in iron ores of Kee- 
watin and Middle Huronian ages was found to 
have the same isotopic composition as in normal 
water, which is significantly different from at- 
mospheric oxygen. From this it is concluded that 
these ores were probably formed by precipitation 
from water, substantiating the general opinion of 
geologists. The excessive density of the water 
prepared from fossil type Silurian ores is attrib- 
uted to the presence of carbonates which are 
shown to have slightly higher proportions of the 
heavier isotope. 

Potassium is made up of three isotopes of 
atomic weights 39, 40, and 41 in relative abun- 
dances of about 8300, 1, and 585, respectively. 
K*° is 8-radioactive* with a half-period of 1.6 
_* Potassium is also y-radioactive with approximately 3 
ys per 100 6’s. The y-rays presumably result from K- 
electron capture by K* see reference 8, p. 538 for discussion. 
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Fic. 10. Autoradiographs of ores. The distribution of 
tungsten in the Spanish wolframite ore (upper left) is 
indicated by the relative brightness in the print prepared 
from the autoradiograph, as shown immediately below. 
The 8-rays from W!*? (7 = 24 hr.) are largely responsible for 
the photographic action. Similarly the 8-rays from Au! 
(T =2.7 days) reveal the distribution of gold in a sample of 
Cripple Creek ore activated by slow neutrons. The bright- 
est areas in the photograph (lower right) correspond to 
native gold, the gray areas to calaverite (gold telluride). 


The dark regions in both lower prints correspond to quartz 
which is activated relatively little. 





+0.3X10° yr. The other isotopes are stable. 
Brewer*® has studied the isotopic ratio, K**/K*, 
in various minerals, in plants, and in animal 
tissues. In general no significant variation has 
been observed for mineral potassium. On the 
other hand, it is claimed that plant tissues 
exhibit appreciable variations and usually have 
the abundance ratio of the soil in which they 
grow, while kelp and other seaweeds show a 
preference for the heavier isotope. Because of the 
widespread distribution of potassium and its 
biologic as well as geologic significance, it would 
be of interest to have these interesting observa- 
tions substantiated by further studies. 

Schumb, Evans, and Leaders*® have shown 
that potassium separated from the Pultusk 
meteorite has the same K*°/K ratio as terrestrial 
potassium within an experimental uncertainty of 
about 3 percent. This research was undertaken 
in an attempt to use potassium as an index of 
the age of meteorites. 


FUTURE APPLICATIONS 


An extension of the helium method of age 
determination to minerals other than magnetite 
already has begun. This should allow a broader 
application of the method particularly if sedi- 


287 





mentary minerals prove reliable. The separation 
and measurement of rock magnetite and other 
minerals will enable a more complete time scale 
to be established after which this physical tool 
can be used to solve many stratigraphic problems. 
Future extension of the lead method should be 
largely based on Pb*°7/Pb?°* ratios. More direct 
comparisons between the helium and lead meth- 
ods are needed since the ultimate success of both 
depends upon their mutual agreement. 

A new era in petroleum prospecting should 
result from an understanding of the origin of oil. 


Nuclear physics may well be of aid in solving this ' 


complex problem. Unless the present emergency 
makes it impossible, a systematic study of the 
genesis of petroleum will begin this year under 
the sponsorship of the American 
Institute. 


Petroleum 


The calibration of y-ray logs would make 
available a wealth of geochemical and geophysical 
data on sedimentary rocks.* While terrestrial 
radioactivity measurements appear incapable of 
casting further light on the thermal state of the 
earth at depths, many important local and near 
surface problems depend upon this information. 
Hence, further geothermal radioactivity measure- 
ments would be of value, particularly if these 
were correlated directly with diffusivity and 
geothermal gradient measurements. The future 
should also see an extension of the application 
of neutrons in geophysics. Already the relative 
scattering coefficients of neutrons appear to be 
useful in well-logging.*® The results of this type 
of survey reveal the relative hydrogen content 
(largely in the form of water) in the formations. 
By extending the neutron induced activity 
method to rocks, an indication of their mineral 
and fluid contents could be obtained. This 
application awaits the development of compact, 
intense neutron sources. 

Further knowledge of the variations in the 
abundance of isotopes in nature and of the 
distribution of radioactive elements in rocks and 
minerals should lead to a better understanding 
of the fundamental processes of geochemistry. 


* Added in proof: A beginning in this direction has just 
been reported by Pontecorvo, Geophys. 7, 90-94 (1942), 
but the sensitivity of his instrument is only sufficient to 
measure sediments containing >5X10-" g Ra/g or its 
equivalent. 
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There is considerable evidence that unknown 
radioactive elements, possibly an entire new 
series,*' may be present in highly differentiated 
terrestrial materials. The origin of helium- 
bearing gases has received only speculative 
explanation. This problem deserves further 
attention, particularly in view of the recently 
discovered light isotope of helium. 

Now that agreement among different investi- 
gators has been attained in oceanographic radio- 
activity, whole new fields have been opened for 
research. The rate of deposition of deep sea 
sediments can be determined from the degree of 
radioactive equilibrium between the geologically 
short-lived elements and their parents. There is 
every reason to suppose that a similar technique 
would be applicable to near shore and continental 
sediments. Artificial radioactive tracers are being 
used and medicine. 
Certain organisms are known to concentrate 


extensively in biology 
radium in preference to its chemical homologs. 
Further study of the mechanism of this and 
related processes should be of value not only in 
biology but in understanding the role of these 
organisms in the formation of sedimentary rocks. 
The radioactivity of potassium affords a means 
of studying its importance in the precipitation 
and consolidation of clays. 

All in all, it seems certain that in the future 
the newest developments of physics, especially 
those exemplified by the under- 
standing of the atomic nucleus, will find ever 
increasing application in the earth sciences. 


progress in 
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Frictional Phenomena. IX 


By ANDREW GEMANT 


The Detroit Edison Company, Detroit, Michigan 


Chapter IX. Applications of Liquid Viscosity 
to Electrical Insulating Liquids, 
Particularly in High 
Voltage Cables 


Abstract 


HIS chapter deals with: (1) mechanical 

processes involving viscosity, particularly 
(a) impregnation, (b) breathing, (c) oil migration 
with respect to high voltage cables; (2) chemical 
processes involving viscosity, such as (a) Voltol 
process, (b) wax formation in cables, (c) poly- 
merization of styrene ; (3) electrical conductivity. 
With regard to the latter, the fundamental 
quantities, besides viscosity, that determine the 
conductivity are discussed, together with the 
various experimental methods whose interpreta- 
tion allows the quantities in question to be 
computed. 


36. Mechanical Processes and Viscosity 


There are three fields in which liquid viscosity 
has to be considered in connection with liquid 
dielectrics. One of these is that of mechanical 
processes, particularly dynamic ones, concerned 
with motions of the liquid. Another is that of 
chemical processes, often accompanied by a 
change of viscosity. A third is that of electrical 
conductivity, controlled among other factors by 
the viscosity of the medium. All these are dealt 
with briefly in this chapter. 

Let us begin with the motion of insulating 
liquids. There are several electrical devices in 
which some kind of liquid motion takes place, 
but according to the nature of this monograph, 
aimed more at a logical and constructive presen- 
tation than at completeness, only one device, 
i.e. high voltage cables, will be considered as an 
illustration. There are several different phases in 
manufacture and operation that are of interest 
in this respect. 

(a) Impregnation.—The insulation of high volt- 
age cables consists of paper strips wound around 
the central copper conductor, the body of paper 
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being afterwards impregnated with oil. The most 
usual impregnation process is carried out by 
winding the cable (copper conductor plus paper) 
around a drum, placing the assembly in a tank 
and, as a first step, evacuating the whole system. 
After that, oil is allowed to enter the tank. 
During this phase of the process oil will pene- 
trate, at least to a large extent, into the pores 
of the paper chiefly by means of capillary forces. 
The paper will be partially impregnated, the 
rest of the pores being filled with gas under 
reduced pressure attained during this phase. 
Finally hydraulic pressure will drive the oil com- 
pletely into the pores and other capillary spaces 
left between the paper strips. Only a minute 
amount (less than 0.5 percent) of gas space 
under atmospheric pressure (so-called voids) will 
be left over. 

The speed of the impregnation process is that 
of a viscous flow. By Eq. (13) a resistance value 
R was defined for any sort of channel, through 
which flow takes place under any driving force. 
For a single circular capillary tube it was 
R=8n/r*. The resistance R is the product of the 
viscosity of the medium (here the oil) and a 
geometric factor, depending upon the configura- 
tion of the channel or channel-system. The 
value of R for papers is determined by apparatus 
originally introduced by Emanueli! and similar to 
that described in connection with sound absorb- 
ing materials in Section 17 (c). The geometric 
factor mentioned is independent of the nature 
of the medium, as long as no change of the 
structure, such as swelling of fibers, takes place. 
The resistance of the paper insulation (Fig. 63) 
is less than that of a single layer times the number 
of layers, because a certain leakage path exists 
along the spaces (so-called butt spaces) left 
between neighboring turns of the same layer and 
along the narrow space left between neighboring 
turns of subsequent layers (see the dotted line in 
Fig. 63). 

The lower the viscosity of the oil, therefore, 
and the more porous the paper, the faster will be 
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the impregnation. However, in the usual cable 
design fairly viscous oils (about 20 poise at room 
temperature) have to be used for several reasons, 
partly discussed below under (c) Migration. The 
porosity, on the other hand, is required to be 
low, since the electrical strength increases with de- 
creasing porosity. These features, then, slow down 
impregnation considerably, for which reason the 
temperature is raised to around 115°C, and the 
viscosity lowered to the order of 0.3 poise. On 
the basis of this low viscosity level, impregnation 
is completed in about 48 hours. It should be 
recalled once more that the driving force during 
the first part of impregnation is chiefly the 


























Fic. 63. Diagram of oil flow through paper insulation. 


capillary force, the external pressure acting only 
in the second phase of the process. The applica- 
bility of the equations of viscous flow is, however, 
independent of the nature of the driving force. 

In recent years a large number of new cable 
constructions have been suggested and adopted. 
In contradistinction to these new types the usual 
one is called, although not very correctly, “‘solid’’- 
type cable, since its insulation, paper completely 
impregnated with a viscous oil, behaves in 
several respects like a solid. One of the new 
types, the Beaver-type cable,? consists of paper 
tapes that are preimpregnated with a very 
viscous oil and thus requires no impregnation 
after taping. The interstices left between the 
individual tapes are then filled with nitrogen gas 
under a pressure of 200 lb. per sq. in. This type 
is mentioned not because of its manifold ad- 
vantages, which do not interest us here, but 
because it does away with the impregnation of 
the whole cable which, as indicated above, is a 
fairly difficult process, particularly if the degree 
of impregnation is required to be as high as 
possible. The impregnation of single paper tapes 
can be completed much faster and more com- 
pletely, owing to the low value of the resist- 
ance R. 

(b) Breathing.—The viscosity of the impreg- 
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nating oil or compound (so-called if mixed with 
resin) has an essential influence not only on the 
impregnation of the cable, but also on its opera- 
tion. During operation the cable undergoes daily 
load cycles, both the amperage and temperature 
varying in a cyclic manner. During a temperature 
rise the oil expands, during a drop it contracts, 
and this motion is called breathing. 

This cyclic process in itself is reversible, thus 
harmless. Unfortunately, another irreversible, 
thus harmful, process follows in its trail. The 
sheath of the cable, made of lead because its 
flexibility and plasticity are desirable in the 
manufacture, handling, and installation of the 
cable, follows the expansion of the oil by a 
plastic expansion, but does not follow the con- 
traction. The reason is that during expansion the 
pressure difference inside and outside the sheath 
is higher than the yield value of lead (see Chapter 
XII) ; thus the lead will flow. During contraction 
of the oil the lead sheath does not follow because 
the pressure difference inside and outside, usually 
around 10, or at the most 15 pounds, is lower 
than the yield value. Hence the contraction of 
the oil in the space inside the sheath that has at 
this point a practically constant volume (the 
expansion coefficient of the lead being around ten 
times smaller than that of the oil) will initiate a 
displacement of the oil through the paper, re- 
sulting in the formation of voids, or aft least 
leaving behind some poorly impregnated layers 
of paper. Because of the ensuing gaseous ioniza- 
tion in these voids these occurrences are highly 
undesirable. 

The actual process during this cooling period 
depends largely upon the viscosity of the oil. 
Cooling starts from the outside, thus a con- 
siderable temperature gradient exists in the radial 
direction. Because of the large temperature 
coefficient of the viscosity of heavy oils and com- 
pounds used in “‘solid’”’ type cables, the outer 
part of the insulation will have a higher viscosity 
than the inner parts near the copper conductor. 

This viscosity gradient will determine to a 
certain extent the location of the displacement 
just mentioned, in that flow through the papers 
within the inner half of the insulation is doubtless 
greatly facilitated as compared with that in the 
outer parts. In agreement with this conclusion 
it is found that the majority of voids and dry 
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spots are located in the inner half of the insula- 
tion. This is undesirable because the field in- 
tensity near the conductor is substantially higher 
than near the sheath, and thus the probability 
of ionization and subsequent deterioration is 
high. 

There are several designs to counteract this 
unwanted occurrence. A recent one is that by 
W. A. Del Mar,’ utilizing two different types of 
paper, a more porous one for the inner layers, 
and a denser one, possessing a higher modulus of 
elasticity, for the outer layers. The compressive 
action of the latter prevents to a certain extent 
void formation in the neighborhood of the 
conductor. 

A more radical departure from the ‘‘solid”’ 
cable is the compression cable,*® the principle of 
which is that the whole cable is placed in a 
steel or a second lead pipe, the interspace con- 
taining gas of 200 lb. per sq. in. This pressure 
makes the lead sheath collapse during the cooling 
period, so that the mechanism of void formation 
as determined by the viscous properties of the oil 
cannot become effective. Usually an oval-shaped 
cross section (Fig. 64) is an added feature of the 
compression cable. When the internal pressure 
increases because of heating, the sheath will not 
expand uniformly as a circular lead sheath would 





Fic. 64. Breathing of a compression type cable. 


do, but will undergo a deformation, assuming a 
shape of lesser eccentricity (dotted line), thus 
accommodating a larger volume of oil. Upon 
cooling, the external pressure will flatten the 
cross section to its original shape, which thus 
appears to hinge on four points like P in the 
figure. 

(c) Migration—A third way in which vis- 
cosity plays a dominating role is in connection 
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with the position of cables, as in hilly country or 
at terminals, where a gradient of elevation exists. 
In such cases a pressure head acts along the 
cable and a flow of oil will take place longitudi- 
nally (not radially as in the previous cases). 
Again the rate of migration will be determined by 


A 








D E 


Fic. 65. Diagram illustrating oil migration in cables, 
according to Scott. 


the viscosity, and thus, since migration is evi- 
dently undesirable, compounds of high viscosity 
at operating temperatures (up to 70°C) are 
selected as impregnants. 

T. R. Scott® among others carried out con- 
siderable work in this field, and the effects of 
migration can best be illustrated in Fig. 65. 
It shows diagrammatically a cable running verti- 
cally at sections A—B and C—D, and horizontally 
at B-C and D-E. The cable is terminated at 
both ends by so-called potheads that contain 
either gas cushions or bellows, allowing a free 
inflow and outflow of compound. The pressure 
heads A—B and C-D cause the oil to migrate 
longitudinally from left to right, involving a 
decrease of pressure in the pothead A and an 
increase at E. 

The flow along B—C is determined by pressure 
head A-B and the resistance of the length B—C; 
similarly the flow in D-E is determined by the 
pressure head C—D and the resistance D—E. If we 
assume that B—C is much longer than D-—E, then 
the resistances of the respective branches are in 
the same ratio and the respective flows are in an 
inverse ratio. The result of this difference, the 
two pressure heads being equal, is an excess of 
the outflow at C over the inflow at the same point 
—in other words, a formation of voids. 

Scott also carried out actual measurements of 
the value of R per centimeter of cable. The flow 
takes place along the conductor and below the 


JOURNAL OF APPLIED PHYSICS 








ire 
C; 
he 
we 
en 
in 
an 


of 


int 





TABLE X. Flow resistance of cable conductors in longi- 
tudinal direction, according to Scott. 





Cross section of stranded 
conductor, sq. cm 


Flow resistance, R per cm 
cable in c.g.s. units 


1.0 35,000 
iS 20,000 
1.6 12,000 
2.0 8,000 


sheath. Values of R for stranded conductors of 
varying cross section are shown in Table X. 

An effective method of preventing longitudinal 
viscous flow of the compound in cables is the 
application of so-called stop joints. It is common 
practice to limit individual lengths of cable to 
about 300 feet and join these lengths together. 
It is possible to use a solid or plastic type of 
insulating material in the joints that will prevent 
a free flow of the compound through the joint. 
Thus migration is restricted to individual sec- 
tions, and cannot assume high values. 

Scott suggests the use of polymeric styrene in 
stop joints. Paper strips are previously impreg- 
nated with the monomer, then polymerized and 
separated. Such papers are then used in pre- 
paring the joint in a cold process, so-called since 
no subsequent polymerization by heat is re- 
quired. The single styrenated paper of 5-mil 
thickness has a breakdown strength of 1500 volts 
per mil. Such tapes are now produced on a com- 
mercial scale. 


37. Chemical Processes and Viscosity 


A second field relating to insulating liquids in 
which viscosity plays an essential role is that of 
certain chemical processes the course of which 
can be followed by noting the increase in vis- 
cosity. The chemical processes are condensation 
and polymerization, the latter having been dis- 
cussed already in Section 30. Both of these 
processes are accompanied by an increase of the 
molecular weight, the former with and the latter 
without an evolution of hydrogen gas. 

Reverting to Section 21, it was shown that 
the constant B of Eq. (68) increases with the 
molar heat of vaporization, thus with the molecu- 
lar weight. It is essentially the same B constant 
that at ordinary temperature determines the 
variation of 7 in a given series of materials, as 
can be seen from Fig. 31, where the height of 
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the lines increases with their slope. Thus it can 
be concluded that the viscosity of a substance will 
increase with an increasing degree of polymeriza- 
tion. This relation, by the way, should not be 
confused with the similar relation of Staudinger 
(Section 30) concerned with colloidal solutions of 
macromolecular substances. In the present case 
we deal with homogeneous liquids, although it 
cannot be denied that in certain instances both 
mechanisms might act together. 

Several interesting examples can be quoted to 
illustrate the significance of the above-mentioned 
relation. 

(a) In the so-called ‘“‘Voltol process’? an in- 
dustrial application of the phenomenon was 
attempted. Lubricating oils (see Chapter X) 
must have a certain minimum viscosity and the 
Voltol process is aimed at producing such oils 
from low viscosity materials by means of an 
electric discharge in hydrogen gas at reduced 
pressure. Thus the molecular weight of mineral 
oils could be raised from 400 to 650. Viscosity of 
around 0.15 poise could be raised to about 3. 
The processes involved seem to be both con- 
densation and polymerization. 

(b) In high voltage cables a similar process 
takes place, wherever voids are present between 
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F1G. 66. Viscosity vs. temperature of hydrocarbons (after 
Sticher and Thomas). (a) Before bombardment; (b) after 
bombardment. (1) Synthetic oil; (2) synthetic turpentine 
polymer; (3) n-hexadecene. 
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layers of impregnated paper. In such voids 
gaseous ionization takes place, and this discharge 
process polymerizes the oil. The result is the 
formation of highly viscous and finally solid 
wax-like products. The present writer has de- 
scribed several methods* for the study of this 


process. A recent interesting study was published 
by Sticher and Thomas,’ attempting a quantita- 
tive evaluation of the chemical changes. Also 
they extend the research to compounds of known 
chemical constitution. From among the many 
results of this research one particular item of 
interest is presented here. Figure 66 shows the 
viscosity as a function of temperature of three 
particular liquids (a) before bombardment (6d) 
after bombardment at 60 cycles for four hours. 
The liquids are specified in the caption to the 
figure. It can be seen that about 50 percent in- 
crease in viscosity has occurred during four 
hours of treatment. In contradiction to the 
Voltol process the changes in high voltage cables 
are not welcome, since they are signs of a 
progressive deterioration of the insulation. Ac- 
cording to Sticher and coworkers®™ and _ to 
Shanklin,'® however, wax formation might lead 
to self-healing and stabilization in certain cases, 
as in the low gas-pressure cable, operating at 
around 15 lb. per sq. in. gauge pressure. In both 
instances quoted, viscosity increases as the 
molecular weight increases, and is a good 
measure of the progress of the chemical reaction. 

(c) A third interesting example is polystyrene, 
of which mention was made in Section 36 in 
connection with Scott’s work on terminals and 
joints. The polymerization of monostyrene has 
become a process of considerable technical im- 
portance, and thus it is of interest that its 
progress can be followed by viscosity measure- 
ments as in the foregoing cases. 

Polymerization of the monomer occurs fairly 
easily, in fact, much too easily for shipping and 
storing of the material, during which time a 
polymerization is unwanted. It was found that 
certain added materials act as inhibitors. The 
monomer then can be stored without any diffi- 
culty. It is of interest to find out how such a 
monomer polymerizes if it is heated to the proper 
temperature. This comparison can be made 
easily by following the changes in viscosity with 
time. Such data of Scott are shown in Fig. 67 for 
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Fic. 67. Viscosity of styrene vs. time of polymerization, 
after Scott. (1) Pure styrene; (2) styrene+0.05 percent 
quinone. 


a heating temperature of 120°C. While curve 1 
refers to the pure monomer, curve 2 refers to 
styrene to which was added 0.05 percent of 
quinone as inhibitor. It can be seen that an 
induction period of 100 minutes elapses before 
the process starts, but after that the rate of the 
reaction is the same as for the pure material. 
The change of viscosity in this case is consider- 
able, 1 : 100 in 100 minutes. 

If too much inhibitor has been added, it is 
necessary to remove it by means of a distillation 
process before polymerization will start. Addition 
of only small amounts is, therefore, necessary if 
the optimum results are to be obtained. 


38. Electrical Conductivity and Viscosity 


The third field relating to insulating liquids in 
connection with viscous processes is their elec- 
trical conductivity. It is a long established fact 
that the term “‘insulating’’ does not mean zero 
conduction of current; it serves only to designate 
a range of conductivity sufficiently low for the 
currents not to cause any serious disturbance 
from the standpoint of the electrical apparatus 
in which the liquid is used. The small current 
across the liquid in these cases is an unwanted 
occurrence only if it increases with time, since 
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it then might reach a level causing interference 
with the chief function of the appliance. Thus the 
concept “‘insulating’’ is purely relative as it 
depends upon the appliance in question. In a 
rough way, however, it is useful to consider 
liquids of a conductivity of an order less than 
10-'° mho/cm as “‘insulating.”’ 

The discussion to be presented in this section 
will be only sketchy, since there are many other 
detailed and complete treatments on the subject. 
First the fundamental equation controlling ionic 
mobility will be discussed, and it will be shown 
how the latter depends upon the viscosity of the 
liquid. The viscosity is, however, only one factor 
in the electrical conductivity since the latter 
depends also upon the concentration, charge, and 
size of the ions. Considering the importance of 
the question, a little digression along these lines 
will be made to point out the experimental 
methods for determining all fundamental factors 
involved. 

In order to fix our ideas, the general case will 
be considered in which at least one of the two 
ions present (cations and anions) are large 
colloidal aggregates. This is indicated in Fig. 68 
for the case of large anions, assumed to contain 
five (generally z) elementary units e of electricity. 
The cations can be large, too, or small as shown 


“ey 


zr, 


Fic. 68. Diagram of colloidal ion. 


in Fig. 68, carrying only one unit each. The 
number of such composite units per cm’ should 
be n. The radii of the two types of ions are rq 
and r., respectively. 

The conductivity o is the total charge flowing 
through 1 cm? during 1 sec. if unit field strength 
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is applied. The charge involved in each unit 
shown in Fig. 68 is ze; thus the charge per cm’ 
is nze. We now define ionic mobilities m, and m, 
such as to give the absolute velocities in unit 
field strength. The two ions having opposite 
signs, their opposite velocities have to be added 
if the total charge flowing across 1 cm? is to be 
computed. Thus a prism of base 1 and height 
(m-+m,) will be displaced across 1 cm? during 
1 sec. The charge contained in this prism is 
nze(m.+m,); thus 


o=nze(m.+mz,). (104) 


It now remains to determine the mobilities. In 
Section 26 Stokes’ equation was discussed (Eq. 
(78)), as giving the viscous resistance of a sphere 
in a liquid. Although it is strictly valid only for 
spheres much larger than the molecules of the 
liquid, Einstein suggested its use for ions as well, 
and this conception proved to be quite fruitful. 
The electric force acting on the anion, for in- 
stance, is ze if the field strength is unity; thus 
the condition for equilibrium will be: 


2e= Orne at; 
hence, the mobility : 
Ma = 2e/Ornf a. (105) 


A similar equation with z=1 will be valid for 
the cations of Fig. 68. In these equations 7 is 
the viscosity of the medium. 

Equations (104) and (105) give the electrical 
conductivity as a function of the viscosity. of 
the liquid and the number, charge, and size of 
the ions. 

The correctness of Eq. (105) as to the in- 
fluence of » has been proved by several means. 
One of the proofs is the so-called Walden rule, 
stating that the mobility of ions is inversely 
proportional to the viscosity of the solvent, a 
rule found empirically before the theoretical 
Eq. (105) was established. Table XI, after data 
of Ulich," gives the mobilities m and the products 
mn for a few ions in different organic solvents. 
While the constancy of the product mn for the 
ion N(C2Hs5)4* is quite remarkable, the agree- 
ment is less close for the other ions, but sufficient 
to show the validity in respect to n of Eq. (105). 

Another way of checking the viscosity relation 
is by varying the temperature. In this regard the 
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TABLE XI. Ionic mobilities at 25°C after Ulich. 


Solvent 1007 N(CeHs)4* N(CHs)a* Na* J 
Methy! 0.545 m 54 65 46 61 
alcohol mn 0.294 0.354 0.251 0.332 
Ethy! 1.087 m 27.2 29.2 21.9 26.3 
alcohol min 0.295 0.317 0.238 0.285 
Acetone 0.316 m 93.1 102.8 68.6 115.9 
mn 0.294 0.325 0.217 0.366 
Furfural 1.494 m 19.6 12.7 28.4 
mn 0,293 0.190 0.424 
Aceto 0.346 m 85 81 109 
nitrile mn 0.294 0.28 0.38 
Epichlor- 1.03 m 28.6 32.8 35.4 
hydrine mn 0.295 0.338 0.365 
Pyridine 0.890 m 35 25 47 
mn 0.294 0.22 0.42 


fact has to be borne in mind that the temperature 
affects not only the viscosity n but also the other 
characteristic quantities, n, r, and z; thus the 
final effect upon o might be very complicated 
indeed. However, in certain simple systems it 
could be proved without doubt that o varies 
with the temperature inversely to , the product 
on being constant. The validity of this relation 
is known from experience with aqueous solutions 
of totally dissociated electrolytes, and has been 
established experimentally for insulating oils by 
J. B. Whitehead” and co-workers. 

On the basis of this relation, several conclu- 
sions can be drawn. One is of a general and 
practical nature, simply stating that in selecting 
insulating oils for cables or other high voltage 
equipment, oils of higher viscosity should be 
given preference, from the standpoint of con- 
ductivity, since, other properties being constant, 
the conductivity will decrease with increasing 
viscosity. Of course, other factors often will 
decide the choice in a contrary direction, for 
instance, mechanical considerations, as shown in 
Section 36, or chemical stability both with regard 
to electronic bombardment and oxidation. The 
chemical stability does not depend directly upon 
viscosity, but since more viscous oils have a 
higher molecular weight, and since this, to a 
certain degree, is an indication of lesser stability, 
there is sometimes an inverse correlation between 
viscosity and chemical stability. 

Another consequence of the o(7) relation is of a 
quantitative nature. If measurements of both o 
and » are taken, for a series of oil samples or 
for a series of temperatures, it is possible to carry 
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out an analysis of the data by considering the 
product (on). In doing this the viscosity effect 
is, so to speak, eliminated and only the other 
effects upon o are left. If, for instance, it is 
assumed that r and z remain constant, the 
change of the product (on) will tell us directly 
how n, the concentration of the particles, changes 
owing to varying electrolytic dissociation. If, on 
the other hand, there is additional information 
to the effect that electrolytic dissociation is 
complete, thus constant, then the variation of 
(on) is a measure of the changes occurring in z 
and r because of ionic association. In the general 
case both the degree of electrolytic dissociation 
and the degree of ionic association will change. 
The reason is that an ionic association (formation 
of complex colloidal ions) reduces the concentra- 
tion of single ions, and thus shifts the electrolytic 
equilibrium 


undissociated molecule@single ions 
in the direction of the upper arrow. 
Such a case can be seen in Fig. 69, referring 
to the equivalent conductance of tetraisoamy!l- 
ammonium nitrate in dioxane-water mixtures at 
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Fic. 69. Equivalent conductance of tetraisoamylam- 
monium nitrate vs. concentration, after Kraus and Fuoss. 
Solvent: water-dioxane mixtures: (1) pure water; (2) 85 
percent dioxane; (3) 96 percent dioxane; (4) pure dioxane. 
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25°C, obtained by Kraus and Fuoss." The 
equivalent conductance A, defined by 


A=10000/c, (106) 


where ¢=concentration in equivalents per liter, 
is useful in comparing solutions of different con- 
centration, since it gives the conductance per 
unit concentration, and thus depends on the 
secondary changes in 1, r, and z, as caused by a 
change in the total concentration. Thus, for a 
totally dissociated electrolyte, A should be inde- 
pendent of the concentration (disregarding the 
so-called Debye effect), whereas for weak elec- 
trolytes, having a finite dissociation constant, A 
should decrease with increasing concentration. 

Figure 69 indicates several curves referring to 
increasing percentage of dioxane and decreasing 
dielectric constant. Whereas A for water as 
solvent is constant, it decreases for a mixture 
containing 85 percent dioxane. Thus it appears 
that the electrolyte is strong in water and weak 
in dioxane. However, the curve for 96 percent 
dioxane goes through a minimum, and that for 
100 percent dioxane increases practically through 
its whole course. The explanation for this in- 
crease, according to Bjerrum,’ as well as to 
Kraus and Fuoss, is the formation of triple ions 
by means of an association of single ions and 
undissociated molecules. This process shifts the 
electrolytic equilibrium towards increased dis- 
sociation, thus increasing the total number of ions. 

The process illustrated by this case is probably 
very general in insulating liquids (characterized 
by low dielectric constants) and an analysis of 
data on a, in which the first step is the building 
of the product (en), will be helpful in discovering 
similar occurrences. The plots of the product 
(on) versus temperature in the quoted work of 
Whitehead often show deviations from a hori- 
zontal line for reasons similar to those just 
discussed. 

Balsbaugh and co-workers,'® conducting very 
extensive researches on deterioration of oils, also 
analyze their recent data, first reducing the con- 
ductivity values to the product (7). 

Caution is required in interpreting o vs. tem- 
perature curves on oils when extended to the 
very low temperatures at which the oil gradually 
solidifies. The oil first becomes a gel, as de- 
scribed in Chapter VIII, in separating into a 
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Fic. 70. Conductivity vs. time obtained from transients 
on two particular oils after J. B. Whitehead. O Experi- 
mental points. Full lines computed after an equation of 
Gemant. 


liquid and a crystalline phase. The electrical 
conductivity of such a heterogeneous system is 
then a sum of the two components, to be calcu- 
lated with due regard to the geometric shape of 
the two phases. In such instances unusual o-tem- 
perature curves are obtained, as shown and 
analyzed by the writer'® who found that the 
conductivity of oils passes through a minimum 
with decreasing temperature. Similar curves were 
obtained by W. Jackson for paraffin wax. 
Several investigators prefer the measurement 
of o by means of alternating voltages instead of 
d.c. voltages. In such cases specific a.c. losses, 
which generally show more or less sharp peaks at 
certain critical frequencies, are superimposed 
upon the d.c. conductivity. If a.c. data are ob- 
tained, it is useful to attempt a separation of 
d.c. and specific a.c. losses, before conclusions 
regarding the nature of ions are drawn. Such 
separation can be made if o vs. frequency is 
measured up to fairly low frequencies, and oa is 
then extrapolated towards zero frequency. In 
homogeneous liquids such precaution is not 
necessary, as specific losses at power and acous- 
tical frequencies are zero; thus the correct d.c. 
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conductivity can be measured by means of 
alternating current. 

Definite deviations from this simple behavior 
occur if the frequency is approaching radio- 
frequencies and higher. In such cases rotation of 
dipolar constituents causes specific losses at 
certain characteristic frequencies. This subject 
will not be pursued here ; more will be said about 
it in connection with solids in Chapter XIV of 
this monograph. It will suffice to point out that 
the characteristic frequency fo is connected with 
the viscosity » of the liquid by the relation 
deduced by Debye 


fo=kT/Arnr’, (107) 


where r=radius of dipolar molecule. Measure- 
ment of fo and knowledge of the viscosity allow 
the radius of the rotating molecules to be 
computed. 

The measured a.c. conductivity differs from 
the d.c. value also in cases of inhomogeneity in 
both a sol and a gel. An interpretation of results 
in these cases might be very involved, but 
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Fic. 71. Activity coefficient vs. concentration of hydro- 
chloric acid in water-dioxane mixtures, after Harned and 
Morrison. (1) Pure water; (2) 20 percent dioxane; (3) 45 
percent dioxane; (4) 70 percent dioxane. 
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enlightening if successful. Inhomogeneous sys- 
tems of this kind were successfully investigated 
in a series of papers by J. 
workers!’ on 


D. Piper and co- 
the conductivity of insulating 
liquids to which electrolytes of known chemical 
constitution and in 
added. 


controlled amounts were 

Besides d.c. and a.c. measurements of con- 
ductivity a third method is sometimes used with 
success, i.e., transient measurement. Carried out 
some time ago by J. B. Whitehead, and recently 
by the Balsbaugh group, it proved very useful. 
The writer's tried to compute absolute mobilities 
from such current-time curves. An equation was 
set up, describing the characteristic decrease of 
conductivity with time on the 
accumulation of ions near the electrodes. Figure 


basis of an 


70 shows two experimental curves obtained by 
Whitehead on two particular oils. The corre- 
sponding full lines are drawn according to the 
equation mentioned, which then allows the mo- 
bilities of the cations and anions to be computed. 
In case of the lower curve the two mobilities 
obtained were : 3.1 10-7 and 0.37 X 1077 cm/sec. 
per volt/cm. Knowing the viscosity (1.1 poise) 
and assuming z= 1, the two radii computed from 
Eq. (105) are 2.51077 and 2.1 X10-* cm. From 
Eq. (104), again assuming z=1, the number of 
charge units becomes 2X10" per cm*. 

Thus, the role of viscosity in the electrical 
conductivity of liquids has_ been 
demonstrated in this section, showing that a 


insulating 


number of other factors, too, are involved in 
establishing a definite conductivity. Because of 
the relatively large numbers of variables, con- 
ductivity and viscosity measurements of the 
type described are not sufficient for the evalua- 
tion of all the data. It will be imperative to use 
entirely different methods as well. 

The ultracentrifuge, so useful in the field of 
synthetic polymers (Section 30) for evaluating 
particle sizes, does not hold out much hope in 
this connection, because of the small number of 
charge carriers. The use of the electron micro- 
scope may be more feasible, because the visibility 
of the particles might allow by some means an 
identification of the ions among other neutral 
particles. The success of the instrument in other 
fields—viruses,'*® for instance—certainly invites 
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its use in the present problem as well. If the 
migration in an electrical field could be observed 
under the instrument—a technique yet to be 
developed, but not essentially impossible—then 
the mobility m and radius 7, and thus the charge 
z, would become known, at least for particles 
above a certain size. 
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Fic. 72. E.m.f.’s of cells containing boundaries copper- 
copper oleate solutions in insulating liquids, after Gemant. 


Another method that in certain instances 
would give directly the last remaining variable, 
namely the number of particles n, is the electro- 
chemical potentiometric method. In aqueous and 
alcoholic solutions this method gives excellent 
results. It measures specifically the concentra- 
tion, more rigorously: activity of certain ions, 
according to the reversible electrode used. In 
the cell 


platinum (hydrogen) /liquid containing 
hydrochloric acid in solution / 
(silver chloride) silver, 


for instance, the hydrogen ion activity can be 
obtained by measuring the e.m.f. of the cell. 
Figure 71 after Harned and Morrison” should 
serve as an indication that precision data can be 
obtained on dioxane-water solvents of dielectric 
constants as low as 17.7. The data deduced from 
e.m.f. measurements on the cell above connect 
the logarithm of the activity coefficient f (a meas- 
ure of the Debye effect) with the concentration. 
The dotted lines refer to coefficients computed 
from equations valid for dilute solutions, the 
agreement with experiment being very good. 
With such data on hand the pH of any HCl 
solution in the dioxane-solvents, in other words 
the number m for the H+ and Cl’ ions, can be 
obtained by means of an e.m.f. measurement. 

The writer* tried to obtain reproducible e.m.f. 
data on liquids of even lower dielectric constant. 
Figure 72 refers to the cell copper/liquid /liquid 
(copper oleate) /copper in xylene-kerosene-phenol 
mixtures (positive if right end of cell is positive). 
An interpretation of these data particularly of 
the reversed (negative) e.m.f.’s has not yet been 
accomplished. A continuation of this research 
with better defined systems is now under way, 
and the transference number (ratio of ‘ionic 
mobilities) is computed from data on concentra- 
tion cells using solvents of a dielectric constant 
as low as 2.7. 
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Bactericidal Radiation 





By A. PAULUS 


Westinghouse Electric and Manufacturing Company, Lamp Division, Pittsburgh, Pennsylvania 


CIENCE is gradually selecting and studying 
various groups of waves in the electromag- 
netic spectrum. These groups are studied for the 
purpose of finding their characteristics and their 
application to human welfare. While the electro- 
magnetic spectrum extends to and beyond the 
infinitely short cosmic waves on one side and 
through the visible and heat regions and radio 
waves of various lengths on the other side, con- 
sideration is here given to a small group of waves 
in the ultraviolet region (Fig. 1). 

Just as visible light is made up of different 
wave-lengths which are separated by a prism into 
the rainbow colors, so the ultraviolet region may 
be separated into groups which have widely 


differing effects. Beyond the visible, the group_ 


of waves of 4000 to about 3300 angstroms are 
chemically active, affect photographic plates and, 
in common with still shorter wave-lengths, cause 
fluorescent and phosphorescent materials to glow. 

A second group measuring 3300 to 2850 ang- 
stroms are biologically effective causing ery- 
themal and tanning reactions, increasing vitamin 
D, and preventing rickets. The effective bac- 
tericidal wave-lengths 2850 to 2000 angstroms 
have been receiving much attention during the 


past few years as the development of suitable 
generators makes possible the practical applica- 
tion of bactericidal radiations. 

Below 2000 angstroms is a group of wave- 
iengths which produce ozone in the air. A small 
proportion of these ozone producing waves, to- 
gether with the more active bactericidal waves, 
produces the best conditions for the destruction 
of bacteria and mold. 


EFFECT OF ULTRAVIOLET ON BACTERIA 


In 1877 Downes and Blunt showed that bac- 
teria were killed by sunlight. The lethal effect 
was attributed to mild and extended oxidation or 
in short to heat. Later the killing effect was 
attributed to the production of hydrogen per- 
oxide in the aqueous media. Both these early 
theories have been disproved by experiments in 
which the lethal effect was produced independ- 
ently of heating or oxygen and where the chemi- 
cal action in the media was negligible. As a 
result of studies by Clark, Norton, and others 
and confirmed by Bayne, Jones, and Van der 
Lingen, it has been generally accepted that the 
bactericidal effect of selected wave-lengths of 
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bactericidal band of the ultraviolet spectrum. 


Fic. 1. The radiant energy spectrum. The Sterilamp is an efficient generator of radiation in the 
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Fic. 2. Progressive stages in 
the killing of bacteria by short- 
wave radiation. At the upper left 
are normal paramecia. At the 
upper right are two paramecia 
after thirty seconds of irradia- 
tion, showing distension just be- 
ginning. Distension continues in 
the third view, followed by burst- 
ing of the cell wall. Death of a 
bacterium is caused by a change 
in the osmotic tension of the cell 
wall, which causes the organism 
to swell and break, as shown in 
the lower right-hand corner. 














ultraviolet is the result of direct photo-chemical 
action on the organisms. The disintegration of 
micro-organisms under radiations, principally of 
2537 angstroms, is shown in Fig. 2, where para- 
mecia which are single cell organisms found in 
stagnant water, are exposed to the radiations 
from a Sterilamp. The picture shows how the 
radiations cause the organism to swell, blister, 
and finally disintegrate. Bacteria are so small 
that the changes which take place because of 
irradiation cannot be shown by an ordinary 
microscope, but they are probably similar to the 
changes in the much larger paramecia. 


ULTRAVIOLET GENERATORS 


The sun is considered the only natural source 
of ultraviolet but only a small part of the short 
wave-lengths reach the earth because of the 
filtering action of the atmosphere. The shortest 
wave-lengths which can penetrate the atmos- 
pheric blanket measure about 2900 angstroms, 
and the amount of this radiation is only about 
one millionth as great as that at 3130 angstroms. 
The ultraviolet radiation from the sun in the 
bactericidal region is great enough to have been 
noted, but it is weak when compared to the out- 
put of many artificial sources. 

Practically all artificial sources of ultraviolet 
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radiation, with the exception of the carbon or 
iron electrode arc lamps, have one element in 
common—mercury vapor. Electrical discharges 
through mercury vapor produce more radiation 
in the ultraviolet region than similar discharges 
through any other known vapor or gas. The 
earliest and best known source of ultraviolet is 
the quartz mercury arc lamp which has an 
envelope of fused quartz and at least one elec- 
trode of mercury. The output of ultraviolet can 
be shifted by varying the temperature and 
pressure of the mercury vapor. Increasing the 
temperature and pressure produces more of the 
longer wave-lengths and the high pressure imer- 
cury lamp becomes an efficient source of visible 
light. If the pressure is kept low, that is, if it can 
be measured by a few millimeters of mercury 
column, the radiation produced by an electric 
current is very largely in the bactericidal region 
with maximum output at 2537 angstrom units. 
These low temperature, low pressure condi- 
tions are Mherent in the Sterilamp, Fig. 3, and 
make it an efficient bactericidal agent, which is 
inexpensive to operate, has long life with small 
change in output, does not heat up the surround- 
ing air, and gives wide distribution of radiation. 
The glass tube which forms the lamp is made of 
glass which cuts off the radiation slightly below 
2000 angstroms, thus limiting the production of 
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ozone so that the concentration does not rise to 
more than one part in three or four million parts 
of air. It that radiation of 2537 
angstroms breaks up ozone and thus acts as an 


is also true 


automatic regulator of ozone concentrations. The 
slight amount of ozone developed is helpful in 
controlling the development of mold and _ bac- 
teria in spaces which are shaded from direct 
radiation. 





Fic. 3. Sterilamps. 


The electrical discharge in a Sterilamp is 
between unheated electrodes, because this per- 
mits simple construction of the lamp and the 
current controlling transformer, permits control 
of the output of the lamp by changing the 
primary voltage of the transformer, makes a 
lamp which will operate satisfactorily under 
varying temperature conditions, and will have a 
long effective life. Deterioration of output comes 
gradually because of changes in the glass known 
as solarization and to a slight extent to blacken- 
ing of the tube by vaporization of the electrode 
material. 


MEASUREMENT OF ULTRAVIOLET 


Much of the recent advance in the application 
of ultraviolet radiation is due to improvements in 
the equipment and methods for measuring the 
different groups of wave-lengths in the ultra- 
violet spectrum. Three photo-cells developed by 
Dr. H. C. Rentschler, head of Westinghouse 
Research Department, cover the principal groups 

of wave-lengths in the ultraviolet spectrum. For 
measurement of the almost complete ultraviolet 
region there is a photo-cell having a thorium 
cathode which is responsive to a spectral range 

from 3750 angstroms to about 200 angstroms, 
with a maximum response at 2600 angstroms. 

Another photo-cell, made up with a tantalum 
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cathode, responds to wave-lengths from 2950 to 
2000 angstroms with maximum sensitivity at 
about 2400 angstroms. This photo-cell is useful 
in measuring the bactericidal region of ultra- 
violet as it is not sensitive to the wave-lengths 
near the visible edge of the ultraviolet spectrum. 
A third photo-cell, with a platinum coated 
cathode, is not affected by radiations above 2000 
angstroms and, therefore, is capable of measuring 
the wave-lengths which produce ozone in the air. 

The ultraviolet radiations in the regions affect- 
ing these photo-cells cause a minute electric 
current to flow through the tube and charge 
a condenser which discharges through a relay 
tube which in turn causes a larger condenser to 
discharge through the coil of a relay which 
operates a counting mechanism with an audible 
click. After discharge, 
charge up and the action is repeated as long as 
ultraviolet light reaches the electrodes and cur- 
rent flows through the photo-cell. As each click 
represents a definite amount of ultraviolet radia- 
tion for which there is no generally accepted unit, 
the meter is called a click meter (Fig. 4). Each 
click may represent any desired quantity of 
radiation by varying the size of the condenser 
in the circuit of the photo-tube. However, the 
standard meter is calibrated so that one click 
represents 220 micro-watt seconds. 


the condensers again 














Fic. 4. The new Westinghouse portable a.c. ultraviolet 
meter. Shown at the left is the photo-tube in its separate 
housing, which fits into the front compartment of the 
carrying case when not in use. 
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Another small indicating meter is an adapta- 
tion of the foot candle meter. A filter which cuts 
out about 95 percent of the visible light but 
permits passage of ultraviolet is placed over the 
window of the meter, and a fluorescent screen is 
placed next to the filter. Ultraviolet light which 
strikes the screen is converted into visible light 
and acts on the light sensitive cell which operates 
the pointer. This meter is very convenient for 
testing ultraviolet sources, such as the Sterilamp, 
which give off little heat but it is not reliable for 
hot sources, as the cell is sensitive to heat as 
well as light, and the meter must be placed very 
close to or in contact with the lamp to get a 
reading. 


APPLICATIONS OF STERILAMPS 


The development of an ultraviolet generator 
which operates at temperatures only a few 
degrees above the surrounding air requires but a 
small amount of power, has a fairly uniform out- 
put and long life, opens up almost unlimited 
possibilities for usefulness. Early applications of 
Sterilamps were for the preservation of meat held 
in refrigerators. The radiations are effective in 
preventing or greatly reducing the development 








Fic. 5. Beef may be made uniformly tender and palatable 
in three days in this room. Weekly capacity of meat 
tenderizing in this room by the Westinghouse Tenderay 
Process is approximately 144,000 pounds. The ultraviolet 
lamps on the ceiling bathe the meat and air in the room 
with invisible rays deadly to mold and bacteria, hence 
prevent spoilage. A relative humidity kept at 90 percent 
prevents dehydration, which causes shrinkage and loss of 
juiciness in meat. A temperature of 60°F speeds up the 
natural chemical reactions which break down tough tissues 
and bring about tenderizing. Those elements, together with 
a scientifically predetermined technique in handling, make 
up the Tenderay Process. 
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Fic. 6. Sterilamps mounted above the operating table 
in this operating room in the Duke Hospital, Durham, 
North Carolina, reduce bacterial count in the air and pro- 
vide a virtually sterile environment on and around the 
operating site. 


of mold and slime during the storage and market- 
ing periods. Installations of Sterilamps in retail 
meat stores reduce trimming losses and, in case 
of failure of the refrigeration equipment, insure 
against meat spoilage until repairs can be made. 
Displays in refrigerated cases are protected from 
contamination by air-borne organisms, and the 
sales appeal is strengthened. The meat packer 
and wholesaler use Sterilamps to keep their 
products free of organisms which cause spoilage 
and thus supply the retail trade with meats 
having low bacterial contamination and, conse- 
quently, of superior keeping qualities and flavor 
(Fig. 5). Protection of meat in this way has made 
economically possible the tenderizing of beef in 
a few days at temperatures of 60°F or more. 
Thus, only about one-tenth of the time is re- 
quired by this new and patented process com- 
pared to the former methods of aging at about 
34°F. 

In hospital operating rooms, Sterilamps are 
mounted over the operating table (Fig. 6) to 
reduce the danger of infection of surgical wounds 
by organisms floating in the air or exhaled by 
the operating staff. The records in more than 
4000 operations in one hospital show that in 
1782 operations performed without radiation there 
resulted 207 infections and six deaths, while 
2463 operations of the same types resulted in 
only six infections and no deaths. Time required 
for healing was also reduced by more than half. 
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In wards and nurseries the use of ultraviolet 
radiation has reduced spread of infectious diseases 
from person to person. Evidence has also been 
accumulated to show that selected radiation in 
schools and offices has reduced or eliminated the 
spread of colds and other diseases which have a 
tendency to become epidemic. 


Sterilamps have been applied in dairy barns, 
poultry houses, breweries, bottlers, bakeries, 
pharmaceutical houses, air conditioning systems, 
and many other places. A new industry has thus 
been developed that has as its primary purpose 
the protection of products and people from the 
unseen dangers ever present in the air. 





Résumés of Recent Research 








Statistical Theory 
of Rubber 


During the past sev- 
veral years considerable 
attention has been paid 
to the behavior of rubber-like materials. Re- 
cently F. T. Wall! developed a statistical theory 
which accounts for some of the properties of 
rubber. This theory postulates that rubber is 
made up of very long chain molecules which, 
because of rotation about the carbon-carbon val- 
ence bonds, can assume various shapes. 

In Fig. 1 are illustrated three possible shapes 
for some relatively short molecules. In that dia- 
gram the x values are the components of the 
distances between the ends of the molecules 
along the direction in which the rubber is to be 
stretched. It is then recognized that the distribu- 
tion of these x values will be subject to the normal 

















p(x) / 





Fic. 2. S: Distribution function for unstretched rubber. 
S’: Distribution function for stretched rubber with 100 
percent elongation. 


1F, T. Wall, J. Chem. Phys. 10, 132 (1942). 
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probability distribution function, represented by 
curve S in Fig. 2. It will be observed that most 
of the molecules have small x values and that 
only a few are in a straightened out condition. 
When a macroscopic piece of rubber is 
stretched, the distribution S will also be stretched 
into a new distribution S’ (Fig. 2). The statistical 
problem which Wall then solved involved the 
calculation of the probability of finding the sys- 
tem with a distribution S’ when its most probable 
distribution is S. This result he related to the 
entropy of stretching in accordance with the 
statistical theory of entropy. Then by well- 
known thermodynamic equations he deduced an 
equation of state for a cylindrical piece of rubber, 


namely: 
NkT sl lo 
aC) 
lo Nlo 1 


where + is the tension of the rubber, N the num- 
ber of molecules, k the Boltzmann constant, 7 
the absolute temperature, /) the length of the 
unstretched rubber, and / the length under a 
tension +. It was observed that this equation 
agrees better with observation than some other 
equations and that it applies to compression as 
well as to extension. 


Industrial Noises 
and Industrial 
Deafness 


The importance of oc- 
cupational deafness be- 
comes obvious when the 
hearing loss among boil- 
ermakers is compared (Fig. 1) to the so-called 
“normal” hearing loss of the average population 
(as measured by the Bell Telephone Laboratories 
at the 1939 World’s Fairs). 

In recent years a survey was conducted in 
France with the purpose of determining if there 
was any relationship between the intensity and 
the frequency distribution of industrial noises 
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F Fic. 1. Comparison between ‘normal’ and industrial 
hearing loss. Permanent hearing loss among boiler- 
makers after more than 15 years of employment. —-—— 
Average hearing loss in 40-49 year age group of the 
general population. Average hearing loss in 50-59 
year age group of the general population. 





and the hearing loss commonly encountered 
among industrial workers.* 

It was possible to show some degree of correla- 
tion between the average noise level and the 
temporary and permanent hearing losses. This is 
illustrated by Fig. 2 which shows furthermore 
how the permanent hearing loss is almost 
completely localized in the region above 1500 
cycles (with the maximum loss occurring at 6000 
cycles), while the most important components 
of industrial noises are usually to be found below 
this frequency limit. The very small permanent 
loss for low frequencies points to a very rapid 
recuperation from temporary losses in this 
region. This is in marked contrast to what 
happens in the higher frequency range where 
some correlation is found to exist between the 
temporary and permanent losses. 

While the mechanism responsible for the high 
frequency deafness as a consequence of exposure 
to low frequency noise is still unknown, the 
phenomenon in itself is significant. Inasmuch as 
the particular high frequency range for which 
the permanent hearing loss will be maximum is 
also the one which is first affected by noise 
deafness, it becomes possible to detect at an 


_ *Some of the results of this investigation were published 
in an article by W. A. Rosenblith, J. Acous. Soc. Am. 13, 
220 (1942). 
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Fic. 2. Hearing loss among various categories of per- 





sonnel after 15 years of employment. Boilermakers 
(average noise level: 90 db). .... Black- and ironsmiths 
(average noise level: 80 db). ---— Machinists (average 


noise level: 75 db). 


early stage the danger of impending deafness 
among industrial workers, aviators, etc. It 
would seem desirable that this localized hearing 
loss caused by industrial noises should be taken 
into account in the design of ear defenders. 

City Noise* The many mechanical 
devices that were in- 
vented during the '20’s, which increased effi- 
ciency and output, so intrigued the general 
public at first that they did not realize the 
damage that was being done to their health by 
the noise inherent in the operation of these 
many mechanical devices. The noise of the auto- 
matic riveter, the compressed air drill, the sub- 
way turnstile, to mention but a few of the 
newer devices, was added to the clamor of the 
elevated, subway, and surface cars, screeching 
brakes, terrifying automobile horns, chain driven 
trucks, and all the other noises encountered in a 
large industrial community. 

As these noises continued to multiply the 
public became conscious of the effects of noise 
upon their nervous systems, with the result that 
complaints from citizens began to pile up in the 
Health Department. However, because there was 
no scientific evidence that noise was harmful, the 
Department of Health had little success when 


* Shirley W. Wynne, J. Acous. Soc. Am. 13, 214 (1942). 
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persons making unnecessary noises were sum- 
moned to Court. 

In 1929 the then Commissioner of Health, 
Dr. Shirley W. Wynne, appointed a committee 
made up of outstanding men in the various pro- 
fessions—physicians, lawyers, automotive engi- 
neers, and acoustical engineers—to make a study 
of noise. The physicians were asked to study the 
effects of noise upon the human system. Experi- 
ments were carried out which proved that inter- 
cranial pressure was increased fourfold, blood 
was raised, heart action was made 
irregular, and hearing was damaged. The engi- 
neers measured the amount of noise made by 
automobiles, elevated and subway trains, and 
suggested practical measures for lessening these 


pressure 





noises. Many of these recommendations have 
been carried out. The lawyers devised a city 
ordinance which would, if enacted into the law, 
make possible the infliction of small fines, and 
yet provide for the expeditious handling of 
violators of the ordinance. Because of lack of 
public sentiment the ordinance was not enacted. 

The results of this study proved that noise is 
detrimental to health, that it interferes with the 
efficiency of workers, and that most noise can be 
prevented. 

The general public lacks an appreciation of the 
effect of noise upon the human organism, and 
until they have been educated and public opinion 
aroused, little progress will be made in abating 
unnecessary noises. 





New Books 








Theory of Gaseous Conduction and Electronics 


By Freperick A. MAXFIELD AND R. RALPH BENE- 
pict. Pp. 483+xiv, Figs. 241, 15x23 cm. McGraw- 
Hill Book Company, Inc., New York, 1941. Price 
$4.50. 

A new book, especially in as rapidly developing a field 
as that of conduction in gases, must be catalogued either 
as suited to the needs of the undergraduate being intro- 
duced to the subject, or as essentially a reference for those 
who have already mastered its fundamentals. This book, 
intended as an undergraduate text, actually succeeds in 
being something of a compromise: there are both explana- 
tions suited to the undergraduate, and tabulations and 
calculation methods which form a convenient reference 
for the expert. 

The thoroughness with which the most diverse aspects 
of the subject are gone into (nine pages of discussion are 
given, for example, of the Fermi-Dirac statistics of elec- 
trons in metals) is praiseworthy, but might discourage 
students not well grounded in mathematics and electricity. 
This difficulty would be heightened by the use of m.k.s. 
units in schools using pre-requisite physics texts based on 
the c.g.s. system. Qualitative explanations are, however, 
also abundant and well-written. The value of the book 
as a text is much enhanced by a large choice of problems. 

The authors have placed most of their emphasis on 
theoretical principles, although recently published experi- 
mental material has not been neglected. In regard to 
applications, even such relatively recent developments as 
double-grid thyratrons, oil-poor circuit breakers, and 
fluorescent lamps are discussed. 

One disadvantage of this book as a reference is that 
sources of material presented are at times difficult to find. 
In setting forth the fundamental principles involved in 
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gaseous discharges and showing how these apply to a 
great variety of commercial devices, however, the authors 
have achieved a high degree of success. The book is 
recommended for both engineers and physicists interested 
in conduction in gases. 
Joun E. WHITE 
University of Pittsburgh 


Temperature Measurement and Control 


By Ropert L. WEBER. Pp. 430+x, Figs. 183, 15} 22 
cm. The Blakiston Company, Philadelphia, 1941. 
Price $4.00. 

This book presents the material on temperature meas- 
urement offered to students with one year of physics at 
The Pennsylvania State College. The volume contains, 
in addition to the descriptions of methods of temperature 
measurements, a sufficiently detailed discussion of the 
basic theory to give the student an adequate background. 
For those who prefer to approach the subject of heat from 
the standpoint of temperature measurements the book 
will serve as a good basic text. References at the end of 
each chapter are a guide to more detailed discussion of 
the material. A set of problems will be found at the end 
of many chapters. 

The book is divided into two parts. The first deals with 
the general theory and the second constitutes a laboratory 
guide. While the reviewer prefers to have the experiments 
dispersed throughout the body of the text so that the 
demarcation between the theory and the laboratory work 
is less obvious, he recognizes that this is just an opinion 
and is likely to cause confusion unless sufficient equipment 
is available so that the two can be coordinated closely. 

The book is a valuable addition to the textbook material 
on the subject. The diagrams are excellent. The type and 
binding are good, but the paper shows too much glare 
for comfort. 

H. B. WAHLIN 


University of Wisconsin 
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Recruiting Continues for Technologists 

About a year ago the Federal Civil Service Commission 
began a recruiting drive for technologists. The number of 
persons on the employment lists now established is not 
sufficient for today’s wartime needs. Accordingly, the 
Commission has recently reissued the announcement of 
these opportunities for Government employment with 
slightly modified requirements. Salaries for the positions 
range from $2000 to $5600 a year. Applications may be 
filed until further notice with the Commission’s Washing- 
ton office. 

By technologist the Commission means a person ex- 
perienced in ‘‘the necessary production, engineering, and 
scientific research work essential for the successful opera- 
tion of an industrial plant, where such plant operation is 
based upon a thorough and expert knowledge of a branch 
of an applied science,’ such as explosives, fuels, plastics, 
rubber, minerals, or textiles. 

For the Junior positions ($2000 a year), applicants will 
no longer have to take a written test, and, as before, no 
written test will be given for the higher positions. The 
maximum age limit for all the grades has been raised to 
60 years. 

Copies of the announcement (No. 188), and the forms for 
applying, may be obtained at first- and second-class post 
offices, or from the Civil Service Commission, Wash- 
ington, D. C. 

* 
Bulletin of Mathematical Biophysics 


The June, 1942 issue of the Bulletin of Mathematical 
Biophysics, edited by N. Rashevsky and published by the 
University of Chicago Press, will contain the following 
articles: 

A Mathematical Analysis of Elongation and Constriction in Cell 
Division, H. D. LANDAHL. 

A Fundamental Form for the Differential Equation of Colonial and 
Organism Growth, NATHAN W. SHOCK AND MANUEL F. MorALEs. 

A Note on the Diffusion of Electrolytes in Cells, ALStoN S. HouseE- 
HOLDER AND ROBERT R. WILLIAMSON. 


Periodic Phenomena in the Interaction of Two Neurons, GEorRGE 
SACHER, 


A Theory of Electrical Polarity in Cells, Ropert R. WILLIAMSON 
AND INGRAM BLOCH. 


An Alternate Approach to the Mathematical Biophysics of Perception 
of Combinations of Musical Tones, N. RASHEVSKY. 


* 
Charles A. Coffin Fellowships 

The General Electric Company has announced the 
winners of the fellowships granted annually by the Charles 
A. Coffin Foundation to enable outstanding young gradu- 
ates to continue or undertake research work in educational 
institutions either in this country or abroad. The fields in 
which awards are made include electricity, physics, and 
physical chemistry. This year’s successful applicants are: 

Boris D. Abramis, University of California, to continue 
at same school the theoretical and experimental study of 
hyper-frequency wave guides; Robert W. Hull, Yale 
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University, to study at Massachusetts Institute of Tech- 
nology the chemical nature of thermionic emitters, such 
as oxide cathodes, and the effects of impurities on their 
properties; Edward R. Kane, Union College, to conduct at 
Massachusetts Institute of Technology an experimental 
study of the absolute temperature scale; Hershel Marko- 
vitz, University of Pittsburgh, to attend Columbia Uni- 
versity to work on the theory of the use of the distribution 
functions in calculating the thermodynamic properties of 
liquids; and Arthur E. Snowdon, Carnegie Institute of 
Technology, to study at the same school the design of 
variable, smooth model transmission line. 


* 


Recent Appointment 

Dr. Louis M. Heil, research associate in science at the 
University of Chicago, has been appointed Head of the 
Department of Physics at Cooper Union with the rank 
of full professor. He succeeds Professor Albert Ball, who 
will retire on June 30. 

Emmett M. Irwin, consulting engineer and chief elec- 
trical engineer in charge of design and construction of all 
controls and drives for the 200-inch telescope being con- 
structed at Palomar Mountain, California, has been placed 
in charge of the newly created Pacific Coast Branch of the 
Waugh Laboratories of New York. 

* 
Necrology 


Sir William Bragg, Director of the Davy-Faraday Re- 
search Laboratory, Fullerian Professor of Chemistry and 
Director of the Royal Institution, London, and President 
of the Royal Society from 1935 to 1940, died on March 13 
in his eightieth year. 


* 
Calendar of Meetings 
May 
11-13 American Institute of Chemical Engineers, Boston, Massa- 


chusetts 
15-16 Acoustical Society of America, Ann Arbor, Michigan 
31-June 5 Society of Automotive Engineers, White Sulphur Springs, 
West Virginia 
June 
8-11 American Society of Mechanical Engineers, Cleveland, Ohio 
8-12 American Medical Association, Atlantic City, New Jersey 
15-20 American Association for the Advancement of Science, Salt 
Lake City, Utah 
16 American Association of Physics Teachers, Salt Lake City, 
Utah 
18-20 National Colloid Symposium, Boulder, Colorado 
22-26 American Institute of Electrical Engineers, Chicago, Illinois 
25-27. American Physical Society, State College, Pennsylvania 
25-27 American Association of Physics Teachers, State College, 
Pennsylvania 
29-July 1 Institute of Radio Engineers, Cleveland, Ohio 
29-July 1 American Society of Heating and Ventilating Engineers, 
St. Paul, Minnesota 
July 
11 American Physical Society, Pacific Coast, Berkeley, California 
22-24 American Society of Civil Engineers, Spokane, Washington 
September 
7-9 Mathematical Association of America, Ithaca, New York 
7-11 American Chemical Society, Buffalo, New York 
8-11 American Mathematical Society, Ithaca, New York 
9-11 American Institute of Electrical Engineers, Vancouver, B. C. 
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On the State of Stress in Thick Bars{ 


M. M. Frocut! anp M. M. LEVEN? 
Carnegie Institute of Technology, Pittsburgh, Pennsylvania 


(Received February 2, 1942) 


The view is now generally held that in bars containing circular holes or other discontinuities 
the state of stress ceases to be two-dimensional when the ratio of the radius of the hole r to 
the thickness of the bar ¢ is less than unity. It is further believed that the factor of stress 
concentration in thick bars is materially smaller than in thin bars. This view also forms the 
basis for the doubts which have arisen regarding the reliability of certain theoretical work on 
stress concentrations produced by very small holes. In this paper results are reported from 
a photoelastic investigation using both elastic and frozen (see reference 4, Chapter 10) stress 
patterns and showing that the present point of view is erroneous. Experimental results show 
that no deviation from a state of plane stress can be found in bars with central circular holes 
when r/t is of the order of magnitude of about } instead of unity. Experiments further show 
that the factors of stress concentration remain essentially the same in thick as in thin bars. 
The present work may also be taken as experimental corroboration of the theoretical value for 
the factor of stress concentration obtained by Neuber for a very small hole in a bar in tension 


or compression. 


N the theory of elasticity two types of two- 
dimensional problems are recognized: one of 
plane stress and one of plane strain. A state of 
plane stress is approached in a very thin bar 
without discontinuities in tension or compression 
and a state of plane strain is approached in the 


Stress Concentration 





0 04 0.8 1.2 1.6 20 
Size of Hole 


. iam Hole 
Ratios Pam ole “Plate 
Fic. 1. Effect of thickness upon factor of stress concen- 
tration in which 2r/D is kept constant (Vose). 





t Presented at the Fourteenth Semi-Annual Eastern 
Photoelasticity Conference, held at Yale University, 
New Haven, Connecticut, December 6, 1941, under the 
auspices of the Civil and Mechanical Engineering De- 
partments. 

1 Associate Professor of Mechanics in charge of Photo- 
elastic Laboratory, Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania. 

2 Mechanician and Research Assistant, Photoelastic 
Laboratory, Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania. 
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central transverse section of an infinitely long 
cylinder subjected to uniform loading along a 
line parallel to its axis. 

When a thin bar in tension has a discontinuity, 
such as a hole, groove, or fillet the state of stress 
may be affected by the discontinuity. It is both 
of theoretical and practical interest to examine 
the effect of the thickness on the state of stress 
resulting from such discontinuities. 

A survey of the literature shows that at least 
one experimental attempt to answer this question 
was made in recent years. In 1935 Vose* published 
a curve purporting to show the effect of the 
thickness of a bar with a circular hole upon the 
factors of stress concentrations, and he concluded 
that there is a pronounced reduction in the factor 
of stress concentration as the thickness increases. 

Vose’s results are shown in Fig. 1. It is seen, 
for example, that for a hole in a bar in tension 
the factor of stress concentration drops from 
nearly 3.00 to 1.90 as 2r/t varies from 2.00 to 
approximately 0.1, the reduction being 37 
percent. 

As far back as 1935 the authors had reason to 

3 See R. W. Vose, ‘‘An application of the interferometer 


Cc al in photoelasticity,” J. App. Mech. 57, A-101 
(1935). 
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Dimensions: 


2r =0.189” 2r/D =0.140 P =298 lbs. 

D =1.355” 2r/t =0.223 Min. area =0.945 sq. in 
d =1.166” SAy =316 Dsi. 
‘=0.811” 2F =106 psi. 


"max 


=8.00; k =8.00/2.98 =2.68; 


NAy =2.98 fringes 
k (from curve) =2.55. 


J. N. Goodier* who in reviewing Neuber’s book 
Kerbspannungslehre states: ‘“‘When the radius 
of the hole, or in notches, the radius of the 
curvature, is small, it is likely in photoelastic 
models to be of the same order as the thickness 
of the plate. The state of stress is then no longer 
even approximately two-dimensional, as Vose 
has recently pointed out.’’ Goodier makes this 
the basis for a likelihood ‘‘to resist the author’s 
[Neuber’s] conclusion that his final formulas 
are satisfactorily confirmed.”’ 

According to this view, it would follow that 
in a }” plate three-dimensional effects exist 
around a central circular hole $ in diameter, 
or to take a more extreme case, in order to 
eliminate three-dimensional stresses around a 
: hole in tension the bar would have to be 
thinner than ;g of an inch. 

The purpose of the present paper is to present 
the results from a recent photoelastic investiga- 
tion of this problem using thicker bars, and the 
newer method of frozen stress patterns, as well 
as elastic stress patterns at room temperature. 


Error = +5 percent. 


Fic. 2. Thick bar with central circular hole subjected to 
pure tension. 











TABLE I. 
Factor 
of 
stress 
Factor con- 
of cen- 
Diam- Thick- stress tra- 
eter Width ness con-_ tion 
of of of cen- for 
I'ype hole bar bar tra- thin 
ot 2r D t tion plates 
stresses in. in. in. 2r/D 2r/t k ko k/ko 
Elastic 
room 
temp. 0.189 0.632 0.640 0.300 0.295 2.32 2.29 1.01 
Elastic 
room 


temp. 0.158 1.540 0.811 0.102 0.195 2.44 2.62 0.93 
Elastic 
room 


temp. 0.189 1.355 0.811 0.140 0.223 2.68 2.55 1.05 
Frozen 


260°F 0.212 1.342 0.805 0.158 0.263 2.43 2.52 0.965 
Frozen 


260°F 0.212) 1.342 0.805 0.158 0.263 2.535 2.52 1.00 


+ From central slice 0273” thick, 
doubt Vose’s conclusion. In our studies of stress 
concentrations carried on at that time, plates 
from 35" to 3” thick were used, and no effect 
of the thickness of the plate on the factors of 
stress concentration was found even for grooves 
ie’ in radius. 

Vose’s conclusions have been accepted by 
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FiG. 2(a). Enlargement of central portion of Fig. 2. 
4J. App. Mech. pp. A-188, 189 (December, 1937). 
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(a) 





Fic. 3. (a). Space effect distributed. (b). Space effect eliminated from left and lower parts. 


THICK RECTANGULAR BAR. ELASTIC PATTERNS 
AT ROOM TEMPERATURES 

Table I shows the dimensions of the bars 
which were pulled in tension by means of axial 
loads. A typical stress pattern is shown in Fig. 2. 
Attention is directed to the clear boundary 
stresses. In view of the relatively great thickness 
of 0.8 inch the high degree of boundary visibility 
around the ;"" diameter hole is indeed remark- 
able. Unless proper precautions are taken it is 
very easy to lose a fringe or two due to shadows 
or space effects. This can happen even in thin 
models with large circular holes, let alone in a 
bar 0.8” thick with a small hole. Inspection of 
Fig. 3 shows one whole fringe was lost because 
of faulty technique in the polariscope. Thus in 
Fig. 3(a) the maximum fringe order is 9 whereas 
in Fig. 3(b) it is 10.° 

Table I shows the resulting factors of stress 
concentration based on elastic patterns at room 
temperature. The values of these factors are seen 
to fluctuate about the corresponding value of the 
factor for a thin plate. The average value is, 
however, in very good agreement with the value 
for a thin plate. These fluctuations indicate 





5’ The proper technique for revealing boundary stresses 
is discussed in detail, Max M. Frocht, Photoelasticity 
(John Wiley & Sons, Inc., New York), Vol. 1, Chapter 11, 
pp. 370, 376. 


310 

















P =21.2 Ibs. 
Min. area =0.931 sq. in. 


Final dimensions: 


r’ =0.106” 


Original dimensions: 


2r =0.189” 


D =1.355” D’ =1.342” Say =22.8 psi. 
d =1.166” d’ =1.157” 2F =3.97 psi. 
¢=0.811” t’ =0.805” nay =5.74 fringes 
2r’/D’ =0.158 k (from 2 dimens. curve) 
2r’ /t’ =0.263 =2.52 


Fic. 4. Thick bar with small circular hole in pure tension 
(frozen stresses viewed through Halowax oil). 
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4 
Final dimensions: 
r’ =0.106” =21.2 lbs. 
D’ =1.342” Sy =22.8 psi. 
d’ =1.157” 2F =S.23 psi. 
t’ =0.612” NAy =4.355 fringes 
Mmax =9-87; k =9.87/4.355 =2.27. 


_Fic. 5. Thick bar with small circular hole in pure tension 
(frozen stresses viewed through Halowax oil). Plate 
thinned down from 0.805” to 0.612”. 


possibly that the pin through which the load 
was applied deflected considerably as the loads 
were increased and bore on, or came in contact 
with, different parts of the bar. 

The above results are obviously at variance 
with those obtained by Vose. The difference is 
in all probability due to the particular method 
employed for determining the boundary stresses. 
It has been pointed out by several observers 
that lateral deformations cannot be used as a 
measure of the sum of the principal stresses in 
thick bars, and Vose seems to have relied 
mainly upon such deformations.’ Additional 
errors are introduced by the fact that the needles 
of the interferometer can never be placed right 
at the boundary. 

The elastic stress patterns by themselves do 
not provide, however, sufficient evidence to 

* This point was stressed by Dr. D. C. Drucker from 


Cornell University in his oral discussion of the present 
paper. 
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warrant a conclusion regarding the distribution 
of the stresses in the bar as a function of thick- 
ness. Strictly interpreted they furnish only data 
for the integrated effect, but they do not spe- 
cifically tell us how the stresses vary from layer 
to layer. In order to answer this question a 
photoelastic stress pattern was frozen into a 
bar 0.8” thick, and this bar was subsequently 
milled down in several stages to give thinner 
bars of 0.612’, 0.392’, and 0.275’, respectively. 
The milling was done equally from both sides 
so that each thinner bar represents a central 
section of the original bar. The corresponding 
stress patterns are shown in Figs. 4-7. Curves 
of (p—gq) in fringes per inch of thickness were 
plotted for several sections. Figure 8 shows one 
set of these curves from which the values of the 
fringe order nm at a particular point could be 
obtained as a function of thickness ¢. The results 
are shown in Fig. 9. Inspection of these curves 
shows that the inner layers carry a somewhat 
greater stress per inch of thickness than the 
outer layers of the bar. The curves are, however, 
straight over the middle half approximately, 
and the maximum deviation from the straight 
line does not exceed 12 percent. 

















Final dimensions: 


r’ =0.106” SAy =22.8 psi. 
D’ =1.342” 2F =8.16 psi. 
d’ =1.157” NAy =2.79 fringes 


’ =0.392” 
Mmax =6-75; k =6.75/2.79 =2.43. 


Fic. 6. Thick bar with small circular hole in pure tension 


(frozen stresses viewed through Halowax oil). Plate 
thinned down from 0.805” to 0.392”. 


311 
































Final dimensions: 


r’ =0.106” P =21.20 Ibs. 


D’ =1.342” Say =22.8 psi. 

d’ =1.157” 2F =11.62 psi. 

t’ =0.275” NAy =1.96 fringes 
Umax =5.00; k =5.00/1.96 =2.54. 


Fic. 7. Thick bar with small circular hole in pure tension 
(frozen stresses viewed through Halowax oil). Plate 
thinned down from 0.805” to 0.275”. 


FACTORS OF STRESS CONCENTRATION 


It should be observed that the frozen stress 
patterns introduce greater difficulties into the 
problem of stress concentration as compared 
with the elastic patterns at room temperature. 
First, the deformation becomes rather large and 
the shape of the transverse hole is materially 
changed; second, the process of annealing or 
freezing produces edge effects which are most 
undesirable as they destroy the boundary 
visibility, which is so essential to the accurate 
determination of stress concentrations. Never- 
theless, the results from the frozen stress pattern 
are interesting and corroborating. Since the 
curves of Fig. 9 show that the greatest stresses 
exist in the central part of the bar we determine 
the maximum stresses around the hole from the 
stress patterns for the inside sections. This 
gives factors of stress concentration of 2.43 and 
2.53 as compared with 2.52 for a thin plate of 
the same face dimensions,’ as shown in Table I. 

? Two-dimensional factors of stress concentration can be 


found in the following: N. C. Riggs and M. M. Frocht, 
Strength of Materia!s, chapter 13 (The Ronald Press 
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In view of the difficulties pointed out above 
these results may be taken to represent good 
support for a conclusion that the factor of stress 
concentration in the thick bar is essentially the 
same as in thin bars having the same ratio 2r/D. 


THICK BAR WITH LARGE RIVET OR PIN 


A bar 0.8” thick and 13” wide was pulled by 
means of a pin 1.00” in diameter and a photo- 
elastic stress pattern was frozen into it. The 
bar was then milled down equally from both 
sides to give several thinner bars. A typical 
frozen stress pattern is shown in Fig. 10. Curves 
of (p—q) were plotted for the different bars 
from which the resulting curve of the fringe 
order n vs. thickness ¢ was obtained. This turned 
out to be a straight line, indicating that the 
small disturbances observed when the hole was 
small are here altogether absent. 


CONCLUSIONS 


1. Experiments show that in thick bars with 
circular holes the inner layers carry a somewhat 
greater stress than the outer layers. This is 
shown by the frozen stress patterns and the 
curves of the fringe order n vs. thickness ¢. 


Y 
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Fic. 8. Curves of (p—gq) in fringes per inch of thickness 
across AC, for slices of different thicknesses. 
Company, New York); M. M. Frocht, ‘Photoelastic 
studies in stress concentrations,” Mech. Eng. J. A. S. M. E. 

(August, 1936). 
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2. The variations in the stresses are, however, 
small and do not exceed 12 percent. 

3. Vose’s conclusions regarding the effect of 
the thickness upon the factors of stress concen- 
tration are definitely not substantiated by our 
work. The curve published by him and shown 
in Fig. 1 is believed to be in error. 

4. Factors of stress concentration for small 
holes in thick bars do not differ by more than 
5 percent from the corresponding factors in a 
thin plate. 
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Fic. 9. Curves of fringe order vs. thickness, based on 
curves such as shown in Fig. 8. To illustrate, the curve 
for point G is obtained by taking the fringe order m; from 
Fig. 8 and multiplying it by the thickness. Thus at point 
G where x/AC=0.7 and t=0.392", n=n,\t=7.2X0.392 
= 2.82. Similarly, at point G for t=0.612”, n=n,t=6.90 
0.612 =4.22. 

5. When the hole becomes large, as in the case 
of the single rivet or pin, experimental data fail 
to show any influence of thickness, indicating 
that the important element is the ratio of the 
dimensions rather than the absolute dimensions 
which determine the state of stress. 

6. For practical design purposes the state of 
stress in thick bars may be assumed to be the 
same as in thin bars and the factors of stress 
concentrations although somewhat reduced ap- 
proach the same values as for thin bars. 
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Original dimensions: Final dimensions: P =21.2 lbs. 
2r =1.004” 2r’ =1.178” Say =77.9 psi. 
D =1.360" D’ =1.330" 2F =5.08 psi. 
2r/D =0.739 d’ =0.340” nay =15.3 fringes 
t’ =0.630” 
2r’/D’ =0.885 


"max =27; k =27/15.3 =1.77; k (from 2-dimens. curve) =1.80. 


Fic. 10. Head of thick bar strained with a single large 
central rivet (frozen stresses viewed through Halowax oil). 
Bar thinned down from 0.800” to 0.630”. 


7. The experiments may be interpreted as 
substantiating Neuber’s theoretical conclusion 
regarding stress concentrations produced by 
very small holes. 
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Measurement of the Initial Inner Permeability of Iron over a Wide 
Radiofrequency Range 


RoBert W. Woops* 
Department of Physics, The University of Chicago, Chicago, Illinois 
(Received January 20, 1942) 


Previous investigations on the variation with frequency of the initial inner permeability uz 
of iron reveals great disparity of results. In this experiment certain refinements of technique 
and of analysis have been introduced which make the accuracy adequate for comparison with 
theory. The method used was to determine the change in resonant frequency of a coil and 
condenser combination caused by insertion of a permeable sample in the field of the coil. This 
change in frequency is interpreted as being caused by change of flux through the coil. The same 
iron wire had been measured previously over a different frequency range by other men working 
in this laboratory. The value of uz was found to be ~60 over the frequency range from 30 kc 
to 40 mc showing no dispersion region between these limits. Results of this method are compared 


with those of other methods. 


I. INTRODUCTION 


HE most commonly used method for 

measuring the initial permeability of ferro- 
magnetic materials at ordinary radiofrequencies 
has been to measure radiofrequency resistance. 
The permeability ue calculated from these data, 
called the outer or resistive permeability, is 
different in value from the inner or inductive 
permeability ur. 

Kreielsheimer! calculated values of yu, from his 
measurements of resistance for the range from 
40 meters to 1000 meters. Several investigators? 
have calculated y, from inductance measure- 
ments for frequencies above 100 megacycles but 
the methods employed have in general not been 
applicable at lower frequencies. Wait* using a 
method due to Belz,‘ obtained some values which 
he called “apparent permeability.’ These values, 
while probably adequate as a check on the 
reported’ anomalous behavior of yw, in the 
100-meter region, was, by reason of inadequate 
analysis, only about one-tenth of the probable 
value of ur. 

In view of the low accuracy of these results 
it was decided to calculate yw, over a more 


* Now at Union College, Lincoln, Nebraska. 

1K. Kreielsheimer, Ann. d. Physik 17, 293 (1933). 

2]. B. Hoag and H. Jones, Phys. Rev. 42, 571 (1932); 
J. L. Glatthart, Phys. Rev. 55, 833 (1939); G. Potapenko 
and R. Sanger, Zeits. f. Physik 104, 779 (1937). 

3G. R. Wait, Phys. Rev. 29, 566 (1927); 32, 967 (1928). 

4M. H. Belz, Phil. Mag. 44, 479 (1922); M. H. Belz, 
Proc. Camb. Phil. Soc. 21, 52 (1922). 

5B. Wedensky and K. Theodortshik, Physik. Zeits. 24, 
216 (1923); W. Kartschagin, Ann. d. Physik 67, 325 (1922). 
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extended range and with more adequate analysis, 


from measurements of inductance. The sample 
of iron was a single strand of the same wire 


used by Glatthart.® 


Il. METHOD 


To measure the permeability uw, at frequency f 
the iron wire sample was placed symmetrically 
along the axis of the coil (Fig. 1) in the frequency 
control tank circuit of a stabilized’ oscillator. 


. 
er 











Fic. 1. Longitudinal section of coil. 


The frequency f was measured. Figure 2 is a 
diagram of the measurement circuit. The 
Franklin oscillator’ was employed because the 
coil in this oscillator has one side grounded, and 
no part of the d.c. grid or plate current flows 
through the coil as in the electron coupled 
oscillator. The other oscillator which provided 
the heterodyne was an electron coupled oscil- 


6 J. L. Glatthart, Phys. Rev. 55, 833 (1939). 

7C. D. Perrine, Radio p. 14 (June, 1939); L. Norton, 
Radio p. 11 (July, 1939); C. E. Berry, Radio p. 50 (Novem- 
ber, 1939); L. Norton, Radio p. 41 (January, 1940); G. F. 
Lampkin, Proc. I.R.E. 149 (March, 1939). 

8 L. Norton, Radio p. 41 (January, 1940). 
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lator. Doubling the frequency of both of these 
oscillators increased their stability and mini- 
mized any tendency toward ‘‘pulling.’’ The beat 
frequency output of the mixer was _ passed 
through a low pass filter, amplified, and applied 
to the horizontal plates of an oscilloscope. The 
Lissajous figure on the screen of the oscilloscope 
enabled the experimenter to set the 80-120 kc 
oscillator to an exact multiple of this beat 
frequency. By zero beating a high harmonic of 
this last-mentioned oscillator against the accu- 
rately calibrated radiofrequency oscillator, it 
was possible to measure the beat frequency to a 
high degree of accuracy. The difference between 
two determinations, once with the iron wire in 
the coil, and once with a copper wire replacement 
of the same dimensions gave the value of df. 
The value of f was measured by direct com- 
parison of the frequency of the test oscillator 
with the calibrated oscillator, again by the zero 
beat method. These are the only measurements 
except the dimensions of the coil, length 2L, 
radius R and of the sample, length 2/=6 cm, 
radius r=0.0088 cm. 

The values of uw, are then calculated from the 
observed dataft by the formula 


ur=[1/Q]{|L2SAdf ]/[f(a—U)]} +P 


where A=7R? and a=7r*; S is a factor which 
allows for radial variation of the magnetic field 
HT of the coil; U is a factor which allows for 
reduction of HZ due to magnetization of the 
sample; Q is a factor which allows for radial 
variation of the instantaneous flux density B 


t In order to keep the value of H constant in all measure- 
ments it seems that a variable resistor should be placed in 
the tank circuit and adjusted so that H be kept constant 
at all frequencies. This was considered and was not done 
for three reasons: (a) The inclusion of a resistor in the 
frequency control tank circuit of an oscillator adversely 
affects the Q of the circuit and decreases the stability of 
operation. (b) Previous measurements of effective permea- 
bility seem to indicate that the permeability does not ap- 
preciably increase above the initial value until the value 
of H exceeds one oersted. (c) Experimentally, at a fixed 
frequency, changes in the capacity of the feedback con- 
denser within the functional limits, which changed the 
voltage across the coil, the current through the coil, and 
consequently the value of H, made no appreciable differ- 
ence in the calculated value of the permeability. The static 
permeability was determined by the ring-ballistic method 
and was found to be 108 at 0.22 oersted. This is comparable 
to the values of effective permeability obtained at radio- 
frequencies, since even at low frequencies, 60 cycles, the 
effective permeability obtained by dynamic methods is only 
a little more than half of the static permeability at initial 
held strengths. The outer permeability was not measured. 
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Fic. 2. Measurement circuit. Values of circuit com- 
ponents have purposely been omitted from the diagram 
because these values of coils, condensers, and chokes had 
to be changed as the frequency was changed. The resistance 
values depend on tube characteristics and may be selected 
from any circuit design tables. 


within the iron sample, due to eddy currents. 
It is a function of f and of uz. P is the analogue 
for copper of Q for iron. 


Ill. DERIVATION OF EQUATION 


From® 
f=[1/2r]V/ {[1/L.C]—[R/2L. ?} 
where L,=effective inductance, C=constant 
tuning capacity, R=constant r-f resistance, 
— 2df/f=(dL./L.|{1—CR?*®/[4L,.—CR?]}} 
=dL,./L, 
since for normal values of C, R, and L,, 
CR?/(4L,.— CR? ]=2.5XK10<1. 
From!® 
L.=L/[(1—4rfLCy] and 4rf?(C+Co)]=1, 
where L=true inductance and C)=distributed 
9A. Hund, Phenomena in High Frequency Systems 
(McGraw-Hill, 1936), p. 92. 


1 A, Hund, High Frequency Measurements (McGraw- 
Hill, 1933), p. 243. 
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capacity, 
dL/L=(dL./L,.\{C/LC+Co]} ~ —2df/f 


if we neglect the effect of distributed capacity 
by limiting C to values which are large with 
respect to C». 

From" N=LI where N=flux through coil 
and J=current in coil, 


dL/L=dN/N—dI/I= —2df/f. 
But 


dN/N=([final flux —initial flux ]/[initial flux ] 
=[((N.+Nm)—(Nit+Nn) ]/(Nit+Nn), 


where NV ,.= flux through copper replacement wire 
at x=0 (Fig. 1)=aH' oP; where H’)=axial field 
= K(I+dI); N;=flux through iron wire = paHoQ; 
where H1)>=KI; N,,=flux around copper wire 
in area (A—a)=[(SA —a)—U(P—1)]H'o; Nx 
= flux around iron wire in area (A —a) =[ (SA —a) 
— U(wQ—1) Jy. The first term of N,, and N, is 
the flux that would be caused in (A —a) by the 
current in the coil. The second term is the 
reduction of flux due to magnetization of the 
sample. The factors S, U, Q, P have been 
defined and will be evaluated later. Substituting 
and simplifying, 


2df /f=dI/I—dN/N 
=((uQ—P)(1+dI/I) }/[(uQ—1)+SA/(a—U)] 


or since dJ/I<1; df/f«1; and 2SA/(a—U)>2; 
u=(1/Q){L2SAdf]/[f(a—U)]+P}. 
Since Q is a function of yu, uw may be evaluated 
by successive approximation. 
IV. EVALUATION OF S 


The flux in the cross section of a coil is obtained 
by integrating the field over the cross section. 
The field is given by Gray” for any point in 
the coil. When x=0 


H=2nnI\2L/\/(R?+L?*) 
+[6R*Ly* ]/[4./(R?+L?)*] 
+[30RL(3R?—4L?)y*] 
[64y/(R?+L*)*}+- =}, 


“L. Page and N. I. Adams, Principles of Electricity 
(D. Van Nostrand Company, 1931), p. 333. 

2A, Gray, Absolute Measurements in Electricity and 
Magnetism (Macmillan, 1927), pp. 210 ff. 
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Fic. 3. Variation of flux in sample wire with frequency. 


R 
N=2n f Hydy =H(AS—aS’) 


=IHT,(AS—a) since S’~1, 
where 


S=14+3R4/[8(R?+L?)?] 
+[5R5(3R?—4L*) ]/[64(R?+L?)*]+---. 
V. EVALUATION OF Q 


Berg" has shown that the differential equation 
0°B/dy?+[1/y jaB/dy=([4rp/p JaB/at 


(where p=resistivity =10' e.m.u. for iron, and 
y=radial distance) gives the condition for radial 
variation of flux density in the iron. We are 
dealing with a sine wave field caused by a sine 
wave current in the coil. At small field strengths 
B= B, cos wt+ Bz sin wt, where the phase angle is 
zero at the surface of the iron and varies radially. 
The effective value also varies radially from the 
surface value Bo. 

The solution of this differential equation is 
exactly analogous to the solution of the similar 
equation for skin effect and yields, using the 
notation employed by Berg," 


B=[Bov/2/(ber? mr+bei? mr) \[ (ber mr ber my 
+ bei mr bei my) cos wt — (bei mr ber my 
— ber mr bei my) sin wt } 
where 
mr =rv/ \ 82 uf /p} =0.78\/ug 


where g=f-10-* and since p=10' and r=0.0088 
cm 


N; (instantaneous) = an f Bydy 
0 
%E. J. Berg, Electrical Engineering, Advanced Course 
(McGraw-Hill, 1916), pp. 151 ff. 
4 Reference 13, pp. 225 ff. 
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and 


N, (effective) =[Bo2ar/m }x/{((ber’ mr)? 
+ (bet’ mr)? \/ ber? mr+ bei? mr |}. 


Using Russell’s approximations'® N;=yH,aQ, 
where for values of mr<2 


O=1-—0.494-10-*(ug)?+0.419-10-4*(ug)4 
— 0.3918 -10~®(ug)®+--- 


and for values of mr>4 


O= 2.585/+/ (ug) —0.581/ (ug) +0.1325// (ug)? 
+0.4175/(ug)?+---. 


These two expressions for Q have been plotted 
(Fig. 3) over the ranges for which they give 
acceptable accuracy and the two curves joined 
by a smooth curve to provide continuity over 
the range where neither approximation is satis- 
factory. Furthermore, since in the derivation 
of these expressions » and p occur only as a 
ratio, the values of P may be read from the 
same graph by using 5.65g instead of yg since 
5.65 is the ratio of the resistivities of iron 
and copper. 


VI. EVALUATION OF U 


Consider a wire of infinite length extending 
along the axis of the coil (Fig. 1). At a point x,'® 


= rf (Bzy—H.)ydy = 4 (uQ— 1), 


0 
= H(wQ—1)[r?/(R?+L?)/(8L) ] 
X | (x+L)//LR?+(«+L)?] 
—(x—L)/V/LR?+(x—L)?]}. 
The surface magnetization in this wire would be 


D’.= —dN,/dx 
= Ho(uQ—1) LP R®vV/(R?+L?*)/(8L) J 
x 11//[R2+ (x—L)? B—1/\/(R2+(x +L)? #}. 


In addition to this magnetization due to the 
non-uniform field, in our sample of finite length 
2/, the flux which would otherwise be in the 
permeable material, beyond x=1 will be “‘re- 
flected” from the ‘‘open” end and leak out 
along the length of the sample in such a manner 


® A, Russell, Phil. Mag. 17, 528 (1909). 
'® Reference 12, pp. 210 ff. 
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that the surface magnetization due to this 
“‘reflected”’ flux is proportional’ to x and to H,; 
i.e., D’,=C\xH, in which C; may be evaluated 


from 
l 


f D" ,dx=Ni; 
0 


(N, for the wire of infinite length). The total 
surface density of magnetization then 


D,=D',+ D",= —dN,/dx+CixH, 
=(uQ—1)Ho-f(x), 
f(x) = C2{1/+/LR?+ (x—L)* B—1/4/LR? 
+ (x+L)? *} +C3{x(x+L)//LR? 
+(x+L)?}—x(x—L)/VCR+(x—L)*]}, 
Co=PR%/(R?+L?)/8L 
and 
C3=CiC2/[r?R?(uQ—1) ]. 


It is apparent that these constants must be 
evaluated for each coil and each sample. It is 
also evident that the surface magnetization will 
be different for each coil and each sample. 
Figure 4 shows the surface magnetization factor 
f(x) for the single sample used in each of the 
several coils. Furthermore f(x) is independent of 
the material and of the frequency, being a 
function of coil and sample dimensions only. 
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Fic. 4. Surface magnetization of sample in one set of coils. 


The reduction of the field of the coil at a 
point (0, y) due to this surface magnetization 


l 


nH =2 f D,xdx//(x?+y?)’ 
0 


(y is not a function of x). With the expression 
for D, as derived above, this integral becomes a 
very complicated one involving the use of 


17]. A. Ewing, Magnetic Induction in Iron and Other 
Materials (Electrician Publishing Company, 1892), pp. 25 ff. 
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TABLE I. Coil data. 


Half 


Coil Radius length 
number Winding cm cm) lurns S U 

| Double layer 1.85 2.855 170 1.031 85.6 
2 Single layer 1.85 2.855 120 1.031 85.6 
3 Single layer 1.85 1.905 75 = 1.084 95.5 
4 Single layer 1.56 2.38 52 1.034 74.8 
5 Single layer 1.56 1.59 33 1.087 81.5 
6 Single layer 1.56 1.19 25 1.149 88.5 
7 Single layer 1.56 1.03 20 =—-1.196 95.5 
8 Single layer 1.56 0.90 17. =: 1.230 99.0 
9 Single layer 1.56 0.45 9 1.395 115.0 

10 Spacedturns 1.61 1.27 17. 1.147 = 108.0 

11 Spacedturns 1.61 1.27 11 61.147 = 108.0 


reduction formulas which give a multiplicity of 
elliptic integrals. Furthermore since subsequent 
to this integration it will be necessary to inte- 
grate with respect to y over the cross section of 
the coil in order to find the flux reduction, the 
complexity of rigorous analysis becomes too 
great. Approximate formulation of D, allows the 
integration to be performed quite readily. 
Approximating f(x) by straight line segments and 
applying the two point formula several times, 
the result is: 


A" =2uQ—1)Ho f [ f(x)xdx//(x?+ y?)* ] 


n 


= 2(uQ—1) Ho is tL finaz) 


n=O 


—fin—1az) Ax } 


nic 
x | [x2dx/4/(x?+y2)*] 
(n—l1)Az 


+[ nf, Naz ("# asi 1 naz | 


nAz 
xf 
(n—1)Az 


n 


[ 2f, nAiz) ~~ (fin— az 
0 


[xdxr/y/(x*+-y*)*]| 


=[2(uQ—1)Ho/Ax]\ 


n 


+finsaz) log [nAx++/(n?Ax?+ y?) }} 
—frazilog [nAx+y ‘(n*Ax?+y?) ] 


+[Ax/s/(n*Ax?+y?) ]}). 


Integrating over the area of the coil 


R 


N”= arf H"' ydy = (nO — 1 TU, 


0 
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U=(24 Ax)(>- 1 2ftnarr— (fin ynartfintnaz) | 


n=0 


‘LR? log (nAx+ \/[n*Ax?+ R?)) 
—4(nAx—/[n?Ax?+ R?])?]}} 
—fraz| R* log (nAx+y/[n*Ax?+ R?*}) 
—4(nAx—/[n?Ax?+ R?])? 

— 2Ax(nAx—\/[n?Ax?+ R*])}). 


It is apparent that U is a function of coil and 
sample dimensions only and is therefore inde- 
pendent of frequency and sample material. It 
has been evaluated for each coil (Table I). 

By evaluating N” for the entire area of the 
coil, this analysis allows for the reduction in 
field inside the sample due to magnetization of 
the sample, since analysis shows it to be of the 
same order as the field reduction outside the 
sample. Actually this internal reduction in flux 
is negligible, due to the small area of the sample. 


VII. CONDITION FOR INITIAL PERMEABILITY 


An upper limit is set to the value of tuning 
capacity by the necessity of keeping the axial 
field of the coil low enough to insure initial 
permeability; i.e., J7<1 oersted. From H 
=0.4rnI cos 6 (Fig. 1) and J=E/wL=Ey(C/L) 


600 
5 6 
Ss 4 Jy 7 
40 -$—4 Ze 


200 — 


Fic. 5. Observed data for one set of coils. 


since w~1/,/(LC) (where E=resonant voltage 
and wl=reactance of coil~impedance) J//» 
=0.42n cos 6Ey\/(C/L) <1 oersted. Oscillograph 
measurements indicate E<5 volts in the oscil- 
lators used. For the longest coil using the finest 
wire this gives C<100L; (C in wyf, L in wh). 
For the shorter coils using larger wire C<900L. 
This allows the use of very high C circuits, 
promoting oscillator stability, and extending the 
usable range of the coil. At low frequencies a 
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multilayer coil is necessary, m increases, and the 
maximum C must come down. There is thus a 
low frequency limit to the application of the 
method (= 50 kc). 


VIII. DISCUSSION OF RESULTS 


The observed data for one set of coils are 
shown in Fig. 5. These values represent the 


‘average of from four to ten individual readings, 


the probable error of which depends on the 
frequency, being greater at higher frequencies 
due to decreased stability of the oscillator. On 
the average the probable error is less than 
2 percent. 

A comparative examination of the results of 
other investigators, listed in the introduction to 
this report, reveals considerable disparity, the 
results being scarcely comparable because of 
differences in sample size and composition, 
difference in frequency range and differences in 
method and analysis. 

In general, at frequencies below 40 mc and 
for the limited range of frequency employed by 
any one investigator, wr as calculated from 
measures of inductance shows an apparent 
variation with frequency. This has been due to 
a use of approximation formulas in the analysis, 
and to the failure to limit the application of 
these formulas to the range of frequency where 
the errors introduced by the approximation 
would be legitimately negligible. The necessity 
for this limitation is made visually apparent in 
Fig. 6 where the dotted portion of each separate 
curve is that part of the curve discarded because 
of the aforementioned limitation. 

The variation of yw, with wave-length is 
shown in Fig. 6. The permeability uz is practi- 
cally constant over the entire radiofrequency 
range from 30 kc to 40 mc. The apparent increase 
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at the high frequency end of the range of each 
coil is due to the neglect of the distributed 
capacity factor which is legitimately negligible 
at lower frequencies. 

In the range between 50 kc and 30 ke the 
apparent increase in wu, with decreasing frequency 





TTTls h lv TT LUM 
(00 8H ft oe ee eee 
8 \ | 

60 +t -+—_# sth t Jiit + 
wm He ~—y TOT} ~ Pitt ++4++4 
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Fic. 6. Variation of ux with wave-length X. 


is due, in part at least, to an increasing axial 
field strength, since even with a multi-layer coil 
the tuning capacity was very large for these 
low frequencies. 

Since coils 10 and 11 were space wound, they 
had very little distributed capacity, allowing a 
check on the variation of wu, to higher frequency 
than was possible with a close wound coil. 
However, the theoretical development of the 
values of S and U was'based on the theory for 
a close wound coil. This explains the disparity 
between the results obtained with coils 10 and 
11 and those obtained with the other coils. 
The results with these coils are therefore quali- 
tative only. 

The difference in the results from coil to coil 
is due to departures of the close wound coil from 
the ideal current sheet and to uncertainty in 
the approximate evaluation of U. 

I wish to express my sincere appreciation to 
Dr. J. Barton Hoag for suggesting the problem, 
and for his encouragement and kindly advice 
during the progress of the investigation. 
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Demagnetizing Factors of Rods 


R. M. BozortH AND D. M. CHAPIN 


Bell Telephone Laboratories, New York, New York 


(Received March 4, 1942) 





A chart is constructed for converting the apparent magnetic permeability to the true permea- 
bility of cylinders of any given ratio of length to diameter. The curves are based on previous 
calculations in which account was taken of the variation of the demagnetizing factor N with 
permeability, but in which the permeability was assumed constant over any one cylinder. The 
flux distribution has been determined experimentally in several cylinders, and as the field acting 
on the specimen is increased from zero, the positions of the effective poles have been found to 
move toward the middle of the rod until (if the permeability is sufficiently high) they are 
located about 0.7 of the distance from the middle to the ends, and then to move toward and 
approach the ends as limits as the permeability declines in high fields. 


HEN using the ballistic galvanometer for 

making magnetic measurements on cylin- 
ders or rods it is common practice to determine 
directly the ‘“apparent”’ permeability yu’, equal 
to the quotient of the flux density B and the 
applied field strength I/’. The B is determined 
at the middle of the rod and II’ is calculated 
simply from the current flowing through the 
magnetizing solenoid. Then yp’ is converted to 
the (true) permeability » by calculation using 
the appropriate demagnetizing factor N, thus 
taking account of the field strength produced at 
the middle of the rod by the rod itself. One 
purpose of this note is to provide a simple 
graphical means of converting true to apparent 
permeability or vice versa when the specimens 
are cylindrical rods of any ratio of length to 
diameter. The chart shown has already found 
considerable use in practical problems dealing 
with such rods. 

Other purposes of this report are to record 
experiments on the distribution of magnetization 
in several rods, to show how the positions of the 
poles change with permeability and to compare 
these findings with theory. 


GRAPHICAL CONVERSION OF wu TO w’ 
The demagnetizing factor N is given by 
H=H'—NI 
or the equivalent 
HT =H —(N/4r)(B-—T), 


in which // is the (true) field strength, J/’ the 
applied field strength equal to that which would 


320 


exist if the specimen were removed, J the 
intensity of magnetization and B the flux density 
measured at the middle of the rod. The N 
considered here may be termed the ballistic 
demagnetizing factor to distinguish it from the 
magnetometric demagnetizing factor which is 
not multiplied by the J at the middle but by 
some kind of an average of the intensity of 
magnetization over the whole bar depending on 
the position of the magnetometer needle. 

The equation given above may be rewritten: 


1 1 N 


u—1 p’—p'/p 4n 

1 N 1 WN 

I~ )- mn 

m An uw 6 Are 
the apparent permeability being defined as 
u'=B/IHI’. In practice sufficient accuracy is 


or 


always obtained by using the simpler equation: 
1/u=(1/u')—(N/4n). 


It is thus possible to plot wu against uw’ for a 
cylindrical rod of given dimensional ratio pro- 
vided N is known, even if N varies with per- 
meability. Such a plot is shown in Fig. 1. The 
values of N have been taken from the compila- 
tions of Stablein, Schlechtweg, Neumann, and 
Warmuth! and are the results of calculations 
_'H. Neumann, Wiss. Veréff. Siemens-Konz. 10, 55-71, 
No. 2 (1931); H. Neumann and K. Warmuth, Wiss. 
Ver6ff. Siemens-Konz. 11, 25-35, No. 2 (1930); F. Stablein 
and H. Schlechtweg, Zeits. f. Physik 95, 630-46 (1935); 


Kk. Warmuth, Arch. f. Elektrotech. 30, 761-79 (1936); 31, 
124-30 (1937) and 33, 747-63 (1939). 
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that agree within experimental error with the 
experimental determinations. The curves of Fig. 
1 are corrected for variation of N with yz, the 
permeability being that corresponding to the B 
existing at the middle of the rod. The original 
calculations of N were based on the assumption 
that uw is constant over the rod and this may in 
certain circumstances cause a definite error in 
the calculated yu vs. pv’ relation. 

The values of N/42 appropriate for .= ~ are 
plotted, in the upper left-hand portion of Fig. 1, 





against the dimensional ratio, m=length/di- 
ameter, of the rod. Since np’ =42/N when p= ~, 
the scale-reading at the right (for yw’) is the 
reciprocal of the scale-reading at the left (for 
N/4r) for any given horizontal line (ordinate). 
Also, the asymptote for any given y vs. p’ curve 
(m=constant) lies on the same horizontal line 
as the value of N/4m that corresponds to this 
value of m. 

The method of using the chart is obvious. 
Interpolation for values of m not given may be 
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: Fic. 1. Relation between apparent permeability uv’ and true permeability » of cylinders having any given ratio m of 
‘ length to diameter. Demagnetizing factors N as defined by H = H’— NT, are shown as a function of m for rods having very 
high permeabilities. The curves may also be used for converting y’ to u for ring specimens in which there is an air gap 
subtending an angle a@ at the center of the ring. 
S 
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Fic. 2. Demagnetizing factors of ellipsoids and cylinders (see the equations and reference 1). 


performed as follows: The appropriate value of 
m (say, 30) is located on the scale at the top 
and the corresponding value of N/42 (0.0027) 
found on the curve of N/4x vs. m. The value of 
N/4rm is noted also for a nearby value of m for 
which a yu vs. pw’ curve is drawn (e.g., m=20). 
Horizontal lines are drawn through the points 
corresponding to these values of m and N/4r, 
and the 45° distance x between these lines is 
determined. The line x of the same length and 
45° inclination, may then be drawn from any 
point on the curve for m=20 and will give a 
point on the curve for m=30, as shown by the 
lines on the chart. This 
method of interpolation is based on the fact 


dotted construction 


that all of the yu vs. uw’ curves are similar, except 
as influenced by the variation of N with yw’, and 
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differ only by a translation along a line inclined 
45° to the axes. The amount of such a displace- 
ment at w=~* is obviously the same, on a 
logarithmic scale, as the displacement made 
when a point on the N/4z vs. m curve is moved 
from one value of m to another; therefore, this 
is the displacement for all parts of the u vs. p’ 
curve. The slight correction in the interpolation 
for variation of N with zis at all times negligible. 

For convenience the demagnetizing factors for 
oblate and prolate ellipsoids of revolution, and 
for rods having finite permeabilities, are plotted 
in Fig. 2. Convenient forms of the equations for 
accurate calculation of N/4m for the ellipsoids 
are as follows: 


(1) For oblate ellipsoid with axes a=c/m, 
b=c, 
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(a) Magnetized parallel to a long axis, b or 


N 1 m: _ (m?—1)! 1 
= are sin | 
Aner 2L(m*?—1)! m m*— 1 


N xwm-—2 
= if 
4r 4m? 


~ 


m> 1. 


(b) Magnetized parallel to the short axis, a, 


N m- m? (m?—1)3 
= - are sin ; 
4er m*—1 = (m*—1): m 
N T 
=1- if m>1. 
An 2m 


(2) For prolate ellipsoid with axes a=mc, 
b=c, 


(a) Magnetized parallel to the long axis, a, 


N 1 | m m+(m?—1)} 
= In —1}, 
dr m?—1\2(m?—1)3 m—(m?—1)}! | 


N 1 
=—(In2m—1) if 
4r m°* 


m> 1. 


(b) Magnetized parallel to a short axis, 6 or c, 
N 1 m* m m+(m*—1 


) 
= — In | 
4r 2\m?—1 (m*—1): m—(m?—1)44 


N 1 In 2m ; 
am (1 — ) if m>1. 
4n 2 m* 


Neumann and Warmuth! give the following 


formula for calculating N/42 for long cylindrical 
rods having infinite permeability: 


(N/42r), + = (4.02 logio m—0.92)/(2m?) m=10. 


It should be remembered that .the demagni- 
tizing factors, Ni, Ne, and N3, in three mutually 
perpendicular directions in a specimen of any 
shape add up to 47: 


N, a Not N3 = dr. 


From this it follows that the curves of Fig. 2 
approach N/4r=1 as an 
approaches zero. Also it is convenient to remem- 
ber that for a thin ring with an air gap that 


asymptote as m™ 


subtends the angle @ at the center of the ring, 
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the demagnetizing factor is given simply by 
N/42r=a/360 


if @ is in degrees. Thus for each rod of dimen- 
sional ratio m, there is a value of N that applies 
also to a ring having an air gap of angular 
width a. In Fig. 1 these values of a are recorded 
in parentheses just below the corresponding 
values of m, and the figure may thus be used 
for ring specimens as well as rods. 


FLUX DISTRIBUTION 


In addition to a knowledge of the true per- 
meability of a rod, as determined from H’ and 
from B as measured in a short search coil 
placed at the middle of a bar, it is often desired 
to determine this » when the search coil covers a 
considerable portion of the rod or even is longer 
than the rod. This can be done if the distribution 
of flux in the rod is known. Flux distribution 
was determined in two rods having widely 
different permeabilities but the same diameter 
(0.97 cm), length (50.6 cm), and dimensional 
ratio (52.4). Both rods were cut from the same 
stock of molybdenum permalloy containing 79 
percent nickel, 4 percent molybdenum, and the 
rest iron; one was heat-treated for high per- 
meability while the other was unannealed, as 
swaged. 

Measurements were made with a fluxmeter 
using a small search coil 2 mm long placed at 
various positions along the rod. The applied 
field of given strength HJ’ was reversed in 
direction to produce the required change in flux. 
The results are shown in Figs. 3 and 4 in which 
the flux density B; at any point along the rod 
is referred to the density B at the middle. 
Actually flux-densities are corrected by sub- 
tracting JJ’ from each of them; they are in 
reality B,—H’ and B—H’. The abscissa //lo is 
the position of the coil on the rod expressed as 
the fraction of the distance from the middle to 
the end. Each curve is marked with the corre- 
sponding value of H’ and corresponding values 
of uw’ and B are given in the insert. 

The distribution of flux in the annealed rod 
when the highest (H7’=5.0, 
uw’ =871) is nearly parabolic. For comparison, 
the dotted line of Fig. 3 is drawn as a parabola 


permeability is 
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cutting the experimental curve at B,;/B=0.5. 
The nearby dashed line is the distribution 
accurately calculated by Wiirschmidt? under the 
assumption that u.= ©; it is nearly parabolic. 
The position of the effective pole, given by 


*B,—H' | 
Iy/lo= f d ’ 
9 B-H’ ly 
may be determined* for each of these theoretical 


curves and for the experimental curves; for the 
parabolic distribution it lies at /,/lo=0.73, 


2]. Wiirschmidt, Theorie des Entmagnetizierung Factors 
und der Scherung von Magnetizierungskurven (Vieweg, 
Braunschweig, 1925), 

* This integral is equal to 


“Bi-H ,l . , 
J, BoH’ ‘i, since I, (Hi—H’')d(1/lo) =0 


- [. (Bi—H) d/o) =0, 


as pointed out to us by Mr. H. J. Williams. In all of our 
experiments B—H and B—H’ (both at the middle of the 
rod) are indistinguishable. 





Wiirschmidt’s distribution gives 0.73, and for 
these experimental curves /,/lo lies between 
0.72 and 0.96. 

The demagnetizing factor of a rod depends 
not only on the dimensional ratio and on the 
permeability, if that is uniform throughout the 
bar, but also on the variation of permeability if, 
as is generally the case, it is different in one 
portion of the rod than in another. It has been 
reported* that the positions of the effective poles 
move toward the middle of the rod and then 
outward again as the magnetization is increased 
from zero, and it has been pointed out* that 
they are nearest the middle when the magnetiza- 
tion is approximately that corresponding to 
maximum permeability. The more recent calcu- 


3L. Holborn, Sitz, Preuss. Akad. Wiss. 159-68 (1898): 
C. G. Lamb, Phil. Mag. 48, 262-71 (1899). 

*G. F. C. Searle and T. G. Bedford, Phil. Trans. Roy. 
Soc. A198, 33-104 (1902). For further discussion of change 
of N with B, see E. Dussler, Ann. d. Physik 86, 66-94 
(1928) and D. Foster, Phil. Mag. 8, 304-13 (1929). 
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Fic. 3. Distribution of flux in rod of annealed permalloy having dimensional ratio m=52.4. Values of B, H’, w’ and 1,/lo 
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(pole position) are indicated in the table. Theoretical distribution (broken line) for «= ~ is as calculated by Wiirschmidt. 
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Fic. 4. Distribution of flux in rod of unannealed permalloy (m=52.4). B, H’, wu’ and 1,/lo are as indicated. Wiirschmidt’s 
curve (broken line) and an experimental curve for the annealed rod (/7=20) are shown for comparison. 


lations of Wiirschmidt establish the fact that a 
high demagnetizing factor and a small interpolar 
distance are associated with high permeability. 
Thus, theory and experiment justify the follow- 
ing description of the movements of the effective 
poles as the applied field is increased from zero. 
The poles appear first at some position, between 
1,/lyb=0.7 and 1,/lp=1, determined by the true 
initial permeability uo of the rod. As HT’ increases 
uw also increases in all ferromagnetic materials 
and the poles, therefore, move toward the middle 
of the bar, approaching the limiting position 
1/ly4+0.7 the more nearly the higher the maxi- 
mum permeability, um. As J/’ increases further 
and y» declines from yu, to 1, the poles move 
toward the ends of the rod and approach these 
ends as limits. The positions of the poles are 
given in the table inserted in Figs. 3 and 4. 
When the permeability is well below the 
highest value the distribution is far from para- 
bolic, as is evident from the curves. In the lowest 
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fields the permeability of the unannealed rod 
is low and uniform, and it is possible to calculate 
the flux distribution and compare with experi- 
ment. Such calculations could be made as 
described by Stablein and Schlechtweg! but they 
are laborious and have not been carried out 
for this report. When the applied field-strength 
is increased the permeability at the middle is 
greater than that at the ends; the poles are then 
nearer the middle than they would be if the 
permeability of the whole bar were that of the 
middle, for the low permeability at the ends 
effectively shortens the rod. That this is so is 
indicated by the data for the annealed rod for 
HI’ =0.5; in this case the poles were already near 
their closest-in positions because of the high 
permeability, and the increase in u» at the middle 
causes the poles to move in still farther to 
l,,/ly=0.72, a minimum value for this rod and a 
value less than that corresponding to Wiir- 
schmidt’s calculated distribution for p= ~. 
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When either rod is subjected to an intense 
fieid, the distribution near the middle is naturally 
quite uniform because the field over a consider- 
able portion of the rod is sufficient to cause 
saturation. The distribution curve for H’=20 
for the annealed rod may be contrasted with 
that for H’=0.5 for the unannealed rod. The 
former curve shows only a very slight dropping 
of the flux in going from the center to //l)=0.4, 
and this is due to the fact that this rod is 
practically saturated while the unannealed rod 
has an almost uniform permeability over its 
entire length. 

Attention may be called to the very slight 
differences in distribution in the unannealed rod 
when H’=0.50 and 33.2. This may be attributed 
to the relatively slight change in permeability 
in different parts of the rod, as compared to the 
large range in permeabilities in the annealed rod. 

It seemed desirable to establish for rods of all 
dimensional ratios the apparent fact that the 
effective pole has a minimum interpolar distance 
given by /,/lo>=0.7. To this end a wire of 70 
permalloy was heat treated in a magnetic field® 
to develop high permeability and its effective 
pole positions were measured. The wire was 
30.4 cm long, 1.01 mm in diameter (m=300), 
its effective permeability yu’ =18,250,* and the 
pole positions were found to be /,/l)=0.70. 
Thus, the minimum pole separation correspond- 
ing to /,/l)>=0.7 may be assumed to apply to all 
cylinders with very high permeability. A rod 
may be considered to have ‘‘very high perme- 
ability”’ for this purpose, if its uw is so high that 
increasing it does not increase the yp’ of the rod 
by more than 10 or 20 percent as indicated by 
the curves of Fig. 1. 

5]. F. Dillinger and R. M. Bozorth, Physics 6, 279-84 
(1935). 

* Using Neumann and Warmuth’s Eq. (12) of reference 
1, N/4xr=0.00005021 for p= ~ and m=300. Combining 
this with our wp’ =18,250 gives »=218,000, a value con- 
sistent with previous results for this material, reference 5. 
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LENGTH OF SEARCH COIL 


The distribution curves can now be used to 
calculate the way in which the flux threading the 
search coil varies with the coil’s length. In 
practical measurement this amounts to being 
able to determine B at the middle of the bar 
from the galvanometer deflection made with a 


‘search coil of any length. The calculation jis 


simply the averaging of B, over the length of 
the search coil, and results are shown in Fig. 5 
e 2 


cose 


Fic. 5. Effect of length of search coil on measurement of 
B. Search coil (of length 2/) and rod (of length 2/5) are 


concentric. B is flux density at middle, B; is the average flux 
density in that part of the rod that is “covered” by the 
search coil. The dashed line is calculated for the parabolic 
distribution of Fig. 3. 


for the two molybdenum permalloy rods in 
fields of IJ’=0.5. Search coils of two lengths, 
one 88 percent of the length of the rod and one 
considerably longer than the rod, were con- 
structed and points on the curves were checked 
as indicated in the figure. The dotted line is 
calculated for the parabolic distribution curve 
shown in Fig. 3. It is evident that for a rod of 
dimensional ratio m=50, an effective 
permeability of u’ =700 (u=3000) is sufficiently 
low to cause an appreciable departure from the 
curve for the 


near 


nearly distribution 
associated with infinite permeability. For the 


upper curve of Fig. 5, uw’ =74 and /,/1)=0.86. 


parabolic 
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The Electrical Oscillations of a Perfectly Conducting Prolate Spheroid* 


RoBert M. RyDERt 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
(Received January 20, 1942) 


The forced oe a perfectly-conducting prolate spheroid of eccentricity nearly 


unity are shown to be 


»ynposable into “‘harmonics”’ corresponding to different modes of 


vibration, each harmonic beifig-quantitatively connected with a certain portion of the im- 
pressed electric field which drives the antenna. The harmonics contribute additively to the 
current and field of the spheroid; each offers a characteristic impedance to the driving field, 
and the properties of the antenna are a composite depending upon the proportions of the 
various harmonics present. The behavior of the harmonics with frequency is discussed qualita- 
tatively; analytical expreésions obtained are useful chiefly at the resonant frequencies of the 
antenna, where the mos##i#nportant harmonic becomes sinusoidal in character. 


s 


HE partial differential equation of the 

electromagnetic field, as is well known, 
can be separated in the coordinates of various 
simple geometrical surfaces. In this way have 
been studied the fields of a prolate spheroid by 
Abraham,' of a paraboloid by the same author,’ 
of an infinite cylinder by Sommerfeld,’ of two 
parallel infinite cylinders by Mie,‘ of the sphere 
by Mie’ and Debye.* Necessarily limited to 
rather simple surfaces (of which the prolate 
spheroid seems best to approximate a finite 
straight wire), the method has compensating 
advantages: perhaps the most important is the 
physical insight obtained into the problem. 

The diameter of a prolate spheroid falls to 
half its maximum value at a distance 0.866a 
from the center, where a is the semi-major axis; 
hence the “equivalent spheroid’ is slightly 
thicker at the center and somewhat longer than 
the cylinder it represents.7* From this rough 
idea of equivalence one might expect that 
measurements on conductors of uniform cross 
section would give resonant wave-lengths some- 
what longer than those calculated by Abraham 
for the spheroid. ‘“‘Exact”’ calculations for the 

_* Condensed from a dissertation presented to the Faculty 
of the Graduate School of Yale University in candidacy for 
the degree of Doctor of Philosophy, 1940. 

+ Now at Bell Telephone Laboratories, New York, New 
York. 

1M. Abraham, Ann. d. Physik 66, 435 (1898). 

*M. Abraham, Ann. d. Physik 2, 32 (1900), 

_* A, Sommerfeld, in Franck and von Mises Die Differen- 
tial- und Integralgleichungen . . . (1925), Volume 2, p. 526. 

'G. Mie, Ann. d. Physik 2, 201 (1900). 

®>G. Mie, Ann. d. Physik 25, 377 (1908). 

’ P. Debye, Ann. d. Physik 30, 56 (1909). 

7 Lord Rayleigh, Phil. Mag. 8, 105 (1904). 

> L. Page and N. I. Adams, Phys. Rev. 53, 819 (1938). 
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free oscillations of a cylinder have been made by 
Oseen® and extended by Hallén,'® whose results 
hold also for curved or even closed coil antennas. 
Like all the above-mentioned authors Hallén’s 
paper discussed only natural frequencies and loga- 
rithmic decrements of free oscillations, leaving 
the steady state alone; still, the experiments of 
Englund" furnish encouraging confirmation. 
Two years later the question was reopened by 
Lindman” who in new measurements 
obtained the very curious result that the 
formulas of Abraham derived for the prolate 
spheroid gave better correspondence than those 
of Hallén, presumably derived for the cylinders 
which he was using. Subsequent discussion™ 
leaves the impression that the mathematical 
approximations of the various methods may 


some 


introduce as much error as the physical approxi- 
mation involved in changing the shape. Accord- 
ingly the forced oscillations of a prolate spheroid 
would appear to be of some practical interest. 

In the extensive literature of the antenna 
there has more often been used the method of 
calculating the field from some assumed current 
distribution in the antenna. Many shapes of 
antenna can be treated; perhaps typical is 
G. W. Pierce’s work on the L-shaped antenna’ 

9C. W. Oseen, Archiv. f. Matem. Astr. och Fysik 9, 
Nr. 12 (1913); Nr. 13 (1914). 

10 E. Hallén, Uppsala Univ. Arsskr., Matem. o. Naturv. 
Nr. 1 (1930). 

1 C, R. Englund, Bell Sys. Tech. J. 7, 404 (1928). 

2K. F, Lindman, Ann. d. Physik 13, 358 (1932). 

13 E, Hallén, Ann. d. Physik 14, 586 (1932). 

4K. F. Lindman, Ann. d. Physik 15, 127 (1932). 

15 G. W. Pierce, Electric Oscillations and Electric Waves 
(1910). 
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which also illustrates the great mathematical 
complexities involved. A discussion of the 
validity of these methods is given by Grosskopf.'* 
From the point of view of the present paper 
they are open to the criticism that usually the 
boundary conditions involved are only approxi- 
mately satisfied, that no quantitative connection 
is obtained between antenna current and driving 
field, and that no explanation appears of why 
sinusoidal current should be a good assumption 
in some cases but not in others. 

Recently L. V. King!’ has developed another 
method with the unfortunate property that the 
fundamental integral equation requires re-solu- 
tion for each new way of driving the antenna. 
More recently still Page and Adams* have 
studied oscillations forced in a perfectly con- 
ducting prolate spheroid by a uniform electric 
field parallel to its axis and varying sinusoidally 
with the time—the problem of the receiving 
antenna. The present paper may be regarded as 
a continuation of the study of oscillations forced 
by a general electric field of that kind. The 
same problem has been studied by Stratton and 
Chu,'* who used rather different 
their analytical the 
equation. 


methods in 


attack on differential 


I. THE DIFFERENTIAL EQUATION 


After the notation of Page and Adams,* if 
xyz are a right-handed set of rectangular co- 
ordinates with the x axis along the long axis of 
the prolate spheroid, and if p?= y?+2*, then may 
be introduced the prolate spheroidal coordinates 
£, n, ¢, forming a right-handed set in the order 
named, and given by the relations 


x = fén, (1) 
p=fi(1—#)(n?—1)}}, (2) 
¢=tan~! (z/y), (3) 
where —1<£<1, 1K9<«%, OS G@<2n. The 


constant f is the semi-interfocal distance of the 

common foci of all the prolate spheroids 7 

=constant= 7. To the remarks given by Page 
16 J. Grosskopf, Hoch. u. Elek: akus. 49, 205 (1937). 
17. V. King, Phil. Trans. Roy. Soc. A236, 381 (1937). 


18. J. Chu and J. A. Stratton, J. App. Phys. 12, 241 
(1941). 
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and Adams on the properties of these coordinates 
we need add only that for very thin spheroids 
(no—1<1), € represents distance measured along 
the spheroid from its center in units such that 
the total length is 2. Also, readers are warned 
that the definitions of € and » are interchanged 
from those used by Stratton and Chu. 

We are interested in fields whose time de- 
pendence is specified by the factor exp (—iw/), 
and in axially symmetric waves, so that the 
only non-vanishing field components are 


fb i] 1 OA 
Ex=i( ) ; (4) 
KJ ef \(1—£&*)(n?—£&) |? 6 


;? On 


; [ph i] | OA 
E,=-1 ) : , (5) 
KJ ef \(n*?—£&)(n?—1)}? 0€ 


S 


1 1 
Hh= A, (6) 
f $(p?—1)(1—-#2) }3 


where A=pll/, satisfies the differential equation 
0°A 


(1—£*)——+ (n?-1) 
ag? On? 


a°A 
+e*(n°—&)A=0 (7) 


as was shown by Abraham.' Heaviside-Lorentz 
symmetrical units are used throughout. Here «x 
and uw are, respectively, the permittivity and 
permeability of the medium surrounding the 
spheroid, while 


€=(w v)f=2nr(f/dr), (8) 


where v=c/(xu)* is the normal wave velocity in 
the medium, c is the velocity of light in vacuum, 
and X is the wave-length. Noteworthy is the 
fact that € is proportional to the frequency v of 
oscillation of the field. Separation of the variables 
in (7) is effected by putting A(£, 7) = X(£) Y(). 


1-#)@X /d?+(a-—€?)X =0, (9) 


(n*—1)d*V/dn?+(—at+en’) Y=0. (10) 


The constant of separation a may have only 
certain eigenvalues such that the electric field is 
continuous on the x axis. Equivalent are the two 


eigenvalue conditions 
E:=0 when 


f= +1, (11) 


X(1)=X(—1)=0. (11a) 
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Fic. 1. Qualitative distribution of current in the ninth harmonic: (a) low frequency; (b) resonant 
frequency; (c) and (d) above resonance. 


Since by hypothesis the spheroid is a _ perfect 
conductor, at its surface holds the boundary 
condition 


E:+E=0, n=np, (12) 


in which E£ is an external exciting or driving field. 

There is also a condition on Eq. (10), in that 
at large distances, in order to represent a 
divergent wave, 


Y(n)—e*" as nx. (13) 


II. QUALITATIVE FEATURES OF THE SOLUTIONS 
Both Eqs. (9) and (10) are of the same form: 
u(x) +f(x)u(x) =0. (14) 


From well-known theorems about the qualitative 
behavior of continuous solutions of such an 
equation" it is possible to discover some proper- 
ties not readily obvious from their series develop- 
ments. For example, take X,(£), the pth eigen- 
solution of (9): it is shown below that X,(£) 
gives the current distribution in the antenna 
for the pth mode of vibration. X,(&) is an 
oscillatory function with +1 nodes, counting 
the two at the ends of the antenna; the other 
p—1 nodes are all confined to the interval 
—(a'/e) <<§E<+(a'/e). At zero frequency (e€=0) 
X ,(£)(1—#)~ is proportional to the associated 
Legendre polynomial P,,(£), and a,=p(p+1). 

19 See for example Franck and v. Mises, Die Differential- 
und Integralgleichungen .. . For an acquaintance with 


these theorems the author is indebted to Professor Einar 
Hille of Yale University. 
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As the frequency parameter ¢€ increases, so does 
a, but more slowly, so that when e=p7/2, 
a,=e. Meanwhile the current loops at the ends 
of the antenna have increased until they equal 
all the other loops in size, and at the “‘resonance 
point” «= px/2, X,(&) is a sinusoid. 

At frequencies higher than resonance the 
character of X, becomes radically different: a 
becomes and remains less than &; the curvature 
of the graph of X,(£) against — is away from 
the é axis for | £| >a’/e, and the loops and nodes 
move in toward the center of the antenna. It is 
evident that it may be difficult to represent X, 
by a simple approximation valid over any large 
frequency range, but near resonance a sinusoid 
may be useful. Qualitative sketches of X9(&) 
(showing only half the spheroid) appear in 
Fig. 1. 

From the Y equation in the form 


9 


i ‘ Yp 9 4 
n?—1 


it is clear that for large n, regardless of p, any 
solution approaches the value Ae‘’+Be-™, 
where A and B are constants. The notation Y 
is reserved for that solution which is asymp- 
totically 


Y,= U,+iV»—cos ent+isin en as n>. (15) 


Since the coefficients in the equation are real, 
both U(») and V(m) satisfy (10a) separately; 
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clearly also they are independent solutions. For imaginary, and V,(1) has always the same sign 


p odd, U, is proportional to X,(n) and is there- as 7?~'. For e<pm/2, | V,(1)|>1; for e=pr/2, 
fore zero at n=1; hence Y,(1)=iV,(1) is pure’ | V,(1)| =1; while for e>pz/2, | V,(1)| <1. 
Ill. ANALYTIC SOLUTIONS OF THE EQUATIONS 
1. Solutions in Integral Form 
The following methods are due to Abraham.'? Let either (9) or (10) be written in the form 


y(v, 


y) =(1—y")o"’(y) + (a— &y*)o(y) =0. (16) 


Here v will denote the solution which in the interval —1< y<1 is holomorphic, and vanishes to the 
first order at the ends of the interval. Now is introduced the transform 


w(x) v2 
= exp (texy)u(y)dy. (17) 
Yi 


1—x? 


Substitution of (17) into (16) and proper determination of the integration limits of the former show 
the following two special forms of w are solutions of (16): 


l 
wily) =(1 -y) exp (zeyt)v(t)dt, (18) 
al 


1+ix 
we(y)=(1 -y) f exp (ieyt)v(t)dt, y>1. (19) 
l 
w,(y) vanishes for y= +1; hence except for a multiplicative constant it must be equal to v(y). If we 
restore the (X, &) notation, 
l 


e,X »(&)=(1 -#)f exp (1e&t) X ,(t)dt, (18a) 
1 


where e, is a function both of the eigenvalue index p and of the frequency parameter e. 
(19) holds only for y>1; we suspect that it may have something to do with Y(). This is indeed 
the case: integration by parts, followed by a few transformations, shows that 


e ——— + [ie(n—s)]_ 


Wo(n) =e'+ X ,(s)ds, (20) 
X ,' (1) 16€,X » (1) ' q—S 
from which is concluded, after Abraham, that 
e€- 
wo(n)=Y,(n), n>1 (19a) 
X,'(1) 


is asymptotic to e**" for large 7; hence it is the solution sought, representing a divergent wave. Further 
reduction gives forms suitable for computation : 


(—1)?t!y 7X ,(n) n+1 . 
Y,(n) =e'"+ : log —G,(n) I, (21) 
1€€y X » (1) n—1 


Y,' (n) =tee*"+ og —G,'(n) 3 ( 


a n+1 2 X,(n) 
X ,' (1) n—1 »?—1X,'(1) 


Ne 
i) 


1€e» 
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where 


ds. (23) 


; ein | X,(n) exp (—ten) —X,(s) exp (—tes) 
G,(9)= f 


r Ff 
X,'(1) 
The point of the manipulations is now clear: Abraham has neatly separated out the parts of the 


expression which give trouble by becoming large when n—1. The logarithmic term in Y,’ is the only 
one not finite at 7=1. For long, thin spheroids values of » near 1 are of chief interest ; and to good 


= es ; 


approximation may be written 


(—1)?tly, 
Y,(n) =e*+ (—G,(1)] if n—-1<1, (21a) 
1€€ » 
c (—1)?t'y, n+1 P 
Y, (n) =tee*+ to -1-6,'()] if »—1<1. (22a) 
1€€ » 7—1 


Since G, and G,’ depend only upon X ,, it is seen that (21) and (22) give the desired solution of (10) 
in terms of the desired solution of (9), which is found below by a perturbation method. 


2. Series Solution of (9) 
Y \_ 
X @+(e+ )x,@=0 (9a) 
1—¢ 


The value y,=0 gives the most important solution of (9a), that is, the sine wave which applies at 
resonance. It seems natural, therefore, to try to solve (9a) as a perturbation of the cosine equation, 
i.e., to find X ,(£) as a trigonometric Fourier series. However, even though y, is small near resonance, 
we cannot well regard this equation as a perturbation of X’’+¢X =0, for the solutions of the latter 
do not satisfy the eigenvalue condition (11a). Instead, we write 


, pr\* Y> 1... 
E r+( ) X »|+| + js p=0, (24) 
2 1—¢ 


wherein 6, = €@ — } p*x* (25) and the second bracket is to be regarded as in the nature of a perturbation. 
Application of the usual methods develops the following expansions: 


Yo=Yp19pHVp29p? + Vp3Op? + * °°; (26) 
X p=X pot X piOptX p29"? + °-:, (27) 
=>. a, sin }ka(1—£), pt+k even, (27a) 
=>. do; ap4.0,/ sin $ka(1—E), p+k even. (27b) 


From use of the above come the following expressions for parameters of the integral formulas (21-23) : 


4 > (1 k)a pk 


e,= , k, podd, (280) 
mw > (—1)%-? ay, 
Sie Dox (1/k) aps 
ep=— , k, peven, (28e) 
1” ar (—1)/ "Rd pi 
1 ‘ 
G,(1)= —_— ; A piLi(€), (29) 
X,'(1) & 
1 
G,'(1)=—1+ > aprM;.(€). (30) 
Xp (1) + 
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In (29) and (30) the new symbols L and M are functions of ¢ only, defined as follows: 


* exp (tet) sin ake 
Lil) J dt, (31) 
0 t 


ke ¢? 1—exp (iet) cos $kat 
M;,(¢)= | dt. (32) 


2 Ho t 


Values for the first few coefficients for (26) and (27) are obtained without too much difficulty from 
the subjoined formulas: 


Yp1= —2/F(2p), 


(26-1) 
8 1 \F(p+k)—F(p—k)}? pt+k even, 
g= — ’ (26-2) 
mw? [ F(2p) | & p* —k? pk 
2 
You™ ae adyx| F(p+k) — Fi (p- k)|- SF (2p) y p2” (26-3) 
[F(2p) }* * 
1 if p=k 
a pok = : (27-0) 
0 if p#k, 
F(p+k)—F(p—k) 
——— , if pt+k even, 
Apik =) T° > F(2p) pb? —k? (27-1) 
a if p+k odd or p=k, 
a 1 1 
Apr = —Apkt > Apim| F(m+k) — F(m—k)| — byp2| F(p+k) — F(p—k)}. (27-2) 
r: (p?— k?) F(2p) m 


In looking over these formulas, one observes that a,j, is zero unless the indices p and k are both 
even or both odd: this condition is required by the symmetry properties of the function X 
discussed above. Secondly, unfortunately already in the first order all the sines of the same symmetry 
as X »o have been introduced. For | p—k| large, however, their coefficients are small. The summations 
indicated in the formulas for the coefficients are infinite series whose general terms are not known 
though in some cases estimates may be found. In practice the finding of higher order terms than 
those given above involves prohibitive difficulty, so that the analytic solutions found are useful 
only in the region where the power series for the various functions of 6, are satisfactorily approxi- 
mated by their first three terms. Near resonance it may be expected that the resonant harmonic 
will be accurately determined: to the other modes we shall have somewhat rough approximations 
but ordinarily the other modes are merely small correction terms. 


IV. PROPERTIES OF THE FORCED OSCILLATIONS 
Let the spheroid be acted upon by an exciting field 
E=i8&(é)e~*** (33) 
parallel to its axis. The boundary condition (12) at the surface of the perfect conductor becomes 
(0A /dn)n=n, =1(k/p)*efno(1 — &*) E(E). (34) 
Now any wave which is axially symmetrical can be expanded in a uniformly convergent series of the 
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—— ie 


eigenfunctions of (7); this particular exciting field is defined by 


ee 


1) A =) p CpX p(€) Y,(n), (35) 
(0A /dn)o = ea CpX »(€) Y,' (no), (36) 
2) wherein the c’s may be formally evaluated by substitution of (36) into (34), multiplication by 
(1—#)"X,,(€), and use of the fact readily seen from (9a) that the X,’s form an orthogonal set over 
the interval —1<&<1 with the weight-factor (1—)-. If N(p) symbolizes the normalizing factor 

m 


of the orthogonal set, then the result is 





K\? : 
1) i( ) af &(E)X ,(E)déE 
) rn 1 (37) 
| t_* r ‘ ’ 
2) N(p) Y>' (no) 
where ; ; 
i UX »(E)X p(E) 0, nA#Ap 
f di=; (38) 
3) x 1-— | N(p), n= p. 
The current in the antenna is given by 
)) 1(€) =2mc(A)n=n0=20€ DY p CpX p(E) Vp(n0) 
} 1 
f &(t)X ,(t)dt (39) 
= 
a => 1ge€ z i 3 X p(§) Y,(n), 
P N(p)Y>, (no) 
where g stands for the constant 
" g=2nc(k/p)*a (40) 
and a=fno is the semi-major axis of the spheroid. 
Under the assumption that the set of functions {(1—£#)-1X,(£)} is complete, we may write 
h | t ] 
, ! 
y ; f S(I)X,(t)dt 
S | X ,(é) Pas, X p(&) (41) 
&(&) = 2 hp —_= 2 — 
: p — £€ Pp N(p) 1— 
1 
| The meaning of Eqs. (33), (39), and (41), which formally solve the antenna problem, may be clarified 


by considering the situation which would exist if there were acting alone the “‘unit’’ driving field 
| & p(t) = X ,(£)/(1—&?). (41a) 
The current in the antenna would then be 


T,(€) =igeX »(E) Vp(no)/ Vp’ (no). (39a) 





Just the pth mode of oscillation would be present—none of the others, whose coefficients all vanish 

by the orthogonality property (38). This result holds regardless of frequency, whether that be right 

for resonance of the pth harmonic or not; but of course the required method of excitation (41a) 
\ varies with the frequency. By (41), then, the driving field is broken up into ‘‘segments”’ each of which 
excites only a single mode in the current. Since the behavior of the spheroid is different for each 
harmonic, properties of the antenna will depend upon the method of excitation : not ordinarily very 
critically for resonant antennas, however, for as shown below the resonant mode usually accounts 
for the major part of the phenomena. 
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Having thus separated the reaction of the antenna into parts, we may proceed to study the be- 
havior of the single mode given by (41a) and (39a). 


1. Properties of the pth Mode 


One of the first things now to be done is to justify the use frequently made in the foregoing of the 
term “‘resonance.’’ The resonant behavior near e= }p7 is due to the factor Y,’(0) in the denominator 
of the expression for the current (39a). By (22a) it is seen that Y,’(0) is of the form 


= not1 
V, (no) =6 log +c, (22b) 
no— 1 


where } and ¢ are not very different in magnitude, though one is real and the other imaginary. Since 
the spheroid is of great eccentricity, the factor log (no+1)/(no—1) is quite large; and 3} is a function 
of frequency which goes through zero at «= 3pm. (22b), then, is a typical ‘‘resonance denominator” 
consisting of two terms, the larger of which vanishes at the point of resonance, leaving the other as 
a small residual so that the fraction is not infinite. Accordingly the current in the pth mode (39a) 
is much greater when ¢= $7 than for neighboring values. 

The sharpness of this resonance is greater for thinner spheroids, whose eccentricity 1/0 is nearer 
unity. A first approximation to the band-width of the antenna is readily computed, and is tabu- 
lated below. 


The phase angle ¢ by which the current leads the driving field can be found in the following 
manner: in the fraction 


T,(&) V,(no) 
= ige(1—é*) 
&p(&) Vy (no) 
(—1)rt'y, 
e'*+ ; [ —G,(1) } (42) 
1€€» 
= ige(1—£*) : 
(—1)Ptly, no+1 
1ee'* + E —1-—G,'(1 | 
1€€ » no— 1 


it is known from (41) that &, is in phase with the driving field if that is all of one phase throughout 
the length of the spheroid. Thus the phase of the current is t¥e same as that of the fraction (42). 
From the qualitative discussion it is known that Y,(), which is Y,(1) in the approximation used 
here, has always the same phase angle as exp (ipa/2) ; hence changes in the phase angle of the current 
depend only on the ratio of the real and imaginary parts of Y,’(mo). Writing 


VY,’ (no) =1"* "(A +7B). (43) 
Then 


tan @=B A. (44) 


Some more properties of ¢ follow from the qualitative discussion : to begin with, A has always the 
same sign, for, considered as a function of , it is the analytic extension of X ,’(n), with a constant of 
proportionality which is real and, although a function of frequency, cannot change sign since it is 
never zero or infinite for finite e«. (Continuity is of course assumed.) Thus A is proportional to X ,’(1) ; 
since the latter is always negative, A cannot change sign. 

A direct consequence is that A is always positive, for at resonance 


A )eapr/2=t-? | + 1€e"* jeep rj2= pa/2>0. (45) 
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Explicitly, the expression for the phase angle is 


Yp no+1 _ oo dca 
| log ————1— ®e{G,'(1)} |+9m[7-"e** ] | 
1 € Cn no—1 | 
g=tan—!"*+ —— — — = =f. (44a) 


Yp 
Rel 1-"ee** |+ gm[.G,'(1) | | 


€\ ep 


For thin spheroids the logarithmic term dominates all the rest, at any rate for the lower harmonics; 
and the fraction has the same sign as y,. At the resonant frequency, y,=0 and the antenna, in this 
approximation, presents a pure resistance load to the driving field. Below the resonant frequency the 
mode is of capacitative character, while above resonance the reactance is inductive. The sharper the 
spheroid, the larger log (no +1)/(mo—1), the faster ¢ changes with frequency, and the sharper the 
resonance. 

The radiation field due to the pth harmonic is given by (4), (5), and (6). At large distances E, is 
: and J/,; the two latter constitute the radiation field: 


s 


negligible compared to E 


€7 0 1 , 
i= — X ,(é) exp [i(en—wt—3pr—¢)], 7 large, (46) 
V,'(mo)| {1 — 2) (yn? -—€2) 3 
K\? Eno 1 7 
1.=( ) X,(é) exp [i(en—wt—3pr—)], 1 large. (47) 
Mm V,'(no)! | (1 —€*)(n?—-1) 33 


In these two equations the asymptotic values Y,(y)~e"” and Y,'(n) ~iee’*” have been used, as 
well as (43) and (44) for Y,’(m). 
The Poynting flux —cE:-//, in the direction of increasing 7 is then at large distances 


ce’no?(k/p)? X »*(—) cos? (en—wt — 3p — ¢) 


—ckh:ll,= , 1 large, (48) 


9 


Vp’ (m0) 2 1—¢ | (mn? —&*)(m?—1)}? 
while the element of area of the large spheroid is 
dS = 2xf?} (n? — €*) (n? — 1) | dé, (49) 


whence comes the mean rate of radiation of power 


* 


e-ny” K\? @! X,7(&)dé 
— bef E: 11 dS = tcf? ( ) | 
Ss YV, (no) : Me 1 1 — §* 
dU K\? €a" 
-| | =xe( ) N(p). 
dt Spy m Vy’ (no)? 


The radiation resistance is defined as the pure resistance which would consume the same amount 
of power as the spheroid when carrying the same current. Since the current varies from point to 
point, the radiation resistance will depend upon the point chosen for reference—usually a current 
loop. Except at resonance, the loops vary in size and the choice most convenient may vary with the 
method of driving the antenna; accordingly the radiation resistance at a point £ is defined as a func- 


[dU /dt \n u\? N(p) 
Rraa= = ( ) , ($1) 
$|I|? KJ 2mwcX »*(E) | Vp(no) ,* 


(50) 


tion of that point : 


¢ 
we 
w 
Jt 
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Combination of (51) and (44a) gives the reactance due to radiation: 
X rad = — Rraa tan ¢. (52) 


It is more usual to give the field strength as a function of antenna current rather than in the form 
(46). Combining with (39a), 


. no 1 I,(é) . 
E:=-—- gilen—wt—}pr) (46a) 


9 1 


g Y (no) (n? — &)? (1 —£*)? 


or, in practical units at large distances, 


u\? 1 59.96 I,(&) 
E;= -( ) ; eee) volts/cm. (46b) 
K Y »( no) r (1 — ¢*)! 


This last equation gives a relation between the current distribution in the antenna and the direc- 
tional characteristics of the field. In polar coordinates, if the field is entirely due to the pth harmonic, 


F4(8) 1 |J,(&) 
=— ‘ (53) 

Ee(0)| sin @\7,(0) 

Unfortunately, if other harmonics are present in appreciable amounts, the simple relation (53) does 

not hold for the total current on account of the weighting factor 1/| Y,(mo), in (46b). 


2. Properties of the pth Harmonic near Resonance 


Formally the solution of the antenna problem is now complete, but the results will be of quantita- 
tive value only if the functions X and Y in (26) and (21) can be computed sufficiently accurately. 
The complexity of some of the preceding expressions is due to the fact that they hold for all fre- 
quencies, and as the qualitative argument shows the functions change very considerably in character 
over the infinite frequency range. Since in practice the harmonics are important near resonance, the 
simple approximate formulas due to the leading (sinusoidal) term only may be used. First, at 


‘ resonance: 
€=3pr; yp=0; X,=sin $pr(1—£); V,=exp (t3prn); (54) 
. K ; 1 r 
cod p= if ) €(l —X p(&) Y,(n). 
M - (no) 
Length 
1 1 
2a =2fno=p 


2 (xu)? }1—(b/a)?}? 
The current in the antenna due to the pth harmonic is 


; ; Y »( no) : pr ~- 
I ,(&) =igeX »(£)— =g sin —-(1—&) (95) 
7» (no) 2 


in phase with &,, and is sinusoidal. In practical units—& measured in volts per centimeter and the 
current in amperes, 
I ,(é) =0.01668(x/u)*a sin 5pr(1—&). (55a) 


At resonance, the field at all distances is given by the sinusoidal functions. In Heaviside-Lorentz 
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52) 


rm 


da) 
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units, 


sin pm(1—£) 

_ — _ eg tore-os bpm) | (56) 
1 (1 — &?) (? — &*) }? 
cos $pr(1—&) 

| (mn? — 1) (mn? — &*) }3 


k\? sin }pr(1—&) 
I= ) ae Psi Rie) ei(hprn wt—kpr) (58) 
u/s \(1—£?)(n?—1)}} 


E;= =a 


E,= — no ei(aprn wt—hpti) | (57) 


or in volts per centimeter 
sin 3pr(1—&) 

ere eens se'torr-we-toe) (56a) 

i(i—F)(9— &)7* 

cos $pmr(1—é) 


{(n?—1)(n?—#)}} 


E;= — ie 


E,=—n eidrrrwt Kot Dn, (57a) 


It seems worth noting that for large distances (56a) becomes the approximate formula in standard 
use for radiation from a free wire in space.*°*! If we change to spherical coordinates, fn=r and 


a sin [}p2(1—cos @) | prr pr 
E,= —- — exp i( —wt—— *)| (56b) 
r sin 6 re 2 


&=cos @: 


The formulation of the present paper has the advantage that a quantitative connection with the 
driving field (41) as well as with the current (39a) is found. 
The radiation resistance for the resonant harmonic is 


u\?} N(p) 
Read = - ° (59) 
kJ 2nc sin® $pr(1—é) 


If the radiation resistance at resonance is referred to a current loop, then X ,?(£) =1, and the radiation 
resistance in ohms is then 
R= (u/«)'29.98F(2p) ohms, (59p) 


*»™ 1—cosx 
F(2p)= { dx. 
x 


wherein 


Since at resonance the current is in phase with the driving field, the phase angle and reactance both 
vanish. 

If the pth harmonic is not precisely at resonance, a complete and exact discussion would require 
the use of all the machinery developed previously, the perturbation expressions (26), (27), (21a), etc. 
It is, however, possible near resonance to get a reasonable approximation without much labor, simply 
by the use of only the first few terms. 

Near resonance, the frequency parameter becomes 


2avf (Ku)! pr? i pr 9, 
mew ) +6,| = 4+—+..., (60) 
C 2 2 pr 
We still take 
X ,(&) =sin 9pr(1—&), (61) 
20 Carter, Hansell, and Lindenblad, Proc. I. R. E. 19, 1773 (1931). 
2F, E. Terman, Radio Engineering, second edition (1937), p. 686. 
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thus assuming sinusoidal current distribution; but the parameter y, is no longer 0; it is 


Yo=Ypi19p= —_ “i 6... (62) 
F(2p) 
When we continue the approximation that all the a’s are zero except a,,=1, there result: 


Si(2pr) 6, 
Y »(no) = y,(1) =i 1— | 


7 (63) 
| F(2p) pr| 
7” (Pm 0, 6, Si(2pr) 6, not 1 | 
VY, (no) =1"! + --— —1 log (64) 
| 2 pr 2 F(2p) F(2p) 


se—1) 


Immediately from (44) and (64) comes the phase angle by which the current leads the driving field: 


0, no+1 pr 1 1 Si(2pr) 7 
tan ¢= — log /| +0,( — )I (65) 
F(2p) no—1 2 br 2 F(2p) 
The radiation resistance (51) is 
u\? F(2p) 
Read -( ) (66) 
K 4rc Y ,(no) . 
or in ohms 
u\? 29.98F(2p) 
R,* ) | 
K Si(2pr) 0, (67) 
1—?2 


F(2p) pr 


when referred to a current maximum. The field follows the current both in magnitude and phase 
(46a). The current (39a) in amperes is 


Si(2pr) 6, 





= 
F(2p) pr 
I ,(£) =0.01668(« /u)*ae sin $pr(1—€&)- ' (68) 
pr 1 1 S1(2pr7) 6, not1 
+6, — —1 log 
2 pr 2 F(2p) F(2p) no— 1 


Since the logarithmic term in the denominator of (68) is large if the frequency is not too close to 
resonance, with increasing slenderness of the antenna y»—1 becomes smaller and the resonance 
sharper. 
Since this current has been evoked by the “unit” driving field (41a), the power consumed takes 
the form 
l 


P= | & »(&)1,(&) -fnodé 


1 
Si(2pr) 
. _ 6, (69) 
K\? pr (2p) 
= 0.01668 ) a°>N(p)R 
u Si(2pr) Op no+1 
ji- 6,—1 log 
2eF(2p) eF (2p) no— 1) 


The frequency-dependent portion is confined to the expression whose real part is to be taken. 
If now in the small correction term involving @, in the numerator we replace px by 2e, considerable 
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62) 


53) 


14) 


~ 
— 


se 








simplification results: 


P=1yaN(p) KR) , ; 
no+1 
6, log (70) 
’ no— 1 
| 1—1 is 
eF(2p) — 30,Si(2pr) 


wherein Jy is the current at a current loop at the resonant frequency. 

It will be noted that even in this first approximation, on account of the correction term in Si(2p7) 
the pass band of the antenna is not quite exactly symmetrical about the resonant frequency. Numeri- 
cally, Table I shows the asymmetry to be quite small in the case of the first mode. The edges of the 


TABLE I. Properties of the first harmonic near resonance. 


log ™*! =10 log M+! a5 log TT ~ 20 log +! 240 
no—1 no —1 ~ no) no —1 
Rp () X ¢ X ¢ X ¢ X 
a1 v/i ohms degrees ohms Io/lor degrees ohms Iy/lo degrees ohms To/lo degrees ohms To/Tor 
—1.2 0.75 51 72.6 160 0.28 
—1.0 0.80 53 69.4 140 0.34 
—0.8 0.84 56 64.7 119 0.42 72.5 180 0.29 
—0.6 0.88 60 57.7 95 0.53 67.2 142 0.38 72.5 190 0.30 
—0).4 0.919 63.7 46.4 —67.0 0.684 57.6 101 0.532 64.6 -—134 0.426 76.6 —268 0.230 
—0.2 0.959 68.04 27.7 —35.66 0.884 38.2 53.5 0.784 46.4 —71.3 0.688 64.5 —142.7 0.429 
—0.1 0.980 70.47 14.7 —18.44 0.966 21.4 —27.65 0.930 27.6 -36.87 0.885 46.3 —73.74 0.690 
—0.05 0.990 71.75 11.1 —14.07 0.981 14.6 —18.76 0.967 27.6 —37.51 0.885 
—0.025 0.995 72.41 14.6 —18.92 0.967 
0 1 73.08 0 0 1.000 0 0 1.000 0 0 1.000 0 0 1.000 
0.025 1.005 73.76 —14.6 19.26 0.968 
0.05 1.010 74.46 —11.1 14.57 0.980 —14.6 19.43 0.967 —27.6 38.86 0.885 
0.1 1.020 75.89 14.7 19.79 0.967 21.4 29.68 0.930 —27.5 39.57 0.885 —46.2 79.14 0.690 
0.2 1.041 78.93 27.5 41.08 0.885 —38.0 61.61 0.787 —46.2 82.16 0.692 —64.4 164.3 0.432 
0.4 1.081 85.8 16.0 88.9 0.689 57.3 133 0.538 —64.3 178 0.432 —76.4 356 0.233 
0.6 1.12 94 §7.2 146 0.53 —66.7 219 0.39 —72.1 292 0.30 
0.8 1.16 104 —64.1 210 0.43 72.1 320 0.30 
1.0 1.20 116 68.7 300 0.35 
1.2 1.25 130 —72.0 400 0.29 


pass band, taken arbitrarily at the half-power points, give a band width in terms of 6, 


no+1 Si(2 pr) 
AO, =2€F (2p) log — , 7 
ny — | not (71) 
4 log 
no— 1 


In terms of frequency, setting v, equal to the resonant frequency gives the simple approximate band 
width 

Av / pr not 17 

a= fF 2p) / 


V, Z 


(72) 


log ’ 
no 1. 


in which the correction term in Si(2pr) has been omitted. The presence of higher harmonics will 
slightly widen the band. 


V. COMPOSITE OSCILLATIONS 


Up to this point has been discussed only the 
behavior of the pth harmonic (or mode of 
vibration) alone; but if any field other than a 
multiple of the special one &,(&) = (1— &)-“XY,(&) 
be acting, then harmonics other than the pth 
will also be excited, by (39). As an example 
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consider a vertical quarter-wave antenna excited 
at the base. If the ground be regarded as per- 
fectly conducting, the image antenna will be 
similarly excited, and the two together may be 
regarded approximately as a single centrally- 
driven prolate spheroid. If we assume the 
original driving e.m.f. to have been Ve~‘*! volts, 
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(41) may be applied: 
&(&)=2V/fmodAt, —AE/2<E<AE/2; 


&(£)=0 elsewhere, (73) 


1 f' 1 2V 
h,=— f &(t)X ,(t)dt= -X (0). (74) 
N(p)4_1 N(p) a 


When this equation is specialized to the funda- 
mental frequency, we find e=$2, X,=sin $2 


xX (1—£), and X,(0)=1. Thus 


h,=1.6409 V/a. ( 


~~ 
wm 
— 


Rather involved calculations give for the third 
mode 


h3= —1.075V, a. (76) 
Higher modes are neglected in this example; the 


same method is applicable, but they are small. 
The current at the center of the spheroid is 


uy! 
(“) I(0) =0.02738V 
K 


— 0.00545 V 


; 


: . (77) 
log (no+1)/(n0—1) 

As expected, the current due to the first mode is 
much the major part of the total; and the 
contribution of the third harmonic is practically 
90° out of phase with the driving field. 

It is apparent that the phase angle of the 
total current is not zero when the first harmonic 
alone is in resonance. If the usual engineering 
criterion of zero phase angle be taken as the 
criterion for since the third har- 
monic is capacitative in character, the frequency 
must be slightly increased so that an inductive 
component of the first mode will appear to 
cancel this capacitative current. (It should be 
noted that this cancellation will not hold 
throughout the spheroid, because the current 
distributions of the two modes differ; accord- 
ingly, the ‘‘resonant”’ frequency depends slightly 
on the driving point chosen.) For the center- 
driven antenna, approximate calculations based 
on (65) and 


“resonance,” 


1,(0) =0.01668(«/u)'a cos ge~‘¥ 
(68a) 


=I) cos ge’, 
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which holds for frequencies near the resonant 
point, show that this cancellation takes place 
when 


—0.2 
tan g= ’ 
log (no+1)/(no—1) 


€ 0.16 
=1+ ; 
31 [log (mo+1)/(me—1) }? 


(78) 


For even the fattest of the antennas considered 
in Table I, this frequency shift is only 0.16 
percent. 

Another meaning of “‘resonance”’ differs from 
both the foregoing: it is maximum total antenna 
current. If we continue to neglect the small 
variation of the non-resonant third harmonic 
with frequency, it is clear that maximum total 
current will occur quite near the resonance of its 
principal component at e=}pz; in fact, the 
frequency for current resonance will be just 
slightly lower than this. The frequency shift for 
this case is —0.15 percent for the fattest antenna 
of Table I. Here again the “resonant”’ frequency 
depends upon the point chosen. Page and 
Adams* have computed this type of resonance 
for various antennas excited by a field whose 
complex magnitude is 
length. 


constant along their 
In practice, optimum antenna performance 
would seem to demand not zero phase angle, 
nor yet maximum current; but instead the 
criterion might well be maximum rate of power 
radiation. In the general case, if the amounts of 
harmonics present are given by the numbers 


1 1 
h,= f sox,(oai 
N(p) 1 


[see (41) |, then the mean rate of radiation of 
power is 


dU K\? 
| = re( ) a*e? >> 
dt Any m P 


obtained by summing (50); cross terms between 
harmonics drop out by orthogonality (38). If 
as heretofore the frequency dependence of the 
non-resonant harmonics is neglected, then maxi- 
mum rate of radiation of power comes at 
«= pz, so that the sense in which we have been 


hy?N(p) . 
in h.l.u. (79) 
Y,'( no) | ’ 


JOURNAL OF APPLIED PHYSICS 





nt 
ce 





ee 


“resonance” for each mode 
applies also to the antenna as a whole. 

To return to the centrally excited half-wave 
antenna, its radiation resistance Re;; is found: 


using the term 


Ry? + R3\ [3|? = Rers| I) ?, (80) 


0.04 
Ress = 73.081 —_-— 
[log (no+1), (no 1) }* 


m b 
x( ) ohms. (81) 
K 


Thus the radiation resistance of the antenna is 
only slightly lowered by the presence of the 
higher despite the smallness of the 
radiation resistance of the latter, which is on 
the order of 0.3 ohm. 

The radiation field at 


mode, 


large distances is 
practically the same as it would be if the third 
harmonic were absent. This at first sight amazing 
result holds in spite of the fact shown by (77) 
that the amplitude of the third mode may be 
percent of the first. Mathematical 
explanation is given by Eq. (46b), wherein the 
proportionality factor | Y3(mo)| is some 17 times 
larger than the corresponding | ¥i(m0)|=1 at 
this frequency. Although the third harmonic 


several 


current may be, say, 5 percent of the first, its 
radiation field is a mere 0.3 percent. No violation 
of physical principle is involved, as can be seen 
qualitatively from the character of the current 
distribution : the two current loops of third mode 
at the ends of the antenna are out of phase with 
the larger central loop, and cancel the field at a 
distance by interference. Alternatively, one may 
say that the explanation is the smallness of the 
third mode’s radiation resistance. 

The foregoing furnishes an explanation of the 
rather surprising success attained in predicting 
the field of an antenna by the oft-used means of 
assuming sinusoidal current distribution. For 
example, Pedersen” has shown that this assump- 
tion gives quite good results for the Copenhagen 
aerial he studied, in spite of the fact that the 
current has been measured and found not to be 
sinusoidally distributed. Physically it is, of 
course, evident that the current in the antenna 
must be the cause of the field; yet the contention 


2 P.O. Pedersen, Ingen. Vidensk. Skr. A. No. 38 (1935). 
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is sometimes advanced that there are ‘‘wattless”’ 
components in the current, in phase quadrature 
with the driving electromotive force, which do 
not contribute to the field. Such a situation is 
exhibited here: the third harmonic leads the 
impressed field by nearly 90°, and also contrib- 
utes practically nothing to the radiation field at 
large distances. It is therefore evident that an 
attempt to predict the radiation field from the 
measured root-mean-square current distribution 
of an antenna may encounter serious difficulty 
—may easily, in fact, go farther astray than a 
bald assumption of sinusoidal current. 


1. Composite Oscillations in General 


In principle, oscillations induced in a thin 
spheroid by a field harmonic in time but of 
arbitrary frequency and spatial distribution may 
be treated by this method. All that is necessary 
is to compute the values of X,(&), Yp(1), Yo’ (no) 
for the relevant harmonics; by (41) break up 
the driving field into corresponding segments; 
whence the resulting numbers hf, give the 
amounts of the various harmonics present, in 
terms of the standard “unit’’ driving field &, 
(47a). In practice, the method is particularly 
suited to resonant antennas where X ,, essentially 
the current distribution, may, for the most 
important harmonic, be approximated by the 
first term of its Fourier series. The analysis of 
paragraph 5.2 can be extended to include the 
second power of 6, but progress beyond that 
point runs into considerable computational 
difficulty which impedes application to non- 
resonant structures. Some entirely different 
approach such as that of S. A. Schelkunoff* 
would be more suitable for these last. The 
present theory may be regarded as a justification 
and expansion of the sinusoidal current assump- 
tion for resonant antennas; particularly for the 
radiation field this assumption is better than 
might, at first, be expected. Our principal 
extension is to relate the amplitude and phase 
of the current to the driving field (39), (44), 
(68a). 

In concluding, a few interesting minor points 
might be mentioned. The first is pertinent to a 
higher 


discussion of the resonances of the 


23S. A. Schelkunoff, Proc. I. R. E. 29, 493 (1941). 
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antenna: the current nodes predicted by the 
simple sinusoidal theory are not found experi- 
mentally but instead are found nearby minima 
in the measured (root-mean-square) current. 
Often the fact that these minima are not zero is 
attributed to losses of a true ohmic resistive 
character, but even for a perfectly conducting 
spheroid the minima do not vanish in general, 
as the small contribution of the non-resonant 
modes remains as a residual. 

Another point is the following. If a vertical 
antenna is three-quarters of a wave-length long, 
it is at the proper length to resonate the third 
harmonic. In practice, such resonance is quite 
easy to establish by the action of a lumped 
electromotive force; but under the influence of 
an incident plane wave acting along the whole 
length of the (receiving) antenna, the resulting 
current distribution is not markedly sinusoidal, 
but seems fairly smoothly distributed with no 
great tendency to show nodes. 

Qualitatively the result may be explained by 
considering the coefficients /, involved in the 
second case. 


1 ol 8S 1 
h,= | &-X,(t)dé~ , 
N(1) « 1 w 2.4 
2 : 3m 8S 1 
h;= J & sin —(1—é)dé= : 
F(6) l 2 wr 10.5 


The evaluation of /; is rough, but indicates that 
the field has something like five times as great a 
tendency to excite the first harmonic as the third. 
The presence of a large amount of first harmonic 
current in phase quadrature could account for 
the concealment of the nodes of the third. For 
the lumped field, the 4, are more nearly equal, 
and the resonant behavior shows up. 

To Leigh and 
Norman I. Adams the author wishes to offer his 
thanks for illuminating discussions. 


Professor Page Professor 


APPENDIX I. F(2p) AND Si(2pzr) 


These functions, used frequently in the work, seem not 
to be available readily without interpolation. The follow- 
ing values were interpolated from the British Association 
Mathematical Tables. 

F(2p) = = 1 —— t j 

Jo f 


=log 2px — Ci(2pr) +0.57722, 


dt 
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Si(2pr) — sin f it 

SU Lpr) = . i dt. (A1-2) 
p F(2p) Si(2px) | p F(2p) Si(2pr) 
1 2.43766 1.41815 6 4.20756 1.54428 
2 3.11435 1.49216 7 4.36153 1.54807 
3 3.51648 1.51803 8 4.49493 1.55091 
} 3.80296 1.53113 9 4.61264 1.55312 
5 4.02553 1.53903 10 4.71794 1.55488 


Higher arguments are found to this accuracy from the 
asymptotic series 


F(2p) =log (2pm) +0.57722+4 : sarees, (ALG) 
pr 
' w | 2! 
niet at r= eee (A\1-2a) 


APPENDIX II 


Radiation resistance at resonance for the pth mode is 
given by 


Read = (u/ «)?29.98 F(2p) ohms. (A2-1) 
p R»(p) (ohms) p R»p(p) (ohms) 
1 73.08 6 126.1 
2 93.37 7 130.8 
3 105.4 8 134.8 
} 114.0 9 138.3 
5 120.7 10 141.4 


It should be emphasized that the tabulated values above 
are not the radiation resistance of the antenna at the corre- 
sponding frequencies; though for thin spheroids, where the 
resonant harmonic dominates, this is a good first approxi- 
mation. The exact radiation resistance of the antenna of 
course depends on the proportions of harmonics present. 

\gain, an ingratiatingly smooth curve can be drawn 
through these points, but between the points the curve 
does not correspond to any harmonic, and only by chance 
would it correspond to some particular method of driving 
the antenna. The radiation resistance of each harmonic 
varies with frequency, but along curves quite different 
from this; see, for example, Table I. 


APPENDIX III 


Properties of the first harmonic near resonance. The 
values of Table I were calculated for four different eccen- 
tricities: 

(1) log (no+1)/(mo—1)=10; b/a=1.35 (10) 
Band width Av/v,=0.155 from (78) 


(2) b/a=1.11 (10) 


2 log (no+1)/(mo—1) = 15; 
Band width Av/v,-=0.103 
(3) log (no+1)/(no— 1) = 20; 
Band width Av/v,=0.078 
(4) log (no+1)/(nmo—1) =40; 
Band width Av/v,=0.039. 


h a=9.08 (10 ° 


b/a=4.12 (10)~° 


Here v, is the resonant frequency, R, radiation resis- 
tance, ¢ the phase angle by which current leads the applied 
field, X the driving point reactance, referred (like R,) toa 
current loop; and J» is the current in the center of the 
antenna, having the value Jo, at resonance. It is assumed 
that a “unit’” amount (47a) of first mode appears. 


The data of Table I are calculated from (65), (67), (68), 
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7) 
}28 
07 
91 
}12 


he 


a) 





and (72), which become for this case 


1 
— 0.2616, log or 
a jo : 
tan ¢= 1 +0.0170, 2 (65-1) 


73.08 


Ro= 70.3700,’ 


(67-1) 
x\} in 

0.01668( ) acos>& 
Mu 


=> - _ 8-1 
 _; 0.26101) mol" niles 


~*7+0.0170, > a —1 
0.2610, ve+1|- 


~ 
Str 


0 To, == 1- = pty ad 
(Lo/Tor) = |1— + T0170, °F e—1| ° (OF!) 
Ap 1.55 
= - 72-1 
Vr | no+ 1 \/ 
og —— 
; no—1 


The sign of 7 is reversed from the usual convention. 
APPENDIX IV 


Numerical constants for the first three odd harmonics. 


rABLe II 
k (2/4) ayu (974/16) = a324 (16/64) = Aaa 
1 0.024055 —0.00119 7 0.00014 33 
3 
5 0.02426 6 0.00018 53 
7 -0.00647 81 0.00031 97 
9 0.00272 95 0.00016 88 
11 0.00141 82 0.00009 543 
13 0.00083 427 0.00005 817 
k (2/4) ask (94/16) = a2 
1 0.00416 22 0.00001 7063 
3 0.02119 8 —0.00041 624 
7 —0.01831 9 0.00001 2531 
9 0.00553 24 0.00014 396 
11 0.00253 19 0.00008 7392 
13 —0.00139 68 0.00005 3792 
15 —0.00086 000 0.00003 5051 
17 —0.00056 989 0.00002 4025 
19 0.00039 824 0.00001 7158 
21 —0.00028 977 0.00001 2657 
k 11k 12 
3 —0.01406 3 0.00039 50 
5 0.00278 57 0.00017 80 
7 —0.00099 229 0.00007 3 
9 —0.00046 255 
h y 1k ‘Y 3k : 
1 0.82046 —0.56875 —0.49683 
x? = x 
2 0.00186 92 0.004281 ( ) -0.00173 04 { : ) 
ni n* 
3 0.00003 0.000037 ( — ) 0.00001 898 ( ~ ) 
16 16 


In computing the properties of the third mode which 
were given in paragraph 6 above, use was made of a nu- 
merical check. Since the solutions at zero frequency are 
known,® the series for the a’s can be checked against the 
exact values for «e=0. They were found to correspond 
within a few percent. Data in Table III were then com- 
puted for the third mode at e=72/2, the resonant point of 
the first mode. 


APPENDIX V. DRIVING POINT IMPEDANCE 


By the method of paragraph V can be found a formula 
for driving-point impedance to replace (52) and gen- 
eralize (51): 
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TABLE III. Properties of the third mode at the resonant point 
of the first. 


a3 = —0.33247 
433 =1 

435 =(0).20609 
a3i =0.07458 
439 =0.03473 
43,11 =0.01898 
43,13 =(0.01148 
73 =10.70 
e3X3/(1) =0.6649 
X3(0) = —1.179 
N(3) = 2.195 

“Ses 

— =104.3i 
17e3 

Y3(1) = —17.12i 


; o+1 
Y¥3'(mo) =104.3i log ; +0.0918 
0 


0.00507: ( : )a 


13(0) = . — for the “‘unit’’ driving field. 
log bl 
no —1 
For the centrally-driven antenna of paragraph 6, 
hs = —1.075V/a ; 
(-7¥ 
1:0) = —0.00545i —=+ 
0 + 1 
log 
i ” es 1 
Rrad 0.323 ( e ) ohms. 
« A 
, (uy no+1 
Xrad = —336.9i ( ) log ohms. 
K no —1 
»+1 
tang =1140 log ” 
no —1 


Zr id = Rrad -= 1X vad 
- (“) 59.960N(p) Y,'(no) 


= ae hms. (51G 
X »?(&) le } p(No) _ \ 


kK 
Here as throughout the sign of the imaginary unit i=(—1)* 
is opposite to the usual convention; R and X are the or- 
dinary reactance and resistance which the driving point 
e.m.f. sees. Using the approximations (63) and (64), the 
impedance presented by the near-resonant harmonic driven 
at a current loop is 


} 
Z, <2, ~iX,« (“) 59.96 F(2p) 
K 


fl _1Si2pr)_. % 1, wt!) 
4 _ ~ —i{— og > 
lpr 2 F(2p) F(2p) * no—1J 
},  Si(2pr) 0, \ 


\” Fp) pwd 


Pr +O, 
, ee ohms. (67G) 
€ 





Letter to the Editor 








Effect of Machining and Aging on the 
Damping Capacity of Steel 
LESTER TARNOPOL, Raytheon Production Corporation, 
Newton, Massachusetts 
AND 


JAMES R. MorGAn, Alabama Dry Dock & Shipbuilding Company, 
Mobile, Alabama 


February 20, 1942 


EVERAL years ago, Sonnemann! found in one instance 
that two bars from the same consignment of standard 
mild steel exhibited widely different damping capacities. 
The damping capacity of one bar was found to possess from 
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8 to 10 times that of the other bar at high values of strain. 
Che specific damping capacity is the percent of the total 


energy absorbed per cycle during free vibration in these 
measurements, 

In order to investigate this phenomenon, we machined 
four torsional damping capacity specimens from com- 
mercial cold rolled S.A.E. 1020 steel. These specimens, in 
the form of solid cylinders with square center and ends, 
were tested in a Cambridge Torsional Damping Recorder. 
Considerable spread in the specific damping capacities at 
high values of maximum fiber stress was found. Whereas, 
at 6000 Ib. per sq. in. (maximum fiber stress) all specimens 
had a specific damping capacity of 2 percent; at 24,000 Ib. 
per sq. in. the specific damping capacities were 7.6, 9, 14.5, 
and 27 percent, respectively. 

Since these four specimens varied in diameter from 
0.3109” to 0.3124”, the need for greater precision was felt. 
In order to eliminate any effect due to differences in the 
dimensions of the specimens, three precision specimens 
were machined from a single bar of cold rolled S.A.E. 1020 
stock steel as shown in Table I. 


Tasie I. 
Specimen No. Average diameter inches 
5 0.3041 
6 0.3041 
7 0.3041 
rasce Il. 


Rockwell hardness 


after 4 hr. at 400°F 
furnace cooled 


Original Rockwell 
Specimen No. hardness 


B93 B93.5 


1 

2 B93 B95 

3 B93 B94 

4 B94 B96 

5 B101.5 B102.5 
6 B101 B102 

7 B103 B103.5 


The three precision machined specimens likewise had 
specific damping capacities of 2 percent at 6000 Ib. per 
sq. in., whereas at 24,000 Ib. per sq. in. the specific damping 
capacities were 4.8, 8.1 and 18 percent, respectively. How- 
ever, after aging at room temperature for several weeks the 
damping capacities of these specimens were found to have 
decreased considerably, bringing them all very close to- 
gether. At 6000 Ib. per sq. in. the specific damping capaci- 
ties became 1 percent while at 24,000 lb. per sq. in. the 
percentages had dropped to 3, 3, and 4. 

All seven specimens were subsequently given a stress 
relieving anneal for four hours at 400°F and furnace cooled. 
After this low temperature anneal, the damping capacities 
of all seven S.A.E. 1020 specimens were very nearly the 
same; i.e., 1 percent at 6000 Ib. per sq. in. and 3 to 4 per- 
cent at 24,000 Ib. per sq. in. After annealing, Rockwell 
hardness tests were made which showed very little change 
from the original values (Table II). 

A similar phenomenon was observed with S.A.E. 2320 
steels containing 3.3 percent nickel, the main difference 
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being that the nickel steels did not reveal such large differ- 
ences in damping capacity among themselves. 

Since the fiber stress increases linearly from zero at the 
cylinder axis of the specimen to a maximum at the outside 
fiber, the damping likewise to a first approximation may 
be considered to increase linearly from zero at the center 
to a maximum at the outside fiber. Therefore the condition 
of the outer volume of fibers would have the most effect 
upon the damping capacity. Since one of the precision 
specimens showed a very high damping capacity immedi- 
ately after machining and subsequently lost its high value 
after aging, it would appear that the high damping ca- 
pacity may have been due to machining which leaves 
stresses in the outside fibers. The original variations in the 
damping capacities of specimens No. 1, No. 2, No. 3, and 
No. 4 might also be partially explained by machining. At 
any rate, the high values of damping capacity were very 
probably due to cold work, since the low temperature stress 
relieving anneal reduced the damping and brought all 
specimens nearly into coincidence. 

The damping resulting from cold work might be par- 
tially explained by encouraged plastic flow resulting in 
stress-strain hysteresis losses. Damping also results from 
heat flow energy losses caused by temperature gradients 
produced during deformation, and by eddy current losses 
due to magnetic effects.?:4 

Machining could encourage plastic flow by producing 
lattice warping in the outer fibers, wherein voids of atomic 
dimensions might exist along with stressed material. These 
stresses plus temperature energy could greatly increase the 
energy of certain atoms. Thus a further stress well below 
the elastic limit of the material would be sufficient to 
plastically move one of these atoms into a neighboring void. 
This might start a process whereby a whole series of atoms 
would flow plastically, following each other in a row, result- 
ing in potential energy losses or damping. 

After aging at room temperature plus the effect of vi- 
bration when testing specimens, or during a low tempera- 
ture anneal, the effects of machining gradually disappear. 
This results in a considerable decrease in damping capacity 
but effects no change in Rockwell hardness. The relief of 
stresses at 400°F takes place in a few hours and is easily 
understood, since strained atoms can readily migrate to 
positions of lower potential energy at this temperature. 

The migration necessary to accomplish relief of strains 
which occurs at room temperature is more difficult to 
explain theoretically. Nevertheless, we have abundant evi- 
dence from other sources that atoms in metals may diffuse 
considerably at room temperature.‘ The Boltzmann tem- 
perature distribution permits an atom to have sufficient 
temperature energy, at some time, to hurdle its potential 
barrier and move to a neighboring position provided the 
atom is already in a state of strain and is adjacent to a 
“lattice void.”” The continuance of this process would 
ultimately produce the aging results indicated by our 
experiments. 

1 Sonnemann, Mitt. d. Woehler Inst. Braunschweig, No. 28 (1936) 

2 E. R. Parker, Trans. Am. Soc. Metals 28 (1940). 


3C. Zener, Phys. Rev. 52, 230 (1937). 
4R. Hultgren and L. Tarnopol, Metals Technology (Jan. 1939). 
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